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NEAR-FIELD  SEISMOGRAMS  FROM  A TWO-DIMENSIONAL  PROPAGATING  DISLOCATION 

ABSTRACT 

Joe  Jean  Litehiser,  Jr. 

The  exact  mathematical  formulation  for  displacements  on  the  surface 
of  an  isotropic,  homogeneous,  perfectly  elastic  half-space  due  to  an 
embedded  line  dislocation  source  is  developed  after  the  method  of 
Cagniard  (1939)  as  modified  by  de  Hoop  (I960).  This  formulation  is  then 
used  to  model  the  observed  strong-motion  accelerometer  records  generated 
at  Pacoima  Dam  during  the  San  Fernando  earthquake  of  February  9,  1971. 

A number  of  model  accelerograms  are  generated  using  several  dis- 
location time  histories,  fault  geometries,  and  fault  rupture  propagation 
velocities.  The  dislocation  time  histories  considered  are  the  simple 
finite  ramp  and  its  first  and  second  order  continuous  counterparts  for 
rise  times  of  0.1,  0.25,  and  0.5  seconds.  The  two  fault  geometries  used 
are  planar,  after  Canitez  and  Toksoz  (1 °72),  and  hinged,  after  Alewine 
and  Jordan  (1973).  Finally,  fault  rupture  propagation  velocities  are 
considered  over  the  approximate  range  of  from  2.0  kilometers/second  to 
2.5  kilometers/second. 

Within  the  limitations  of  the  mathematical  idealizations  used  it 
is  found  that  the  faulting  process  of  the  San  Fernando  earthquake  is 
best  represented  by  a dislocation  whose  time  history  at  any  point  on 
the  fault  is  that  of  a second  order  continuous  smoothed  finite  ramp  of 
0.25  second  rise  time  propagating  at  approximately  2.3  kilometers/second 
over  a hinged  fault  surface.  The  average  dislocation  amplitude  as  con- 
strained by  the  observed  Pacoima  Dam  peak  acceleration  in  that  part  of 
the  record  for  which  a model  motion  match  is  successful  is  between  one 
and  two  meters. 
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One  fundamental  difficulty  with  all  models  tried  is  the  failure  to 
match  both  vertical  and  horizontal  observed  motions  simultaneously.  In 
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general,  the  vertical  component  is  more  successfully  duplicated, 
problem  is  not  resolved  by  this  study. 
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Chapter  1 

INTRODUCTION  AND  REVIEW  OF  PREVIOUS  STUDIES 
1.1  Introduction 

The  San  Fernando  earthauake  of  February  9,  1971  was  the  most 
damaging  earthquake  in  the  United  States  since  the  Alaskan  earthquake 
of  1964.  Although  only  a moderate  event  of  magnitude  6.5,  sixty-four 
lives  were  lost  and  property  damage  in  excess  of  five  hundred  fifty 
million  dollars  occurred.  Both  of  these  numbers  could  very  easily  have 
been  much  larger  as  the  near  failure  of  the  lower  Van  Norman  Dam  con- 
vincingly indicates. 

This  earthquake  was  also  important  from  a scientific  viewpoint. 
Occurring  in  a metropolitan  region,  it  provided  a valuable  test  of  the 
adequacy  of  modern  engineering  design  practices.  Many  structures,  most 
notably  Olive  View  Hospital  and  many  freeway  overpasses,  did  not  perform 
up  to  expectations . There  is  obviously  still  much  to  learn  in  the  area 
of  earthquake-resistant  design.  Seismological ly , this  earthquake  pre- 
sents many  unique  opportunities.  The  faulting  episode  with  which  it  was 
associated  had  a large  thrust  component.  This  is  unusual  in  California 
where  the  vast  majority  of  seismic  activity  is  related  to  strike-slip 
movement  along  the  San  Andreas  fault  and  its  branches.  Most  important, 
however,  is  the  vast  number  of  high-quality,  strong-motion  records  ob- 
tained from  this  single  event--more,  indeed,  than  in  all  previous  earth- 
quakes combined. 

One  record  in  particular  of  these,  that  generated  at  Pacoima  Dam, 
is  especially  significant  because  it  presents  a recording  of  strong 
ground  motion  at  a site  directly  above  the  region  of  faulting.  The  en- 
deavor here  is  to  seek  a description  in  some  detail  of  the  faulting 
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process  associated  with  the  San  Fernando  earthquake  of  February  9,  1971 
as  evidenced  by  the  strong-motion  record  written  at  Pacoima  Dam. 

This  is  an  inverse  problem  of  some  complexity  which  will  depend 
for  success  upon  the  applicability  of  the  medium  transfer  function  used 
as  well  as  upon  the  judicious  use  of  the  work  of  many  previous  investi- 
gators on  the  subject  of  seismic  sources  in  general  and  on  the  source  of 
of  the  San  Fernando  earthquake  in  particular. 

In  this  chapter  we  shall  briefly  review  that  small  portion  of  the 
literature  that  is  of  immediate  significance  to  the  evolution  of  the 
ultimate  conclusions  of  this  thesis.  This  review  is  by  no  means  exhaus- 
tive but  is,  it  is  hoped,  adequate  to  indicate  the  original  impetus  for 
this  research. 

1.2  Review  of  Previous  Studies 

Previous  solutions  to  the  inverse  problem  posed  by  the  large  and 
varied  data  set  of  the  San  Fernando  Valley  earthquake  have  incorporated 
mathematical  idealizations  of  many  forms.  Indeed,  several  important 
papers  discuss  the  problem  in  almost  purely  qualitative  terms,  and  many 
early  papers  take  a comprehensive  approach  using  data  of  several  differ- 
ent forms  to  compliment  one  another,  using  them  collectively  to  generate 
a complete  conception  of  the  faulting  process  of  the  event.  For  the 
purposes  of  this  introduction  it  is  useful  to  classify  previous  investi- 
gations as  being  essentially  quantitative  or  qualitative,  by  the  domin- 
ant frequency  range  of  ground  motion  considered,  and  by  the  distance  from 
the  source  region  of  the  data  analyzed.  Such  a classif ication  scheme  is 
largely  superficial  but  by  considering  first  those  papers  employing  long 
period  or  teleseismic  data  and  then  those  studies  using  increasingly 
higher  frequency,  closer  observations  it  is  hoped  that  the  method  of 
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this  thesis  may  be  seen  as  a logical  continuation  of  the  evolution  of 
the  understanding  of  the  seismic  source. 

An  important  attempt  aimed  at  satisfying  several  teleseismic  and 
static  manifestations  of  the  San  Fernando  earthquake  of  February  9,  1971 
is  undertaken  in  the  paper  of  Canitez  and  Toksoz  (1972).  These  authors 
wish  to  determine  this  event's  source  mechanism  characteristics  by 
analysis  of  the  polarity  of  teleseismic  P-  and  S-waves,  of  the  spectra 
of  surface  waves  (also  at  teleseismic  distances)  and  of  the  static  dis- 
placements, strains,  and  tilts  in  the  source  area. 

The  groundwork  necessary  to  the  calculation  of  near-field  static 
fields  has  been  profitably  researched  in  a number  of  papers  (Steketee, 

1958;  Chinnery,  1963;  Maruyama,  1964;  Mansinha  and  Smylie,  1971).  Using 
the  formulation  of  Maruyama  (1964)  which  may  be  integrated  numerically 
over  a fault  surface  composed  of  discrete  dislocations  of  arbitrary 
orientation  in  a half-space,  Canitez  and  Toksoz  calculate  the  pertinent 
static  fields.  The  solution  is  constrained  so  that  the  fault  surface 
is  planar  and  rectangular  with  the  sides  of  the  rectangle  being  parallel 
or  perpendicular  to  the  surface  of  the  earth.  A constant  amplitude  of 
dislocation  over  the  fault  surface  was  also  assumed.  Within  these  con- 
straints the  authors  find  that  the  static  strain  steps  recorded  at 
Isabella,  the  Nevada  Test  Site,  and  San  Diego  (Jungels  and  Anderson, 

1971;  Berger,  1971)  as  well  as  the  vertical  displacements  near  the  sur- 
face fault  trace  (Burford  et  al , 1971)  may  generally  be  fit  by  a single 
theoretical  model  whose  parameters  appear  in  Table  1-1. 

A dynamic  source  model  is  also  calculated  by  the  authors  on  the 
basis  of  polarity  of  first  P-wave  arrivals  at  one  hundred  twenty-eight 
stations  and  from  the  spectra  of  the  vertical  components  of  Rayleigh 
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wave  motion  at  six  stations.  In  the  first  motion  polarity  studies, 
angle  of  incidence  tables  of  Nuttli  (1969)  are  used  for  computing  the 
reduced  distances  using  a focal  depth  of  15  kilometers.  In  the  Rayleigh 
wave  study  the  average  continental  crust  and  upper  mantle  model  of 
Canitez  and  Toksoz  (1971)  is  used.  For  attenuation  corrections  constant 
Q values  are  used  in  two  intervals:  Q = 175  for  periods  of  from  10  to 

28  seconds,  and  Q = 200  for  longer  periods.  These  are  chosen  as  being 
in  the  mid-range  of  measured  Rayleigh  wave  Q values  of  Tsai  and  Aki 
(1969)  and  Solomon  (1971).  The  source  parameters  obtained  from  this 
part  of  the  study  are  also  summarized  in  Table  1-1  and  are  in  substan- 
tial agreement  with  the  statically  constrained  parameters. 

The  authors  feel  that,  using  extreme  values  allowable  for  Q,  the 
depth  of  the  initial  rupture  may  vary  by  ± 4 kilometers  without  des- 
troying the  distinctive  character  of  the  surface  wave  spectra  at 
several  stations. 

For  the  purposes  of  the  present  study  the  most  important  parame- 
ters derived  in  this  paper  are  the  dip  of  the  fault  plane  at  52  degrees 
and  the  focal  depth  of  approximately  14  kilometers.  The  primary  datum 
constraining  the  dip  angle  is  the  fault  plane  solution  derived  from 
first  motion  polarity  values.  As  Canitez  and  Toksoz  note,  this  datum 
is  most  pertinent  to  that  part  of  the  fault  where  rupture  is  initiated, 
inferred  to  be  at  the  lower  end  of  the  fault  plane.  The  focal  depth  of 
14  kilometers  is  supported  by  all  dynamic  evidence  within  the  uncer- 
tainty mentioned  above. 

One  difficulty  with  the  Canitez  and  Toksoz  study  is  the  teleseis- 
mic  distances  at  which  all  dynamic  data  are  evaluated.  Some  spherically 
symmetric  model  of  the  earth  is  implicit  in  all  teleseismic  data 
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inversion  schemes.  Any  departure  of  the  real  earth  from  this  ideal 
model  will  result  in  a biased  determination  of  the  source  parameters 
and  the  longer  the  wave  transmission  path  the  greater  will  be  the  like- 
lihood of  such  a departure.  Furthermore,  far-field  data  can  give  no 
direct  measure  of  the  ability  of  a source  model  to  explain  near-field 
ground  motions  which  are  of  great  theoretical  and  practical  importance. 
For  these  reasons  some  benefit  is  to  be  expected  from  the  consideration 
of  near-field  data. 

This  is  done  by  Mikumo  (1973)  as  part  of  another  comprehensive 
approach  to  source  parameter  determination  for  the  San  Fernando  earth- 
quake. Mikumo  carries  out  a detailed  investigation  of  the  source  mech- 
anism using  vertical  and  horizontal  displacement,  strain,  and  tilt  data, 
dynamic  qround  motions  in  the  near-field,  and  the  spatial  distribution 
of  aftershocks  and  features  of  the  surface  fault  break. 

The  study  of  the  static  displacements  undertaken  by  Mikumo  is 
essentially  identical  to  that  of  Canitez  and  Toksoz  (1972).  A slightly 
more  general  fault  geometry  is  considered  based  on  evidence  of  the 
spatial  distribution  of  aftershocks  (Allen  et  al,  1971;  Wesson  et  al, 
1971;  Hanks  et  al , 1971;  Allen  et  al , 1972;  Bolt  and  Gopalakrishnan, 
1975),  of  focal  mechanism  solutions  (Whitcomb,  1971;  Canitez  and  Toksoz, 
1972),  and  surface  faulting  (Kamb  et  al , 1971;  U.  S.  Geological  Survey 
Staff,  1971).  The  fault  surface  is  still  planar  and  approximately  rec- 
tangular, but  it  is  tilted  so  that  it  intersects  the  ground  surface 
obliquely.  This  has  the  effect  of  truncating  one  corner  of  the  fault 
plane  at  the  ground  surface,  the  strike  of  this  truncation  being  ap- 
proximately that  of  the  surface  expression  of  the  observed  fault. 

Mikumo  also  wishes  to  explore  the  effects  of  allowing  the 
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amplitude  of  dislocation  to  vary  over  the  fault  surface.  Since  pre- 
earthquake stresses  and  frictional  stress  during  faulting  may  not  be 
uniformly  distributed  with  depth  (Chinnery  and  Petrak,  1968),  several 
schemes  of  dislocation  amplitude  variation  are  considered  and  their 
resulting  static  surface  displacements  compared  with  those  observed 
near  the  fault  break.  In  all  cases  dislocation  amplitude  is  assumed  to 
decrease  with  depth,  either  linearly  or  exponentially  attaining  zero 
displacement  near  the  hypocenter.  Several  of  the  exponential  decay 
schemes  are  found  to  best  fit  the  observations.  Thus,  it  is  found  that 
static  displacements  may  be  fit  by  a roughly  rectangular  fault  of 
approximately  14  kilometer  length  and  19  kilometer  width,  intersecting 
the  ground  somewhat  obliquely,  with  the  dislocation  amplitude  having  a 
maximum  of  3.5  to  4.5  meters  at  the  surface  and  decreasing  exponentially 
with  depth. 

Mikumo's  dynamic  displacement  model  analysis  is  of  a much  differ- 
ent type  than  that  of  Canitez  and  Toksoz  (1972).  Here  an  attempt  is 
made  to  calculate  the  theoretical  dynamic  displacements  appropriate  to 
the  source-receiver  geometries  at  the  Pacoima  Dam  and  Pasadena  Seismo- 
logical  Laboratory  strong-motion  recording  sites  using  the  pertinent 
fault  parameters  developed  in  the  static  part  of  the  study  and  the  near- 
field displacement  formulation  of  Maruyama  (1963)  and  Haskell  (1969) 
which  apply  to  moving  dislocations  in  an  infinite  homogeneous  medium. 

It  is  reasonable  to  suppose  that  such  a transfer  function  may  success- 
fully be  used  to  model  direct  P and  S body  wave  motion  (with  the  excep- 
tion of  amplitude)  although  waves  that  depend  explicitly  on  the  exist- 
ence of  a surface  and  on  more  complex  medium  geometries  will  not  be 
adequately  modeled. 
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Some  parameters  necessary  to  the  consideration  of  a dynamic  model 
are,  of  course,  unrecoverable  from  static  data.  Such  parameters  are  the 
time  history  of  the  dislocation  at  any  given  point  on  the  fault  and  the 
fault  propagation  velocity.  Mikumo  assumes  a Haskell  (1969)  finite  ramp 
dislocation  time  history  with  rise  time  between  0.2  and  4 seconds.  Sev- 
eral different  modes  of  fault  rupture  are  considered  corresponding  to 
several  orientations  of  the  fault  slip  vector,  and  rupture  propagation 
velocity  is  allowed  the  discrete  values  of  2.0,  2.5,  or  3.0  kilometer/ 
second.  The  P-wave  velocity  in  the  epicentral  region  is  taken  to  be 
6.10  kilometer/second  from  seismic  refraction  studies  by  Wesson  and 
Gibbs  (1971).  The  S-wave  velocity  is  then  approximated  by  3.5  kilometer/ 
second  assuming  the  customary  Poisson's  ratio  of  0.25. 

The  somewhat  complex  modes  of  fault  propagation  chosen  by  Mikumo 
require  numerical  integration  of  the  appropriate  equations.  There  is, 
of  course,  nothing  in  the  Haskell  or  Maruyama  formulation  to  preclude 
contributions  to  the  ground  displacements  from  all  frequencies  includ- 
ing zero  frequency.  However,  since  he  wishes  to  compare  his  results 
with  the  data  of  Trifunac  and  Hudson  (1971)  which  represents  the  band 
limited  approximation  to  the  real  ground  displacement  between  essen- 
tially 0.04  and  15  seconds,  Mikumo  approximates  this  band  limiting 
process  by  applying  a boxcar  bandpass  filter  with  cutoff  frequencies  of 
25  and  0.05  Hertz  to  his  raw  ground  displacements.  In  this  way  direct 
comparison  of  data  and  theoretical  displacements  are  made. 

It  is  concluded  on  the  basis  of  such  comparisons  that  fair  agree- 
ment in  waveform,  arrival  time,  and  absolute  amplitude  is  obtained  for 
propagation  velocities  of  2.5  kilometer/second  or  slightly  less  and  for 
rise  times  between  0.6  and  1.2  seconds.  Generally  satisfactory 
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agreement  is  obtained  between  the  theoretical  and  observed  ground  dis- 
placements at  the  Pasadena  Seismological  Laboratory  using  similar  fault 
motion  parameters.  Appropriate  fault  parameters  derived  in  this  study 
are  summarized  in  Table  1-1. 

Trifunac  (1974)  undertakes  a very  similar  dynamic  motion  study  of 
the  implications  to  the  source  mechanism  of  various  San  Fernando  earth- 
quake strong-motion  records.  He,  however,  considers  the  strong-motion 
records  of  five  stations  (Pacoima  Dam;  8244  Orion  Blvd.,  1st  floor;  Jet 
Propulsion  Laboratory  basement;  Palmdale  Fire  Station;  Castaic  Old 
Ridge  Route)  and  allows  considerably  more  complexity  in  his  fault  geo- 
metry. The  fault  remains  planar  but  now  is  broken  into  irregular  rec- 
tangular subareas  whose  dislocations  are  allowed  to  differ  from  each 
other  in  both  amplitude  and  direction  within  the  plane  of  the  fault. 

The  medium  transfer  function  is  again  that  of  a homogeneous,  isotropic 
whole-space  as  developed  by  Haskell  (1969).  The  inversion  procedure 
itself,  however,  is  more  formalized  in  that,  having  specified  the  final 
fault  shape,  the  best  values  for  fault  propagation  velocity,  slip 
direction  and  amplitude,  and  the  rise  time  are  determined  in  the  least 
squares  sense. 

In  this  study  Trifunac,  as  Mikumo  before  him,  uses  near-field 
displacements  as  the  data  to  be  inverted  rather  than  velocities  or 
accelerations.  This  is  done  because  fewer  inadequacies  in  the  whole- 
space  medium  transfer  function  are  likely  to  manifest  themselves  at  the 
lower  frequencies  apparent  in  displacements  as  opposed  to  velocity  and 
acceleration  records.  This  contention  is  supported  by  several  recent 
studies  (e.g.,  Reimer,  1973;  Trifunac,  1973;  Boore,  1972;  Smith,  1975). 

The  conclusions  reached  as  a basis  of  this  study  are  listed  in  Table  1-1. 
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All  studies  considered  so  far  treat  the  fault  as  planar  with  a 
single  dip  angle  for  all  segments  of  the  fault.  There  was  growing 
evidence,  however,  that  the  fault  might  be  more  accurately  represented 
by  a curved  surface  with  dip  decreasing  as  successively  shallower  seg- 
ments of  the  fault  were  considered.  For  example,  fault  plane  solu- 
tions derived  from  first  motion  studies  indicated  an  initial  fault  dip 
angle  of  approximately  52  degrees  (Canitez  and  Toksoz,  1972;  Mikumo, 
1973)  while  surface  and  near-surface  fault  expression  indicated  a fault 
plane  of  much  shallower  dip  angle  (Kamb  et  al , 1971).  This  is  sugges- 
tive of  a curved  fault  surface  and  this  interpretation  is  consistent 
with  the  aftershock  data.  Static  dislocation  models  that  incorporated 
the  curved  or  hinged  fault  plane  concept  (Jungels,  1973;  Alewine  and 
Jordan,  1973)  were  also  shown  to  be  compatible  with  the  observed  static 
surface  displacements. 

Turnbull  and  Battis  (1973)  in  a continuation  of  earlier  work  be- 
gan by  Tsai  and  Patton  (1972;  1973)  brought  the  hinged  fault  model 
under  consideration  for  the  near-field  dynamic  analysis  of  the  Pacoima 
Dam  strong-motion  record.  In  this  study  comparison  is  made  between  the 
computed  Pacoima  Dam  velocity  waveforms  (Trifunac  and  Hudson,  1971)  and 
the  whole-space  theoretical  velocity  waveforms  computed  from  the  Haskell 
(1969)  formulation  as  in  previous  papers  considered  above.  The  hinqed 
fault  effect  is  obtained  by  linearly  summing  the  contributions  from 
several  planar  Haskell -type  faults  with  delay  times  properly  arranged 
to  simulate  continuous  fault  propagation  from  one  segment  to  the  next. 

Besides  orientation,  such  parameters  as  rupture  velocity,  dis- 
location amplitude,  and  fault  width  are  also  conceptually  allowed  to 
vary  between  contiguous  fault  segments  to  generate  a source  model  whose 


complexity  can  match  that  of  the  observed  velocity  record.  As  in  the 
study  of  Trifunac  (1974)  the  inclusion  of  such  complexity  is  an  advan- 
tage that  exacts  its  price.  The  final  solution  obtained  takes  advantage 
of  the  complete  generality  of  the  above  formulation  to  only  a limited 
extent.  This  solution  consists  of  two  planar  segments,  the  lower  dip- 
ping at  60  degress  and  the  upper  at  33  degrees.  The  lower  segment  is 
further  divided  into  two  co-planar  segments  while  the  upper  segment  is 
subdivided  into  three  co-planar  segments.  The  ranges  of  fault  width, 
longitudinal  (up  dip)  dislocation,  and  rupture  propagation  velocity 
among  these  segments  are  12.0  to  16.0  kilometers,  100  to  250  centi- 
meters, and  3.2  to  2.85  kilometer/second  respectively.  The  functional 
form  of  the  source,  as  in  all  previous  dynamic  studies,  is  the  Haskell 
(1969)  finite  ramp  whose  rise  time  is  found  to  be  0.6  seconds  for  all 
segments . 

Turnbull  and  Battis  feel  that  the  fit  of  their  generated  wave- 
forms to  those  observed  is  superior  to  previous  efforts  in  at  least  one 
respect:  the  initial  motion  of  the  large  pulse  on  the  vertical  record 

now  has  the  currect  polarity  with  their  fault  model.  This  is  purely  a 
result  of  the  geometry  of  the  above  model.  They  also  conclude  that  the 
rupture  velocities  of  different  segments  of  the  fault  are  not  well  con- 
strained by  this  approach  within  a moderate  range. 

It  must  be  noted  that  even  with  the  added  complexity  of  the  seg- 
mented fault  the  complexity  of  the  motion  after  the  large  S-wave  pulse 
is  not  well  represented.  This  could  be  due  to  the  failure  of  the 
Haskell  transfer  function  used  to  explicitly  consider  the  presence  of 
the  ground  surface,  any  layering,  or  surface  topography.  Alternatively, 
it  could  be  due  to  an  insufficient  number  of  fault  segments.  This 


possibility  is  not  explored  since  if  rupture  velocity,  rise  time,  dip 
angle,  dislocation  amplitude,  and  fault  width  are  allowed  to  vary  for 
each  of  a large  number  of  segments  the  inverse  problem  general  1 treated 
(even  within  narrow  limits  for  all  parameters)  soon  becomes  unmanageable. 

An  alternative  approach  to  extracting  the  sourch  mechanism  infor- 
mation contained  in  the  Pacoima  Dam  strong-motion  record  may  be  charac- 
terized as  semi-quanti tati ve . In  this  approach  an  attempt  is  made  to 
extract  implications  from  the  most  fundamental  data  of  seismology-- 
phase  arrival  times.  Two  papers  of  this  type  are  those  of  Bolt  and 
Gopalakrishnan  (1975)  and  Bolt  (1972).  In  these  papers  the  recorded 
components  are  resolved  along  directions  parallel  to  the  strike  of  the 
fault  plane,  down  dip  in  the  fault  plane,  and  normal  to  the  fault  plane. 
This  is  done  to  facilitate  the  analysis  in  terms  of  elastic  wave  theory. 
The  fault  is  modeled  by  a planar  surface  dipping  at  45  degrees.  Assum- 
ing P-  and  S-wave  velocities  of  5.5  and  3.3  kilometer/second,  respec- 
tively, the  salient  features  of  the  Pacoima  records  are  explained  in 
terms  of  direct  P-  and  S-wave  arrivals  from  the  hypocenter  and  a later 
S arrival  from  a point  on  the  fault  closes  to  the  Pacoima  Dam  site. 

Still  later  arrivals  are  conjectured  to  be  due  to  S-wave  arrivals  from 
more  distant  points  on  the  fault  surface.  The  contribution  of  surface 
waves  is  not  considered  important  in  the  near-field  since  their  far-field 
dominance  is  due  to  the  slower  attenuation  characteristic  of  surface 
waves  relative  to  the  attentuation  of  body  waves. 

Such  an  approach  is  very  appealing  in  some  respects.  Primarily,  it 
is  tempting  to  apply  intuition  and  knowledge  of  elastic  wave  behavior, 
gained  almost  exclusively  from  far-field  observations,  to  a near-field 
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problem.  However,  conclusions  reached  from  such  analysis  must  be  re- 
garded as  more  circumstantial  than  those  derived  from  more  straight- 
forward inversion  techniques.  This  contention  is  supported  by  a similar 
study  'undertaken  by  Hanks  (1974). 

In  this  study  a hinged  fault  plane  similar  to  several  considered 
before  is  favored  and  P-  and  S-wave  velocities  of  5.6  and  3.2  kilometer/ 
second,  respectively,  are  used.  As  in  the  previous  study  by  Turnbull  and 
Battis  (1973)  the  fault  geometry  is  found  to  conform  to  initial  motion 
polarity  requirements.  It  is  also  noted  that  the  orientation  of  the 
deeper  part  of  the  fault  is  such  that  the  P-wave  excitation  is  near  a 
minimum  while  the  S-wave  excitation  is  near  a maximum,  a circumstance 
that  explains  qualitatively  the  relative  compressional  wave-shear  wave 
amplitudes  apparent  on  all  components  of  the  Pacoima  Dam  strong-motion 
record.  As  in  the  previous  studies  of  Bolt  (1972)  and  Bolt  and 
Gopalakrishnan  (1975)  explanations  for  several  important  phase  arrivals 
on  this  record  are  sought  which  are  understandable  in  terms  of  elastic 
wave  arrivals  from  some  part  of  the  fault.  Besides  P-  and  S-waves  from 
the  initial  rupture  three  other  arrivals  are  investigated.  The  last  of 
these,  arriving  at  a time  of  approximately  7.3  seconds  after  the  com- 
pressional wave,  is  interpreted  as  the  breakout  phase  or  that  surface 
wave  phase  associated  with  the  rupture  of  the  earth's  surface.  No  very 
convincing  explanation  is  found  for  the  two  phases  between  this  breakout 
phase  and  the  initial  S-wave  arrival  although  the  second  of  these  is 
more  impressive  looking  on  the  acceleration  and  velocity  records  than  is 
the  motion  attributed  to  the  surface  rupture.  Based  primarily  on  the 
interpretation  of  this  breakout  phase  arrival  and  the  fault  geometry  a 
minimum  fault  propagation  velocity  of  3.8  kilometer/second  is  calculated. 


As  demonstrated  by  a comparison  of  the  conclusions  reached  by  Bolt 
and  Bolt  and  Gopalakrishnan  with  those  derived  by  Hanks,  small  uncertain- 
ties in  model  parameters  such  as  seismic  wave  velocities  or  fault  geome- 
try may  have  profound  implications  to  the  faulting  process  as  evidenced 
by  the  strona-motion  data.  Both  analyses  depend  primarily  on  the  inter- 
pretation of  three  arrival  times--ini tial  compressional  and  shear  wave 
arrivals,  and  one  later  phase.  In  the  case  of  Bolt  this  later  phase 
arrives  at  approximately  3.9  seconds  after  the  onset  of  the  recorded 
motion  and  is  explained  in  terms  of  an  S-wave  arrival  from  the  point  of 
closest  approach  of  the  fault  surface  to  the  Pacoima  Dam  Station.  In 
Hanks'  study,  as  we  have  seen,  a later  arrival  is  identified  as  the  sur- 
face wave  generated  when  the  fault  rupture  reaches  the  surface.  The 
phase  arrives  at  over  7 seconds  after  the  P-wave.  In  both  cases  the 
proper  identification  and  interpretation  of  this  later  arrival  provides 
a very  important  constraint  on  the  fault  rupture  propagation  velocity. 

The  inherent  difficulty  in  this  type  of  study  is  made  apparent  by  noting 
that  the  S-wave  arrival  time  chosen  by  Bolt  and  Bolt  and  Gopalakrishnan 
is  approximately  0.7  seconds  earlier  than  that  favored  by  Hanks.  This 
is  a significant  difference  over  the  short  ray  paths  pertinent  to  a 
near-field  study. 

We  have  now  briefly  discussed  a number  of  important  contributions 
to  the  understanding  of  the  faulting  process  associated  with  the 
February  9,  1971  San  Fernando  earthquake.  It  is  obvious  that  a great 
deal  has  been  learned  and  a qreat  deal  remains  to  be  learned.  For  ex- 
ample, little  has  been  determined  regarding  the  details  of  the  faulting 
episode  at  any  one  point  on  the  fault.  Is  the  time  history  of  faulting 
adequately  represented  by  a finite  ramp  displacement  function?  Is  there 
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hope  of  recreating  an  acceleration  record  from  some  fairly  simple  medium 
transfer  function  or  must  substantial  complexity  be  called  upon?  It  is 
the  purpose  of  this  thesis  to  find  the  answers  to  these  questions. 

As  we  have  seen,  previous  solutions  to  the  inverse  problem  posed 
by  the  large  amount  of  data  generated  by  this  earthquake  have  incorpor- 
ated mathematical  idealizations  of  many  forms.  Ultimately,  however,  the 
transfer  function  which  models  the  effect  of  the  real  earth  is  not  well 
known  even  over  the  short  distances  characteristic  of  the  near-field. 

The  ideal  near-field  transfer  function  would  be  one  that  could  adequately 
take  into  consideration  the  effects  of  a layered  half-space  composed  of 
material  not  necessarily  perfectly  elastic  and  perhaps  not  even  iso- 
tropic. It  would  also  allow  explicit  consideration  of  geometric  inhomo- 
geneities of  a fairly  extreme  nature  such  as  surface  topography  and 
material  differences  across  a dipping  fault  surface. 

The  medium  transfer  function  used  in  all  the  above  dynamic  near- 
field studies,  that  of  Haskell  (1969)  and  Maruyama  (1963),  may  be  im- 
proved in  a more  modest  fashion,  however,  by  generalizing  it  to  include 
explicitly  the  effects  of  the  ground  surface.  This  is  done  in  Chapter  2 
where  the  displacements  at  the  surface  due  to  line  forces  and  disloca- 
tions are  formulated  exactly.  This  work,  which  was  carried  out  by 
Dr.  Lane  Johnson  preparatory  to  his  more  general  three-dimensional  formu- 
lation (Johnson,  1974)  of  the  same  problem,  will  allow  the  straight- 

9 

forward  inclusion  of  Rayleigh  waves  in  the  model  interpretation  of  the 
recorded  data  as  well  as  one  body  wave  phase  whose  existence  depends  on 
the  surface--a  wave  that  leaves  the  source  as  an  S-wave  and  is  reflected 

* along  the  surface  as  a P-wave.  The  two-dimensional  nature  of  the  formu- 

J lation  of  Chapter  2 requires  that  the  data  considered  be  close  enough  to 
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the  source  that  this  source  may  be  treated  as  two-dimensional.  That  is, 
the  station  must  be  close  enough  to  the  source  so  that  the  distance  be- 
tween source  and  receiver  is  small  relative  to  the  dimension  of  the 
source.  This  requirement  is  fulfilled  by  the  position  of  Pacoima  Dam 
relative  to  the  zone  of  energy  release  of  the  San  Fernando  earthquake 
and  so  it  is  these  records  that  are  used  as  the  observational  basis  for 
comparison  with  generated  model  accelerograms. 

In  Chapter  3 those  additional  parameters  necessary  to  the  genera- 
tion of  a model  accelerogram  are  considered.  They  include  the  instrument 
response  of  the  accelerometer  through  which  the  ground  motion  is  record- 
ed, the  time  history  of  the  source  at  any  point  on  the  model  fault  sur- 
face, the  geometry  of  some  two-dimensional  fault  relative  to  the  location 
on  the  surface  at  which  acceleration  is  to  be  evaluated,  the  fault  rup- 
ture propagation  velocity,  and  the  amplitude  of  the  final  displacements 
of  the  sides  of  the  fault  relative  to  each  other  as  a function  of  posi- 
tion along  the  fault.  As  we  shall  see  in  Chapter  3 it  is  necessary  to 
consider  ground  motions  at  frequencies  high  relative  to  the  rise  time  of 
the  source  time  history  in  order  to  have  any  opportunity  to  distinguish 
the  functional  smoothness  of  this  time  history.  From  previous  work  we 
have  seen  that  derived  rise  times  of  approximately  0.5  seconds  are  com- 
monly found  so  that,  by  the  above  contention,  ground  motions  of  2 Hertz 
or  higher  frequency  must  be  considered.  Such  high  frequencies  are  most 
characteristic  of  accelerometer  records  and  this  is  one  reason  why  these 
records  are  chosen  as  the  data  base  rather  than  their  velocity  or  dis- 
placement equivalents  which  may  be  derived  from  them.  This  study  is  the 
first  quantitative  attempt  to  explain  an  actual  observed  accelerometer 
record.  With  the  exceptions  of  basic  review  in  sections  3.2  and  3.5  all 
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work  in  Chapter  3 is  original  as  is  the  work  in  succeeding  chapters. 

In  Chapter  4 the  results  of  the  previous  two  chapters  are  combined 
to  allow  the  generation  of  a series  of  model  accelerograms  which  may  then 
be  compared  with  the  recorded  accelerations  of  the  Pacoima  Dam  strong- 
motion  site.  The  degree  and  manner  of  the  success  or  failure  of  such 
model  accelerograms  to  match  the  observations  are  used  to  draw  inferences 
about  the  faulting  process  of  the  San  Fernando  earthquake. 

Finally,  in  Chapter  5 conclusions  are  drawn  and  areas  for  further 
exploration  briefly  mentioned. 
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Chapter  2 

The  Green's  Functions  of  the  Problem 
2.1  Some  Basic  Considerations 

Consider  the  general  linear  differential  equation 
Ly(x)  = f(x) 

where  L is  a known  linear  differential  operator,  f(x)  is  some  known  func- 
tion and  y(x)  is  to  be  determined.  The  inverse  of  a differential  operator 
is  an  integral  operator  which  we  may  symbolize  by 

L''f(x)  = /G(x,fjf  (C)a£ 

provided  that  the  kernal  G(x,^)  can  be  chosen  so  that 
LL'1f(x)  = f(x)  . 

Thus  we  wish 

LL*1  f (x)  = L/G(x,e)f(C)dC  = /rLG(x,5)]f(Ud5  = f(x)  . 

This  condition  can  be  easily  fulfilled  if 

LG(x,£)  = <S(x-0  2.1-1 

where  6(x-£;)  is  the  Dirac  delta  function  defined  such  that 

«(x-C)  = f0’  X * K 
t°°»  x = E, 

and 

/6(x-C) f(C)dC  = f(x) 

over  any  interval  that  includes  x.  When  G(x,f;)  has  been  found  (where 
G(x,C)  is  the  general  solution  to  2.1-1)  the  solution  to  the  original 
equation  is 

y(x)  = L " 1 Ly ( x ) = L"1 f (x ) = /G(x,C)f(Od?  . 

The  particular  kernal  with  this  property  is  known  as  the  Green's  func- 
tion for  the  general  linear  differential  operator  L. 

For  the  problem  to  be  attempted  here  the  linear  differential 
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operator  of  interest  is  that  appropriate  to  the  equation  of  motion  for 
a homogeneous,  isotropic,  perfectly  elastic  solid.  Initially,  however, 
the  more  general  case  of  an  inhomogeneous,  anisotropic  elastic  solid  will 
be  considered.  The  equation  of  motion  then  is  (using  the  standard  nota- 
tion as  in  Burridge  and  Knopoff,  1964) 

[ c i j pq (x ) up<q(*’t)l,j  - p(x)  u i ( x , t ) =-fi(x,t)  2.1-2 

where 

x = x-j^i  + X2$2  + x^  is  the  Cartesian  position  vector, 
u.j(x,t)  is  the  ith  component  of  displacement  at  x and  at  time  t, 
f.(x,t)  is  the  i^  component  of  the  body  force  at  x and  t, 

C i i (x)  are  the  elastic  constants  of  the  medium,  and 

' J r H 


Up,q^x,t)  = 9T  Up(x’l) 


U.-(x,t)  = u.  (x,t)  . 

at 

The  summation  convention  is  in  effect  so  all  repeated  indices  are  to  be 

summed  from  1 to  3.  It  may  be  shown  from  conservation  principles  that 

even  in  the  most  general  elastic  solid 

C. • = C..  = C . . . 2.1-3 

ijpq  J l pq  pq-ij 

Now  let  Vj(x,t)  be  another  motion  due  to  a body  force  g^(x,t).  Then  in 
an  inhomogeneous,  anisotropic  elastic  solid  we  may  write 

lCiinn(X)Vn  " p(x)v.(x,t)  = - 9 - (x  , t)  2.1-4 

These  equations  (2.1-2  and  2.1-4)  are  satisfied  within  a volume  V bounded 
by  the  surface  S for  all  time. 

Next  we  wish  to  require  that  f^(x,t)  and  g^(x,t)  vanish  for  some 
time  t < - T where  T is  a constant.  Then  from  the  principle  of  causality 
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v.j(x,t)  and  iij(x,t)  must  also  vanish  for  t < - T.  In  2.1-4  replace  t by 
- t so  that 

[Cijpq(x)vp,q(**t)]»j  ' p(x)v(x,t)  = - g.(x,t)  2.1-5 

where 

vi  (x,t)  = v • ( x » - 1) 
g^.t)  = g^x.-t) 

and  g.(x,t)  = (x,t)  = 0 for  t > T . 

We  now  multiply  2.1-2  by  (x,t)  and  2.1-5  by  u^(x,t),  subtract  the  pro- 
ducts and  integrate  over  the  volume  V for  all  time  t. 

-^dti{[Cijpq  Up5q^’j  Vi  " ^Cijpq  VP»q^’j  Ui 


[pli^v.  - p v^.]  }dV  = 
-T^lUiSi  - V.f.JdV 


2.1-6 


We  may  immediately  dispose  of  the  second  term  on  the  left  of  equation 
2.1-6  since 

_/°°dt/[pUi ^ - pv.jU.jJdV  = 

-<£>/P  rui v.  - v.u-ldV  = 

• t = °° 

i I 

V t = - 00 

But  u^  vanishes  for  t < - T,  and  v^  vanishes  for  t > T so  the  integrand 
is  identically  zero. 

Next  using  the  expressions 

(C.  . u I , . v.  = I C . 

1 ijpq  p>q”j  i 1 i 


U V.  , . 

ijpq  P.q  1 J 


- [C . . u | v.  . 
ijpq  p ,q  i,j 


and  similarly 


Li 


I" 


» V*4  * . , *Vr*V'r  • 


> ^ * * 4 


J 
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1Cijpq  Vp,qJ’j  Ui  “ lCijpq  Vp,q  Ui 1 ’ j ' 1 C i j pq  Vp,q]ui,j 


We  may  write  the  first  term  in  2.1-6  as 


f°dtf{  [C . . u ],.  v.  - [C..  v u.}dV  = 

y 1 ijpq  p.qJ  j i t jpq  P,qJ’j  r 

/°°dt/{  [v.C.  . u - u.C.  . v 1 , . - 

-«  v i upq  p»q  i ijpq  p.qJ  j 

c. . fv.  .u  - u.  .v  ] } d v 

i jpq 1 i,j  p»q  i,j  p.qJ 


2.1-7 


The  last  term  of  this  expression  vanishes  because  C..  = C ..as  indi- 

ijpq  pqij 

cated  in  equation  2.1-3.  That  is,  since  ij  and  pq  are  repeated  and 

therefore  dummy  indices  under  the  summation  convention,  they  may  be 

interchanged  at  will  so  that 

C..  v.  .u  -C..  v u.  .= 

ijpq  i,j  p,q  ijpq  P,q  i ,J 

C.  • v.  .u  - c . ."v.  .u 

-ijpq  i,j  p,q  pqij  i ,J  p,q 


C. • v.  .u  - C . • v.  .u  =0 
-ijpq  i,j  p,q  ijpq  l ,J  p,q 


2.1-8 


From  these  considerations  and  from  equations  2.1-6  and  2.1-7  we  have 
-fdWVi  Cijpqup>q  - UiCijpqVp^j.jJdV 


=_/  dt/[uigi  - v.f.]dV 


2.1-9 


We  may  now  use  the  divergence  theorem  on  the  left  side  of  this  expression 
to  derive 

oo  — — 

/ dt/n  • I v.C • • u - u-C • • v ldS 

s J 1 1Jpq  p’q  1 1Jpq  p’q 

= jfdt/ru^  - v^f^dV  2.1*10 

In  this  expression  n.  is  the  outward  pointing  normal  to  the  surface  S. 

J 

Equation  2.1-9  is  a form  of  the  reciprocity  theorem  of  Betti  and  Rayleigh. 

We  may  now  show  how  to  solve  for  u^  in  terms  of  its  Green's  func- 
tion. As  suggested  by  equation  2.1-1  let 


x.  ^ h ...  . * \ 
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gi (x,t)  = 6ni  6(x-x ' ) 6 ( t- t 1 ) 


2.1-11 


where  6 . is  the  Kronecker  delta 
m 


ni 


o,i7n 
1 ,i=n 


Substituting  into  equation  2.1-10  and  integrating  we  have 

oo  — >■  -y  > 

un(x',t')  = dt/Gin(x,-t;x,,-t')fi(x,t)dV 


+ -rdt{nj{Gin(^'t^'’-t,)CiiDo(?)un.0(^t) 


upq  ' p.q 


' Ui(^’t)Cijpq{^)Gpn,q(^’_t^,*"t,)}dS 


2.1-12 


In  this  expression  we  have  written  G..  instead  of  v^  to  symbolize  the 
solution  to  equation  2.1-4  with  g^  defined  as  in  equation  2.1-11.  This 
is  in  conformity  to  the  notation  of  equation  2.1-1.  G.jn(x,t;x 1 , t * ) is 
the  displacement  in  the  i direction  at  (x,t)  due  to  an  impulsive  point 
force  in  the  n direction  at  (x’,t').  If  in  equation  2.1-12  ui  and  G^ 
satisfy  the  same  homogeneous  boundary  conditions  on  S and  if  we  require 
f.  to  be  an  impulsive  point  force  we  may  derive  the  general  form  for  the 
Green's  function  reciprocity  theorem  for  the  differential  equation  2.1-2 
Gj j(x,t;x' » t * ) = Gjj (x1 ,-t1 ;x,-t)  2.1-13 

Interchanging  the  variables  x with  x'  and  t with  t'  and  using  equation 
2.1-13  we  may  write  equation  2.1-12  as 

un(x,t)  = _/"dt'/Gni(x\t;x,,t,)fi(x,,t')dV' 


+ -^Td^ ' ^ Gn i ( * ’ * ' » ^ ' )Ci  iDn^'  ^ un .n^ * ' ’ ^ ^ 


i JPq  P,q 


Ujlx'.t'JCjj  (x')G  2.1-U 


ijpq  ' np ,q ' 

where  integration  is  now  with  respect  to  the  primed  variables  and  where 
Gnp  ^q(x  ’ t;x  ’ ’ ' ) 3 9-Fq-Gnp(X’t;x,’t') 


• ^ 
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Of  course  the  summation  convention  still  applies.  Equation  2.1-14  is 
the  generalized  representation  theorem  of  de  Hoop  and  Knopoff.  This  ex- 
pression will  give  us  the  displacements  at  any  point  x at  time  t due  to 
forces  within  or  at  the  boundary  of  volume  V'. 

Let  us  now  find  u^  for  prescribed  discontinuities  in  the  displace- 
ment and  its  derivatives  across  a surface  E embedded  in  V.  Let  Vj  be 
the  unit  normal  to  E and  let  ( u - ( x * , t ' ) ] and  [ un  (x',t')]  be  the  dis- 

continuities  in  u-  and  u across  E.  Assume  u-  and  G . to  satisfy  the 
i p }C|  i n l 

same  homogeneous  boundary  conditions  on  S and  apply  equation  2.1-14  to 
the  region  bounded  internally  by  I and  externally  by  S.  G . does  not 
have  prescribed  discontinuities  on  E.  Then  the  surface  integrals  over  S 
are  replaced  by  integrals  over  E only.  Let  us  further  require  that  there 
be  no  body  forces  in  V and  no  stress  discontinuities  across  E so  that 
f_j(x',t')  = 0 for  x'  in  V,  and 

lun  n(x' ,t')]  = 0 for  x'  on  E 

r »M 

Equation  2.1-14  now  becomes 

un(x,t)  = jfdt'/Vjll  u^x*  , t:' ) lci  jpq(x' ) • 

GnnJx.bxM'lldr  2.1-15 

nP 

where  x’  is  made  to  conform  to  the  dislocation  surface  E and  the  spatial 
integration  is  over  this  surface.  It  is  obvious  that  to  proceed  further 
we  will  need  the  solution  for  the  Green's  functions  themselves,  or  more 
specifically  for  the  purposes  of  evaluating  equation  2.1-15,  wo  will 
need  the  spatial  derivatives  of  the  Green's  functions.  These  will  be 
found  below. 

2.2  The  Green's  Functions 

Let  us  consider  a line  source  buried  in  a half  space  as  shown  in 


Figure  2.2(1).  Let  us  further  restrict  the  forces  acting  along  this 
line  to  be  in  a direction  perpendicular  to  this  buried  line  source.  If, 
as  in  Figure  2.2(1)  the  line  is  oriented  in  the  y direction  at  x'  =0 
and  at  a depth  h below  the  surface  of  the  half  space  the  forces  are  con- 
strained to  act  in  the  xz  plane. 

Then  from  the  conservation  of  linear  momentum  and  assuming  that  the 
half  space  is  homogeneous,  isotropic,  and  perfectly  elastic  we  may  write 

p -^5  u(x,z,t)  = f(x,z,t)  + pV2u(x,z,t) 

ar 

+ ( A+p)\)(^*u(x,z,t) ) 2.2-1 

where 

u(x,z,t)  = v(x,a,t)  x + W(x,z,t)z, 
f(x,z,t)  = [e(t)x  + g(t)z]6(x)6(z-h), 
p is  the  density  of  the  medium, 
p and  A are  Lame's  elastic  constants, 
v is  the  Laplacian  operator, 

'?  is  the  gradient  operator,  and 
V-  is  the  divergence  operator. 

Carrying  out  the  operations  in  equation  2.2-1  and  resolving  into  compon- 
ents we  may  write  a matrix  equation  of  the  form 
Ax  = b 

where 


f(u2u)  72  + "72  - 7721  {(u")3'x'a¥ 
ax  az  at 


,2  ,2  ,2  ,2 

f(A+n)  Ua+2m)  7 + i1-?  ‘ p— p) 

' az  ax  at 


x = 


v( x,z,t) 
w(x,z,t) 


b = 


-e(t)  6(x)  iS(z-h) 
-g(  t.)  <S ( x ) (S(z-h) 


. 


. ♦*>*  . 


IT  • r '• 


26 


Now  let  us  take  the  Laplace  transform  of  this  matrix  equation  in  the 
order  t-*-s,  x-»-C*  and  The  Laplace  transform  with  respect  to  t is  the 

ordinary  one-sided  Laplace  transform  while  those  with  respect  to  x and  z 
are  the  generalized  two-sided  Laplace  transforms.  In  later  development 
we  will  require  s to  be  real  and  positive.  With  this  requirement  the  re- 
maining two  transform  variables  must  be  restricted  as  outlined  below  in 
order  to  assure  convergence  of  the  transform  integrals. 

If  s is  real  and  positive  with 


cx 


1/2 


B = 


[H] 


1/2 


K = - 

a a 


K - — 

Ke  3 


then  in  taking  the  x->-£  transform  we  find  that  we  must  have 
I Refill  < minimum  (K  ,K„) 

CX  p 

where  |Rp{^}|  is  the  absolute  value  of  the  real  part  of  £.  If  we 
further  define 

, ,1/2  , 9 1/2 

‘ [*«  -e  1 • ^ -s2] 

we  see  that  the  above  restriction  on  Z,  requires  that 


Ref v 1 > 0 

a - 


RefVg}  > 0 


Then  in  taking  the  z-><;  transform  we  must  require  that 
| Ref U I < minimum  [Refv^l ,Refvp} | 

When  the  transforms  have  been  taken  under  the  above  conditions  and  in 
the  given  order  we  know  that  the  Laplace  transform  of  A exists  and  we 
may  write 

Ax  = b 2.2-2 


where  A = AfAi  and  similarly  for  the  other  terms  and  where 


* 
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A = 


x = 


{(A+2u)£2  + uc2-  ps2)  (A+u)£A 
(A+u)a£  {(A+2p)i;2  + y^2-ps2}  I 


V(C,4,s) 

WU.A.s) 


b = 


-E(s)e_ch 


-G(s)e 


-Ah 


where  E(s)  = L{e(t)}  = /°°e( t)e  Stdt  and  similarly  for  V,  W,  and  G. 

0 

From  elementary  matrix  theory  it  may  be  shown  that  the  system  of 
algebraic  equations  symbolized  by  the  matrix  equation  2.2-2  may  be  solved 
if  the  determinant  of  the  2x2  coefficient  matrix  A is  not  zero.  That  is 

2.2-3 


x = A_1b 

where,  from  equation  2.2-2 

A'1  = 1 
- D 


(uC2  + (A+2p)c2-pS2}  - (A+u)AC 
-(A+p)CA  {(A+2q)C2  + MA2-PS2} 


2.2-4 


Here  D is  the  determinant  of  the  matrix  A and  may  be  written 

D = |(A+2uk2  + UC2-pS2]  [(A+2yk2  + M^2-PS2]- [( A+u)cej2 
We  wish  to  write  this  in  such  a way  that  it  will  be  meaningful  when  the 
inverse  Laplace  transform  of  the  variable  A is  taken.  With  some  alge- 
braic manipulation  D may  be  rewritten  as 

D = p(A+2p)[(A2-va2)(A2-v62)]  2.2-5 

We  have  now  solved  the  problem  in  the  transform  domain.  It  is  now  neces- 
sary to  return  to  the  real  domain. 

2.21  The  c->z  Inverse  Transform  and  Boundary  Conditions 

Let  us  consider  the  V(f,r,s)  term  of  the  matrix  equation  2.2-3.  This  may 

be  wri tten  out  as 


VU.a.s)  = ’ {(tjf2+(A+2u)A2-ps2)(-E(s)e 


-Ah, 


+ (A+p)r,-.G(s  )e 


- Ah , 


2.2-6 


..  ■ «■»* 


. . 


-*  tr  • 
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Then  since 


L [ F ( C ) ] = f(z)  = f 

-i°° 


we  have  from  equation  2.2-6 


VU. z,s)  = 2^-  { [J-{(u?2+(A+2m)?2-ps2)(-E(s)) 

-ioo 

+ (A+p)UG(s)}e^(h‘z)]ds 


2.2-7 


From  this  expression  it  is  seen  how  intimately  the  determinant  con- 
trols our  problem  since  it  specifies  the  poles  whose  residues  determine 
the  value  of  this  integral.  From  equation  2.2-5  it  is  clear  that  the 
integrand  of  equation  2.2-7  has  simple  poles  at 


C = ± \> 


C ~ ± Vr. 


There  are  two  cases  to  consider.  For  z<h  we  have  h-z>0  and  for  z>h  we 
have  h-z<0.  Thus  for  the  integral  to  exist  in  equation  2.2-7  we  must 
close  the  contour  in  the  right  half  plane  of  5 space  for  z<h.  This  im- 
plies residues  from  the  poles  at  v and  v„.  Note  that  this  contour  is 
closed  clockwise  so  that  an  additional  minus  sign  is  necessary.  For  z>h 
we  close  the  contour  in  the  left  half  plane  picking  up  the  residues  from 
the  poles  at  -v^  and  -v^.  These  two  cases  are  shown  in  Figure  2.2(2). 


For  z<h  we  have  the  two  residues 


-vjh-z) 


ResK=vJ  = C1|r2E(s)+CvaG(s)l-e- 


where  Cj  is  a constant  such  that 


C = t (K  2-K  21 

L1  2pf A+2u)  K6  ' 


Similarly 


-rv a 
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-vp(h-z) 


Res [r=v  | = -C-j  [v  2E(s)-Cv„G(s)]! 


From  equations  2.2-7,  2.2-8,  and  2.2-9  we  have 


VU.z.s)  = 2^f  (-2iri){Res[t*wiJ(l  + Res[r.=vf 


2.2-9 


-va(h-z) 

V(£,z,s)  = C1{(F2E(s)+FvrYG(s))? v 


+ (vB2E(s)-FvgG(s))^ 


-v6(h-Z) 


-} , z<h 


This  is  one  quarter  of  the  complete  solution  in  (C,z,s)  space.  Now  we 
wish  to  find  W(£,z,s)  for  z<h,  and  V(£,z,s),  W(£,z,s)  for  z>h.  The  es- 
sentials of  the  derivation  of  these  quantities  are  the  same  as  those  used 
above.  Therefore,  we  write  without  proof  the  whole  solution  below.  In 
this  matrix  equation 

'l  ,h>z 

*9"<h-z)  ' -,,h<z 


V(£,z,s) | ' If2  fA>asgn(h-z) 

-v  ti- 
nt' 

e 

W(r,z,s)j  1 Jfvasgn(h-z)  v^2 

Va 

h * 

uB2  -FvpSgn(h-z)  e’Vh_Z'' 

E(s) 

j -FVpSgn(h-z)  F2  “P 

G(s) 

Having  solved  for  this  particular  solution  we  wish  to  find  the  compli- 
mentary solution  to  the  homogeneous  equivalent  to  equation  2.2-2  in 
( ,z,s)  space.  That  is,  we  wish  to  find  the  solution  to  the  set  of 
equations 


» Vv*'.  — - * **>  , ■ »*»»»  , . •»  ‘ X 
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{ (x+2u)c2  + u— ~ - PS2}  (A+y)c|r 

SI  * 

V(5,z,s) 

0 

2 

(A+u)^  ((A+2y)-^  + uC2  - ps2} 

3Z 

W(?.z,s) 

0 

These  equations  may  be  shown  to  have  solutions  which  remain  finite  as  z 
approaches  +°°  of  the  form 


-V(1Z  -VpZ 

V(£;,z,s)  = Af,e  + Bv,,e 

-v  z -vRz 
W(S,z,s)  = - Avae  a + BCe  K 


2.2-13 


That  these  are  the  only  nontrivial  solutions  to  equation  2.2-12  may  be 
seen  from  direct  substitution.  A and  B are  constants  in  the  sense  that 
they  are  independent  of  z.  Otherwise  A and  B are  not  restricted.  Later 
algebra  will  be  greatly  simplified  if  we  write 

-v^h 

A = C, [A  E(s)  + A G(s) ]— 

a a 


-v„h 


B ? C]  lBeE(s)  + BgG(s)|-e~ 


(3 


Then  the  general  solution  to  the  problem  may  be  written  in  matrix  form  as 


|V(f.,z,s) : 
W(f.,z,s 


= C, 


I"2  sgn(h-z) 
L I- v sgn(h-z)  o 


-v  |h-z) 

„ a' 


A 


g ! e 


-\  A -v  A 
, e < g 


’ ' jv  2 - '.v  sgn(h-z) 

2 

< sgn(h-z) 


-vp I h-z  | 


\>n 


V:-Bg  j e VV  (h+z) 
B_  f.B 


e • g | ‘ J G(s)1 


| E ( s ) | 


2.2-14 
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where  Ag,  Be,  A^,  and  are  constants  in  the  same  sense  as  A and  B and 
must  be  determined  from  boundary  conditions.  Since  the  solution  is  now 
in  (C,z,s)  space  let  us  briefly  examine  boundary  conditions  in  this 
space . 

As  is  usually  done,  let  the  stress  on  a surface  whose  normal  is  in 
the  i direction  due  to  a force  in  the  j direction  be  Thus,  if  we 

wish  to  specify  the  stresses  along  a surface  of  constant  z (or  one  whose 
normal  is  in  the  z direction)  we  need  to  find  r , t , and  r . 

L X Lj  L L 

Let  us  further  define  the  strain,  e^,  as  that  strain  in  the  i 
direction  due  to  a displacement  in  the  j direction.  Then,  for  the  type 
of  medium  under  consideration  it  may  be  shown  that 


, !)u.  au, 

e.  . = I ( — 1 + — J-) 

1 j 2 v3Xj  3x- ; 


where,  as  before,  u.  is  the  displacement  in  the  i direction.  Then 


Tij  ‘ “to"  * 2,,e'j 


2.2- 


where  0 


ekk  = ell  + e22  + and  is  the  Kronecker  delta  as  previ 


'33  ij 

ously  defined.  In  previous  notation  u-j  = u(x,z,t)  and  u^  = w(s,z,t). 
Then 


exx  ■ H • eyy’0=ezy  • ezz  ’ Tz  "(«'z't) 

ezx  = 2 ^iz  u^x’z’t)  + w(s,z,t)) 

and  from  equation  2.2-15 

Tzx  = (£  u(x,z,t)  + |x  w(s,z,t)) 


zy 


= 0 


I zz  = u(x’z’t)  + 


+2i0  ' w(s ,z,t) 
<1  z 
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From  equation  2.2-18  it  is  easily  seen  that  for  t (£,o,s)  and  t (f,o,s) 

Z X z z 

to  be  zero  for  arbitrary  F(s)  and  G(s)  we  must  have  the  four  conditions 
2f;U-Ae)e  + ^ (vg  -C)(vg-Be)e  6 = 0 

(vB  -r)(C+Ae)e  a - 2C(Vg+Be)e  =0 

B 

2ava-A?)e"V“h  - (vg2-f,2)(r+Bg)e'VBh  = 0 
f (vgV)(va+A  )e'V«h  + 2C(f,-B  )e'Vf?h  = 0 

a a 3 

Solving  these  two  sets  of  two  simultaneous  equations  in  two  unknowns  we 
find 


Ae  = ‘ R * 4^W4VB(V-^)e  “ " 1 


2 ,-2x2 


2 ,2  <VvB,h 


2 ,2x2  ,,2.  ^r2,„  2_f2)e(VvB,h| 

f r,  2 ,2x2  .,2  , 2 r2xJVV«)l' 


a — oc  r/  2 r2\ 2 /| r2  j-/irw  2 r2 \ 

Ag  " - T t(v3  ^ } " ^ V3  4C  VB  )e 


R 1 v 8 ' ' ’ a 

2_r2l2  . /.A 


5g  = - R ^B  “«  > - 4f-  Vb’4Vb(v8  -r'  )e 

where  R = (vp2-f2)2  + 4C2vavg. 

We  may  now  return  to  equation  2.2-14,  set  z=0  and  substitute  for  A , A , 
Be,  and  Bg  to  solve  for  the  general  solution  at  the  surface  in  (C,z,s) 
space.  That  is,  we  have  now  solved  for  a line  source  in  a half-space 
with  the  correct  boundary  condition  in  the  transform  space.  From  equa- 
tion 2.2-14  and  the  above  expressions  we  have 


V(f.,o,s) 

r 

- 

zr~\.  2fv,< 

w(f ,o,s) 

R 

. , ? . 2,  , 2 .2. 
T(vp  ) vt(vg  -f  ) 

-v  h 

a 


r 
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vg(vg2-C2)  -S(v62-A 


-2^v  v„ 
^ ot  6 


9,2 

2£ 


-v 


E(s) 

G(s) 


2.2-19 


2.22  The  £+x  and  s-*-t  Inverse  Transforms 


In  taking  these  last  two  inverse  transforms  we  shall  be  using  the  method 
of  Cagniard  (1939)  as  modified  by  de  Hoop  (1960).  This  method  involves 
the  playing  of  a pair  of  inverse  transforms  against  one  another  so  that 
the  evaluation  of  either  of  the  inverse  transforms  becomes  unnecessary. 
That  is,  we  wish  to  manipulate  two  integrals  in  such  a manner  that  we  may 
find  an  expression  that  is  recognized  as  the  inverse  transform  of  a trans- 
form so  that 

v(x,o,t)  = {V(C,o,s)}^x)s^t 

= i'1{i{v(x,o,t)}t_+s)s^t 

Thus,  it  is  necessary  to  find  under  what  conditions 
r-l 


L {vu,  o,s)}  = l{v(  x,o , t) ) t^s 


2.2-20 


If  such  conditions  can  be  found  the  argument  of  the  Laplace  transform  on 
the  right  hand  side  of  eguation  2.2-20  will  be  the  desired  answer. 

To  begin  with,  let  us  write  the  £->x  inverse  transform  of  equation 
2.2-19  explicitly  as 

. i»,  9 -v  h oo  — v^h 

V(x,o,s)  = ^2ttT  ( ^ vBe  +VB^VB 

-v  h oo  ~^nh 

+ [2eVpe  “ -C(vp2-C2)e  L ]G(s)}e-xdC 


or 


V(x, o,s)  = / l{[2f,2vflE(s)+2f,v3vaG(s)le 


£x-v  h 


g2iri_^R 


o o p o f.x-v„h 

+.  |v3(v62-f2)E(s)-i;(v6  -C2)G(s)le 


}df 


Li 


\ V'W*  *.  ^ * % - 


In 


2 1 

where  we  have  used  the  fact  that  2Kp  Cj  = — . Similarly 

i o o oo  £x-v,  h 

W ( x , o , s ) = y27n  [ t £ ( vg  "5  ) E ( s ) +va  ( vp  - £.  ) G ( s ) J e 

o £x-vfih 

- [2CvavBE(s)-zrvaG(s)]e  * }d^ 

Let  us  now  consider  the  problem  in  polar  coordinates  (r,0)  where 

op  -lx 

x = rsin0,  h = rcose  and  conversely  r = /x  +h  , 6 = tan  (pj-).  Also  let 
E,  = sy,  va  = sna>  and  vp  = snp  where  from  previous  requirements  on  vq  and 
Vg  we  have 

na=/-^-Y2  , Refn}>0 

a 

Hg  = / -L-  - y2  , Refng}  > 0 

Making  these  substitutions  into  the  above  expressions  for  V(x,o,s)  and 
W(x,o,s)  we  have 

i o s(yx-n  h) 

V(x ,o,s)  = -2^-  / ^([2y  nBE(s)+2YnanpG(s)le 

00  oo  s(Yx~nRh) 

+ [VV"Y )&($)]*  P 

1 i°°i  o o op  s(YX-nh) 

W(x,o,s)  = — { pt fY(nB^-Y  ;E(s)+n  (nB^-Y2)G(s)le 

p s(YX-nf.h) 

- f2YnangE(s)-2Y^naG(s)]e  )dy  2.2-21 

where  now  R = (np2-Y2)2  + 4y2nanB. 

Let  us  now  consider  a contour  with  the  following  properties  in  the 
complex  y plane:  this  contour  is  symmetric  about  the  real  y axis  (Imlyl 

= 0)  and  one  part  of  it  is  the  line  Refyl  = C where  C is  a positive  con- 
stant. Such  a contour  is  shown  in  Figure  2.2(3).  We  know  from  Cauchy's 
Integral  Theorem  that  if  this  contour  encloses  no  s ingul ari ties  the  total 
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contribution  of  an  integral  along  this  contour  is  identically  zero. 
Furthermore,  if  all  singularities  are  arbitrarily  close  to  the  real  y 
axis  the  contributions  to  the  integral  along  the  segments  of  the  contour 
at  y = ± i°°  are  also  zero.  Therefore,  if  s is  real  and  positive  and  using 
the  notation  of  Figure  2.2(3)  we  have  from  Cauchey's  Integral  Theorum 


loo 


/ V(sy,o,s)eSYXdy  = / Vfsy.o.s^  yAd 

-ioo  Viy 


Finally,  since  on  segment  Im{y)  > 0 and  on  ' Im{y}  < 0,  and  since 

F1  and  r1 1 are  symmetric  mirror  images  of  each  other  across  the  lm{y)  = 0 

axis  we  have  for  any  general  function  of  y,  f(y) 

/ f(y)dy  = / { f (y)+f*(y) 1 dy  = 2i  / Im{ f (y) } dy 

ri+rl'  ri  F1 

where  f*(y)  is  the  complex  conjugate  of  f(y).  With  these  considerations 

we  may  alter  equations  2.2-21  to  read 


1 i°°  i o s(yx-n  h) 

V(x,o,s)  = ^ / Im{i{[2y^neE(s)+2ynangG(s))e 


+ fn3(n32‘T2)E(s)'T(r|g2‘Y2)G(s^e 


s(yx-Hph) 


} Idy 


W(x,0,s)  = / Im{l{  [y(ng2-y2)E(s) 


+ na(np2-y2)G(s) |e 


s(yx-nQh) 


? ....  s(YX'n3h) 


- |2ynrinpE(s)-2y  n(iG(s)  Je 


}}dy 


2.2-22 


Since  our  general  aim  now  is  to  make  these  integrals  look  like  the  for- 
ward Laplace  transform  z,{f(t)f^.  s let  us  define  two  new  variables 

t =-yx  + nh 
a u a 


= - V + n3h 


HHHbLau 
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k 


or  Tc  = - Yc*  + nch  where  c may  be  either  a or  6.  If  we  now  substitute 


for  nc  in  terms  of  y and  c and  solve  for  y we  find 


Tcsin0  ,2.2  2 cose 

r + /r  /c  "Tc  ~F~ 


2 , r 

Tc  -7 


Yc  = 


t s i n 0 o a « A 

-V-  * '^c  -p  ^ 


t. 2 > ^ 2.2-23 

c c2 


And  from  these  expressions,  we  may  derive 


-n. 


dYc  = 2 2 — 

/r  /c-x. 


dx 

2 c 


T‘ 


/7J-'rV  dTc 
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We  now  come  to  the  crux  of  the  problem.  It  is  wished  to  find  a 
path  of  integration  in  the  complex  y plane  such  that  xq  and  x^  will  be 
positive  and  real  for  each  term  of  the  integrand  of  equation  set  2.2-22. 
For  example  let  us  consider  the  term 


1 2y  nBE(s)+2YVBG(s)  "STa 


}dy  2.2-25 


Vl(x’°’s)  = Im{-  R 

The  singularities  of  the  integrand  are  of  prime  importance  since  we  must 
be  able  to  construct  a contour  segment  that  has  the  properties  required 
of  r-|  in  Figure  2.2(3)  and  along  which  xc  (x  in  the  case  of  the  above 
expression)  is  real  and  positive.  Specifically,  this  contour  segment 
must  be  such  as  to  enclose  no  singularities.  In  expressions  such  as  the 
one  above  y = yc  where  c is  the  velocity  subscript  of  xc  in  the  same  ex- 
pression. R,  as  noted  above,  is  the  Rayleigh's  equation  for  this  problem 
which  entered  naturally  when  we  considered  the  free  surface  boundary  con- 
ditions. Thus,  it  is  easily  shown  that  there  are  singularities  at  y = 


i 1/m  y = t 1/R,  and  y = t 1 /cp  where  cr  is  the  solution  to  Rayleigh's 


equation.  The  first  two  sets  of  singularities  are  branch  points  while 


»■*» 
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the  last  set  are  simple  poles.  Since  for  a realizable  elastic  solid 
a > 3 > cr  we  may  draw  the  pertinent  singularities  as  in  Figure  2.2(4). 
The  values  used  in  this  figure  are  a = 5.58  km/sec  and  6 = 2.90  km/sec. 

A branch  line  is  drawn  between  the  two  branch  points  - 1/a  and  - 1/6. 

It  may  be  seen  from  equation  2.2-23  that  for  integrals  involving  xc  such 
as  equation  2.2-25  y is  real  and  negative  for  xc  < r/c  and  will  become 
complex  beyond  this  point.  Notice  also  that  for  very  large  xc  y 
xc/r  [-sin6+icos9] . Thus  y asymptotically  approaches  this  line  in  the 
complex  y plane.  This  path  of  integration  over  the  variable  y satisfies 
all  the  above  conditions  on  and  xc  and  has  the  proper  integration 
limits.  It  is  also  shown  in  Figure  2.2(4)  for  the  above  elastic  wave 
velocities  and  for  x = h = 10  km. 

We  now  wish  to  change  the  variable  of  integration  from  y to  xc  for 
each  term  of  the  equation  set  2.2-22  individually.  Let  us  consider  first 
equation  2.2-25.  This  term  involves  changing  integration  over  y to  that 
over  x^.  The  limits  of  integration  may  be  easily  found.  From  equation 
2.2-23  we  see  that 


lim  x _ rcos0 

a = 

n a 

y-HD 

2 2 1/2 

and  from  the  identities  x = - yx+n  h,  n = [1/a  -y  1 , x = rsine,  and 

a a a ' 1 

h = rcosO  as  well  as  the  asymptotic  approach  of  y to  the  line  x /r[-sin9+ 


icosO]  for  large  x we  have 


L.1" 


lim  x^  = lim  {rcosOf^-p  - {-^-(-sinO+icosO) 

y->_ oof  ioo  T c/r^1”  a*" 


- rsinO[-p(-sin0+icos0)] } 


v"-  fras“i7tI7e,2ei,/2  - sW») 


A 


Fig.  2.2(4)  Contour  for  which  the  application  of  the 

Cagniard-de  Hoop  method  is  allowed.  The  singulari- 
ties of  the  problem  are  also  shown.  In  this  figure 
a=5.58  km/sec,  8=2.90  km/sec,  and  0=45°. 
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= J/JL,  {rcose  -f  el9  - irsinO  y-  e1(J} 
a' 

lim  r i 0 -iG-. 

= . It  p e \ = °° 

t /r+°°  1 or 
a 

where  we  have  used  Euler's  formula 
e10  = cos0  + i sine 

For  0 £ tt/2  we  have  r cos6/a  £ r/a.  This  and  equations  2.2-23  and  2.2-24 
suggest  that  the  integrals  be  separated  into  several  terms  to  avoid  am- 
biguity. Using  the  above  limits  of  integration  and  equation  2.2-24  we 
may  write  equation  2.2-25  as 


r 

a 


Vx'°-s)=ro  ’ 


Im{ 


2y  ngE(s)+2yriangG(s) 


r cose 


‘ST. 


-n  e 

a 


/r2/«2-T  2 

a 


}dT 


+ -i-  / Im{- 

W r 

a 


2y  ngE(s)+2YnangG(s) 
R 


-si 

in  e a 

^V}dTa 

a 


But  for  r cos6/a  < Ta  < r/a  we  have  seen  that  y is  real  and,  furthermore, 
y < 1/a  < 1/6  in  this  range.  Then  nn>  n^*  and  R are  also  real  so  that 
the  integrand  of  the  first  term  in  the  above  expression  is  identically 
zero.  We  are  left  with 


i °°  2y2n  nRE(s)+2yn  2nRG(s)  -st 
V-|  (x,o,s ) = ~~  f Re{ g e a}dT„ 


R/t  2-r2/a2 
a 
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where  we  have  used 

/Imfiz}  dz  - /Re f z ) dz 

We  now  make  one  final  concession  to  generality.  Let  the  functional 
form  of  the  forcing  function  be  the  Dirac  delta  as  previously  defined. 
Then  E(s)  = G(s)  = Mfi(t)}  = 1.  The  problem  has  now  been  maneuvered 


into  the  form  indicated  in  equation  2.2-20.  That  is,  with  the  judicious 
use  of  the  Heaviside  unit  step  function  defined  as 

I"  0 , x < 0 


H(x)  = 


11  , x > 0 


we  have. 


V*-0-s>  'ml 


H(ta- r/a) 

A 2-r2/a2 

a 


2y2Vs 

[Re(  -T-:  ) 


+ Re{- 


2yn „ n„  -st 


'a  '8 


}]e  adx 


R ' J~  ' ‘a 

But  this  is  just  the  definition  of  the  t s forward  Laplace  transform  so 
that  we  have  finally  that 

The  physical  interpretation  of  this  expression  is  straightforward. 
v-|(x,o,t)  is  the  vertical  displacement  at  (x,o,t)  due  to  the  arrival  of 
compressional  waves  from  a point  source  at  (o,h,o). 

There  is  one  additional  complication  encountered  in  taking  dig 
integrals  as  we  shall  now  see. 

Let  us  now  consider  the  second  term  in  the  expression  for  V(x, o,s) 
as  it  appears  in  equation  2.2-22.  This  is 

V2(x,o,s)  =^r|  Im{l(fn3(n62-y2)E(s)  - 


y(n32-Y2)G(s)]e 


s ( yx-riph ) 


) Idy 


2.2-27 


We  may  now  proceed  essentially  as  before.  Now,  however,  the  appropriate 
c in  the  expressions  2.2-23  and  2.2-24  is  8.  The  limits  of  integration 
are  found  to  be 


. 


. . \ 
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1 im  _ r cos6 

Y-0  Tg  6 


. T = °o 

y->-co-f  7 oo  g 

and  we  have  from  equations  2.2-23,  2.2-24,  2.2-27,  and  the  above  consi- 
derations 

r ~st 

1 [nft(nft2-Y2)E(s)-Y(nR2-Y2)G(s) ] - nRe  6 

V?(x,o,s)  = — / Im{ — - — n 2}dTR 

2 ^rcose  R /r2/B2-x  2 B 

~~T~  3 

i oo  nR(nR2-Y2)E(s)-y(nR2-Y2)G(s)  inRe 

+ -77  / Im{ — - — 5 2 V ?7  9>dTft  2.2-21 

r R A 2-r2/g2  6 

6 3 

It  is  easily  shown  that  for  integration  over  the  variable  x^  the  argument 
of  the  integrand  of  the  first  integral  in  the  above  expression  is  not 
real  for  all  x^  between  r cos9/B  and  r/B.  That  is,  the  first  integral 
may  no  longer  be  disregarded.  Since  a > B for  any  elastic  solid  it  fol- 
lows that  1/a  < 1/B.  Then  for  some  6 between  0 and  tt/2  we  must  have 
r cos0/g  < r/a.  Let  us  denote  that  angle  for  which  the  above  inequality 
is  first  true  as  0cr.  For  relative  source-receiver  geometries  such  that 
0 _>  0cr  it  may  be  shown  that  there  is  a time,  t = tcr  at  which  y = - 1/a. 
At  greater  times  y and  ng  remain  real  but  na  becomes  complex.  It  may  be 
further  shown  that 

r cosO  . r 
8 - cr  - g 


and  that 


t = r sin0  + r cos0 


where  cos6cr  = /l-B2/a2 

Thus  equation  2.2-28  may  be  written 
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V2(x’°’s)  = w / lm{" 

C r 


6 nQ(nQ  -Y^)E(s)-y(ne  -Y^)G(s)  -n0  e 


/rV- t/"’6 


i oc  nR(nft2-Y2)E(s)-y(nfi2-Y2)G(s)  inRe 

+ — / Im{— — ■= ^ ^ ^ — p}diR 

vu  r R A 2-r2/32  B 

3 6 

As  before  we  may  use  the  Heaviside  unit  step  function  and  the  i in  the 
second  integral  to  write 

! oo  H( TR-r/i3)  nR2(nfi2-Y2)E(s)-Ynft(nfi2-y2)G(s)  -str 

V2(x,o,s)  = -W  2-Ke(-g— ^ 5 " ~ >e  'JdT 

2 7tMo/T„2-r2/82  R 

P 

, oo  H(r/B-TR)  nc2(nft2-Y2)E(s)-YnR(nR2-Y2)G(s)  -str 

- ><W R ~ 16  H 

Finally,  we  require  e(t)  = g(t)  = 6(t)  so  that  E(s)  = G(s)  = 1.  Then  the 
above  expression  again  becomes  the  definition  of  the  forward  Laplace 
transform  and  we  have 

v2(x,0,t)  . ^^|Re(^,Re(^): 

2 7111  /t  -rZ/B2  R 


cr  /r2/S2-t2 


VV’A  T..rYnB(r,B2'y2). 


The  physical  interpretation  is  again  straightforward.  v2(x,o,t)  is  the 
vertical  displacement  at  (x,o,t)  due  to  the  arrival  of  the  shear  wave 
(first  term)  and  a wave  that  is  reflected  along  the  surface  of  the  half 
space  (second  term).  This  reflected  arrival  leaves  the  source  at  (o,h,o) 
as  a shear  wave  and  is  reflected  as  a congressional  wave.  This  arrival 
can  only  exist  for  those  source-recei ver  geometries  for  which  t £ r/3 
which  we  have  seen  is  equivalent  to  those  geometries  for  which  0 > 0 . 


Thus,  the  physical  significance  of  Qcr  is  also  clear.  For  any  given  a 
and  6,  it  is  the  angle  of  incidence  on  the  half-space  surface  of  a shear 
wave  such  that  a compressi onal  wave  is  reflected,  by  Snell's  Law,  along 
the  surface  of  the  half-space. 

We  may  now  solve  similarly  for  W(x,o,s).  Doing  this  and  separat- 
ing each  solution  into  those  parts  due  to  point  forces  in  the  horizontal 
and  vertical  directions  we  have  the  desired  Green's  functions.  That  is, 
we  derive  four  equations  of  the  form 
G^j  (x,o,t,o,h,o) 

where  this  is  the  displacement  in  the  i direction  at  (x,o,t)  due  to  an 
impulseive  point  force  in  the  j direction  at  the  point  (o,h,o).  These 
four  equations  are  written  below. 


GS  = 1 fH(t-r/a)^  Rc{2y_2VBuH(t-r/e)  p^V^2"^ 
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6 _ 1 ,H(t-r/a) 


zz 


"“[/t2-r2/c2 


2,  2 2v 

(nR  -y  ) 
Ret- a B 


o 2 

2y  n„nc 
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H(t-r/g)  R,tT  a 6 

o o o — R 


/t2-r2/B2 


- h( t-t  ) ^2Y-p-a-n^I  2.2-; 

cr  /r2/8  -t2  K 

These  functions  are  plotted  below  in  Figures  2.2(5)  and  2.2(6)  for  the 
case  x = 5 km,  z = 0.5  km,  a = 5.58  km/sec,  8 = 3.26  km/sec,  and  p = 

2.81  x 10"  dynes/cm  . This  implies  a material  density  of  2.64  gm/cm  . 
This  general  source  geometry  was  chosen  so  that  6 > 0cr  and,  thus,  all 
phases  possible  of  this  medium  are  present.  The  shallow  source  also 
assures  that  the  Rayleigh  wave  pulse  is  well  developed.  The  abscissa 
scale  in  each  of  these  plots  is  1 sec/in.  Following  a criterion  that 
will  be  more  pertinent  for  the  model  accelerograms  to  be  presented  in 
Chapter  4 than  it  is  here,  the  ordinant  scale  is  such  that  the  maximum 
value  of  either  component  of  a given  figure  is  the  zero  to  peak  full 
scale  value  for  the  ordinate  scale  of  that  figure.  Thus,  the  ordinate 
scale  will  generally  vary  from  figure  to  figure  although  it  is  always 
the  same  for  both  components  of  any  given  figure. 

The  arrival  at  t = tcr  of  the  reflected  SP  phase  is  visible,  al- 

* r 

though  small,  on  the  G and  G traces  of  Figure  2.2(5)  as  a small  step 

XX  X z 

in  displacement  approximately  0.1  second  after  the  compressional , or  P- 

2 ? 2 -1/2 

wave,  arrival.  The  [t  -r  /c  ] ' amplitude  decrease  after  the  arrival 

of  a phase  is  also  quite  noticeable,  especially  for  the  G°x  and  G^z  time 
histories.  This  form  of  amplitude  degeneration  is,  of  course,  charac- 
teristic of  the  two-dimensional  problem. 

2.23  The  Spatial  Derivatives  of  the  Green's  Functions 

As  has  been  shown  in  equation  2.1-12,  it  is  the  spatial  deriva- 
tives of  the  Green's  functions  developed  in  the  previous  section  and  not 
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bottom  - for  a = 5.58  km/sec  and  6 = 3.26  km/sec. 
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the  Green's  functions  themselves  that  are  needed  if  we  are  to  calculate 
displacements  due  to  dislocations  on  some  surface  E.  Straightforward 
chain  rule  differentiation  of  the  expressions  in  equation  set  2.2-29  is 
certainly  possible  but  is  also  sufficiently  monotonous  and  difficult  to 
suggest  looking  for  an  alternative  method.  Another  analytic  procedure 
for  evaluating  these  derivatives  does,  in  fact,  exist. 

Let  us  return  to  the  expressions  2.2-22  writing  them  in  conformity 
with  the  notation  of  equation  2.2-29  as 


i i“  i o s(yx-n  h) 
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Because  of  the  invarience  of  the  problem  with  respect  to  translation  in 
the  x direction  it  may  be  shown  that 


A G = tt-  G 2.2-30 

3z  3h 
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where  the  primed  coordinates  are  those  of  the  source.  In  further  devel- 
opment we  will  consider  only  one  of  the  eight  desired  derivations  ex- 
plicitly. The  others  will  follow  virtually  by  inspection.  Let  us 
differentiate  Gxx  with  respect  to  x'.  Then  from  the  previous  expression 
for  Gxx  and  the  relationships  of  equation  2.2-30  we  have 


TTP  0 


s(\x-n  h) 
e a 


2 2 

nR(nR  -y  ) s(Yx-nRh) 

+ ^ e ‘ ]E(s)}d> 

This  differentiation  is  trivial  because  y,  arJ  therefore  na>  n^,  and  R 
are  all  independent  of  x (or  h)  at  this  stage.  The  dependence  of  these 
complex  quantities  upon  x and  h does  not  enter  the  problem  until  we  re- 
quire -t  = y -n  h to  be  such  that  t is  real  and  positive.  Thus, 
c x c c 
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Now,  just  as  before,  we  require  -x  = YX-nch  deriving  expressions  for  y 
and  dY/dxc  identical  to  equations  2.2-23  and  2.2-24,  respectively.  The 
limits  of  integration  are  also  changed  as  before  in  changing  from  y to 
xc  for  the  variable  of  integration.  Then  we  have 
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As  before  the  first  integral  is  zero  and  the  lower  limit  of  integration 

in  the  second  integral  may  be  changed  from  rcosG/3  to  t where  tcr  is 

as  defined  previously.  Using  the  Heaviside  step  function  for  conven- 

'STc 

ience  of  manipulation  and  noting  that  e is  real  and  positive  we  may 
wri  te 
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From  this  point  there  are  several  logical  ways  to  proceed.  Most  funda- 
mentally we  may  recognize  that  s is  real  and  positive  by  definition  and 
may  be  regarded  as  a constant  for  the  purposes  of  integration  over  x . 
Thus,  removing  s from  the  integrals  and  again  requiring  e(t)  = g(t)  = 

6 ( t ) we  have  an  expression  of  the  general  form 

G^x(x,o,s;o,h,o)  = s L{f(t)}  2.2-32 

But  from  the  theory  of  Laplace  transforms 
s L{f ( t) } = f(t)} 
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or,  taking  the  inverse  transform  of  both  sides  of  equation  2.2-32  we 
have 


G$x(x.o.t;o,h,o)  = ^ f(t)  2.2-33 

where,  as  before,  f(t)  is  just  the  sum  of  the  integrands  of  equation 
2.2-31  with  sE(s)  omitted.  Thus, 
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2.2-34 


Alternatively,  we  may  require  e(t)  = g(t)  = H(t)  where  H(t)  is  the  Heavi- 
side unit  step  function.  Then  E(s)  = G(s)  = 1/s  and  sE(s)  = sG(s)  = 1. 
Then 


ax1-  Gxx(x,o,s;o,h,o)  = 


or 

g|r  GjJx(x,o,t;o,h,o)  = f(t)  2.2-35 

U 

Here  Gxx  is  the  displacement  at  (x,o,t)  due  to  a point  force  whose  time 

X Li 

functional  dependency  is  H(t)  at  (o,h,o).  With  G and  G as  defined 
above,  equations  2.2-33  and  2.2-35  are  seen  to  be  special  cases  of  the 
more  general  identity 


Gl]j(x,o,t;o,h,o)  = G*l(x,o,t;o,h,o) 

The  seven  other  desired  quantities  may  be  derived  from  equation  2.2-29 
and  the  above  considerations  virtually  by  inspection.  They  are 
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In  each  case  the  equivalent  expressions  for  3/3 x ' GV,  or  a/a z ' G. . may 

* J 1 J 

be  found  by  deleting  the  differentiation  with  respect  to  time  in  the 
above  equations.  As  we  shall  see  in  the  next  chapter  the  general  ground 
motion  or  instrumental  record  appropriate  for  more  sophisticated  source 
time  histories  may  be  derived  from  the  convolution  of  the  above  Green's 
functions  with  a source  and  instrumental  term.  In  such  a case  it  is 
always  possible  to  modify  the  source  or  instrument  response  terms  in  a 
straightforward  way  so  that  G^  terms  may  be  used  to  obtain  results 
equivalent  to  those  derived  from  the  use  of  G-j  expressions  as  required 


w*V«*  ’ ^ -•  f I . , » -> •*  • • \ 


by  equation  2.1-15.  Thus,  the  differentiation  implied  in  equations 
2.2-34  and  2.2-36  is,  in  practice,  never  necessary. 

The  development  of  these  source  and  instrument  terms  is  the  prin- 
cipal topic  of  the  next  chapter. 
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Chapter  3 

Development  of  Several  Source-Instrument  Terms 

3.1  Introduction 

Having  solved  the  problem  of  an  implusive  line  source  in  a half 
space  we  are  ready  to  apply  this  solution  to  some  practical  question. 
Therefore,  it  is  of  interest  to  find  an  earthquake  whose  ground  motion 
was  recorded  at  a site  approximately  satisfying  the  conditions  under 
which  the  mathematical  solution  was  derived. 

As  we  have  seen  these  requirements  are  fairly  well  met  by  the 
Pacoima  Dam  strong-motion  station  for  the  source  fa ul t- geometry  of  the 
San  Fernando  earthquake  of  February  9,  1971.  However,  before  we  can 
attempt  to  understand  the  strong-motion  records  written  at  Pacoima  Dam 
for  the  San  Fernando  earthquake  we  must  consider  all  the  parameters 
necessary  to  the  generation  of  a model  acceleration  record.  These  will 
include  the  instrument  response  of  the  accelerometer,  the  time  history 
of  the  source  at  any  point  on  the  model  fault,  the  geometry  of  some 
two-dimensional  fault  relative  to  the  location  on  the  surface  of  the 
half  space  at  which  we  wish  to  evaluate  acceleration,  the  propagation 
velocity  of  the  time  history  along  this  fault  surface,  and  the  ampli- 
tude of  the  final  displacements  of  the  sides  of  the  fault  relative  to 
each  other  as  a function  of  position  along  the  fault. 

It  is  the  purpose  of  this  chapter  to  discuss  these  parameters  in 
a general  way.  The  actual  synthesis  of  artificial  accelerograms  will 
be  deferred  to  Chapter  4. 

3.2  The  Strong-Motion  Instruments 

The  mathematical  theory  developed  in  Chapter  2 will  give  the 
ground  displacement  at  some  surface  point  which  has  been  fixed 
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geometrically  relative  to  a line  source  at  some  depth  in  the  half  space. 
Even  when  a collection  of  these  line  sources  have  been  manipulated  so 
as  to  model  the  desired  history  of  motion  on  some  model  fault  surface 
we  must  convert  the  model  motion  at  the  desired  recording  site  into  the 
equivalent  accelerogram  that  would  be  written  by  some  specified  instru- 
ment at  this  site.  This  record  can  then  be  directly  compared  to  the 
actual  strong-motion  record  of  the  test  earthquake.  This  comparison  is, 
of  course,  the  criterion  by  which  the  success  or  failure  of  a model  is 
judged . 

' To  see  how  we  may  include  the  instrument  in  the  synthesis  of  the 
accelerogram  for  a given  model,  let  us  consider  the  theory  of 
forced  and  unforced  single  degree  of  freedom  harmonic  oscillators. 

First  consider  the  damped,  unforced  simple  harmonic  system  repre- 
sented by  the  equation 

x(t)  + 2nx(t)  + mn2x(t)  = 0 3.2-1 

In  this  expression  a dot  stands  for  differentiation  with  respect  to 
time,  ui  is  the  undamped  or  natural  angular  frequency  of  the  oscillator, 
and  n is  a constant  related  to  the  damping  of  the  system.  Where  x(o) 
and  x(o)  are  the  initial  displacement  and  velocity  of  the  system,  res- 
pectively, the  well-known  solution  to  the  differential  equation  3.2-1 
is : 

x(t)  = e”r|tfx(o)  cos(v^>n2-n2  t) 

+ xMtnxiol  sin(/2  2 t)}  3.2.2 

/ c c n 

/a)  -n 

Equations  3.2-1  and  3.2-2  are  useful  for  demonstrating  the  significance 
of  the  two  constants  n and  wn  and  they  are  also  helpful  in  relating  the 
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many  constant  pairs  equivalent  to  n and  w . 

For  example,  let  us  take  advantage  of  the  periodic  nature  of 
equation  3.2-2  and  define  e to  be  the  ratio  of  two  successive  peaks  in 
amplitude  of  x(t).  This  ratio  will  obviously  be  a function  of  the 
degree  of  damping  of  the  system  and  we  see  from  equation  3.2-2  that  the 
period  of  harmonic  motion  is  also  a function  of  damping.  Let  us  call 
the  damped  angular  frequency  and  the  equivalent  damped  Deriod  T^, 


wd  = n ‘n 


2 2 _ 2m 

-n  T 


cos(  vwn  -q  t)  = cos  w^t  = cos 

si n( t)  = sin  u^t  = sin  ^(t+T^) 


and  it  follows  that 
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So  defined,  c is  known  as  the  damping  ratio,  n is  the  damping  coeffi- 
cient. Another  common  way  of  relating  to  wn  is  as  below: 


_ / 2 2 X 2,  2 . X 2 

wd  ~ n _r|  wn/l-ri  / = tan/l-p 


where  obviously  p = n/m  . This  new  constant  is  useful  because  it  is 
the  percent  critical  damping  since  for  p = 1 u> . = 0 and  this  is  the 
criterion  for  critical  damping.  We  may  define  a new  constant,  u,  such 
that  this  is  the  ratio  of  the  damped  to  the  undamped  angular  frequency. 
Thus  we  have 


3.2-5 


■'ff  ■ . , i'.v «.-,» 


60 


We  may  use  these  relationships  to  find  percent  critical  damping,  o,  in 
terms  of  the  damping  ratio,  c.  That  is,  from  above  considerations 

( = e2"P/P  = e2Tip//l-p2 

Taking  the  natural  logarithm  of  both  sides  and  performing  some  algebra 
we  find 
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l_nc 

[ (2tt)2+(  lnc)2j 1/2 


3.2-6 


Using  these  formulas  we  may  find  any  of  the  common  constants  (oun,  w. , 
n,  p,  u,  e)  once  we  know  any  two  of  them.  Since  all  of  these  constants 
are  in  popular  use,  these  relationships  will  prove  useful  later. 

Now  let  us  continue  the  analysis  to  consider  the  forced  case.  The 
equation  of  motion  now  becomes: 


x + 2gx  + = - u 3.2-7 

where  u(t)  is  the  ground  displacement  and  a dot  represents  differentia- 
tion with  respect  to  time,  as  before.  The  response  of  a damped  single 
degree  of  freedom  oscillator  to  general  ground  motion  is  well  known. 
Housner  (1941)  has  presented  the  results  for  an  undamped  oscillator 
using  a method  given  by  Lord  Rayleigh  (Theory  of  Sound,  v.  1_  p.74).  A 
somewhat  more  straightforward  method  using  variation  of  parameters  and 
applicable  to  damped  as  well  as  undamped  systems  is  outlined  by  Bullen 
(1965,  p . 1 45 ) . From  these  earlier  studies  we  may  write: 
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where  w.  and  n are  as  defined  in  equation  3.2-4.  The  asterisk  denotes 


convolution. 
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It  is  often  convenient  to  define  instrument  response  in  terms  of 
ground  velocity  or  ground  displacement.  Integrating  equation  3.2-8  by 
parts  and  writing  the  solution  in  terms  of  a convolution  we  have: 


e‘nt 

x ( t ) = u ( o ) sin  wdt 
d 

+ [ u ( t ) * e'nt(nsin  o>dt  - a>d  cos  ^t)] 


e'nt 
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Finally,  integrating  once  more  by  parts  we  find: 
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wd  dt2  d 

In  these  expressions  u(o)  and  u(o)  are  the  initial  velocity  and  displace- 
ment of  the  forcing  function  respectively. 

Finally,  we  may  take  the  Laplace  transform  of  equation  3.2-7  to 
get  the  frequency  domain  equivalent  of  equation  3.2-10.  This  is 


X(s)  = 
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3.2-11 
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where  X(s)  = x(t)  etc.  This  form  of  the  solution  is  useful  if  u(t) 
is  such  that 


| (s+n)2+wd2l 


U(s) 


has  a recognizable  inverse  transform.  Obviously  where  u(o)  = u(o)  = 0 
the  problem  is  somewhat  simplified. 

3.3  The  Time  History  of  a Point  Source 

As  before  we  will  adopt  the  convention  of  using  the  term  "point 
source"  when  talking  about  an  infinite  line  source  since  this  is  equiva- 
lent to  a two-dimensional  point  source.  Let  us  consider  for  the  moment 
one  such  point  source. 

When  we  first  considered  this  problem  it  was  natural  to  talk  about 
forces  acting  at  this  point.  These  forces  produced  motion  at  the  sur- 
face of  the  half  space  which  we  found  we  could  determine.  The  problem 
was  then  to  find  a way  to  model  motion  at  the  surface  due  to  a disloca- 
tion at  the  same  source  point  because  it  was  desired  to  simulate  sur- 
face motion  due  to  faulting.  The  point  dislocation,  of  course,  has  long 
been  though  of  as  the  fundamental  mathematical  unit  of  fault  model 
motion. 

We  know,  however,  that  this  jump  from  one  position  to  another, 
this  process  whereby  two  originally  coincident  points  find  themselves 
some  finite  distance  apart  with  no  passage  of  time,  is  unlikely  to  be 
an  accurate  description  of  what  happens  on  a fault  surface  during  a 
seismic  event.  Knowing  this  does  not  help  tell  us  what  actually  does 
happen . 

In  this  section  we  will  examine  three  source  terms  that  will  get 
us  from  the  original  to  the  final  state  more  smoothly. 
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3.31  The  Source  as  a Finite  Ramp  Function 

The  simplest  way  to  require  a finite  time  for  finite  displacement 
at  the  source  is  to  use  the  truncated  ramp  function. 

s(t)  = | [H ( t )-H( t-T ) ] + H(t-T)  3.3-1 

as  seen  in  Figure  3.3(1)  motion  begins  at  t = 0,  proceeds  linearly  until 
t = T,  and  then  stops.  T is  called  the  rise  time.  The  amplitude  of  the 
final  displacement  is  normalized  to  one  by  the  form  used  in  equation 
3.3-1.  We  may  easily  find  the  Laplace  transform  of  the  source  term 
equation  3.3-1.  This  is: 

L{ s ( t) } = S(s)  = (l-e‘sT)  3.3-2 

Ts^ 

Now  as  before  let  u(t)  represent  the  ground  motion  at  a given  lo- 
cation and  let  x(t)  be  the  motion  displayed  on  the  record  of  the  instru- 
ment placed  at  the  same  location.  If  now  we  define  i(t)  to  be  the  time 
domain  equivalent  of  the  instrument  transfer  function  we  have: 

x(t)  = i (t)  * u ( t ) 3.3-3 

Let  us  further  resolve  u(t)  into  a geometrical  term  and  a source 
term,  G(t)  and  s(t)  respectively.  For  convenience  we  define  G(t)  to  be 
the  ground  motion  at  some  given  location  due  to  a point  dislocation  that 
turns  on  and  off  like  the  Dirac  delta  function  6 ( t ) . In  short,  G(t)  is 
the  dislocation  Green's  function  for  our  problem.  Another  useful  quan- 
tity is  a similar  G(t)  giving  ground  motion  at  the  same  location  due  to 
a point  dislocation  that  turns  on  like  the  Heaviside  unit  step  function. 
We  will  call  the  first  quantity  G'(t)  and  the  second  GH(t). 

These  two  quantities  are  closely  related  to  each  other  since: 

c/(t)  = ^ GH(t) 


. . \ 
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The  only  reason  to  distinguish  between  them  at  this  point  will  now 
be  made  clear.  In  equation  3.3-3  we  may  now  write 

x ( t ) = i (t)  * IG6( t)  * s ( t ) ] 3.3-4 

The  precise  manner  in  which  our  point  along  the  fault  moves  is  clear  in 
this  expression  since  it  is  just  s(t).  However,  the  more  useful  Green's 
function  to  use  in  this  equation  is  G (t)  since  it  is  this  value  that 
may  be  found  analytically  for  a two-dimensional  dislocation  without  the 
cumbersome  necessity  of  taking  a time  derivative.  Thus,  having  made  the 

{•  U 

distinction  between  G (t)  and  G (t)  we  may  still  keep  a simple  mental 
picture  of  what  the  fault  is  doing  and  write 

x(t)  = i(t)  * GH( t)  * s ( t) ] 
from  which  it  is  a simple  step  to  find: 

x(t)  = i(t)  * [GH( t)  * ^ s(t)]  3.3-5 

for  the  simple  sources  we  will  investigate  here  it  will  always  be 
possible  to  find  the  quantity 

i ( t ) * eft  s ( t ) 3.3-6 

analytically.  Thus,  we  need  only  take  one  numerical  convolution  to 
arrive  at  the  final  model  record  derived  from  a point  dislocation  with 
some  straightforward  source  term. 

The  part  of  equation  3.3-5  represented  by  equation  3.3-6  can  be 
easily  solved.  We  have  already  seen  the  transform  of  i(t). 


/,( i(t)l  = 


(s+n)2+b? 


Also  where  s(t)  is  as  defined  in  equation  3.3-1 


L(dt  s(t) 


Jr  • \ « -f.  * V 
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k 


/2  2 

where  b = w . = /u>  - n and: 


L{ i(t)  * s(t)}  = 


(s+n)2+b2  Ts 


Jr  d-e'sT) 


= T< 


(s+n)2+b2 


0-e"sT)} 


3.3-7 


From  tables  we  may  find  the  inverse  transform  of  this  expression 


to  be: 


i ( t) 


s(t)  = 


,-nt 


dt  bT 

-n(t-T) 


[b  cos  bt-  sin  bt|  H(t) 


bT 


[b  cos  b(t-T)  - nsin  b(t-T)]  H(t-T) 


3.3-8 

A plot  of  this  function  for  T = .25  second  appears  in  Figure  3.3(2).  It 
is  obvious  from  equation  3.3-8  and  from  its  plot  that  discontinuities 
exist  at  t = 0 and  at  t = T.  It  will  be  seen  later  that  these  discon- 
tinuities introduce  severe  difficulties  into  the  problem  of  generating 
a model  accelerogram.  For  the  moment  further  discussion  of  this  prob- 
lem will  be  deferred  until  two  progressi vely  smoother  finite  ramp  func- 
tions can  be  proposed. 

3.32  The  Source  as  a Finite  Ramp  Function  with  Continuous 
First  Derivatives 


Consider  the  source  term: 


s ( t ) = l [l-cos(^)l  (H(t)-H(t-T) ) + H(t-T) 


3.3-9 


We  see  that  displacement  represented  by  this  equation  starts  at  a value 
of  zero  at  t = 0 and  stops  at  a value  of  one  at  t _>  T in  such  a way  that 
this  displacement  and  its  first  derivative  are  continuous  for  all  t > 0. 
A plot  of  this  source  term  with  T = .25  appears  in  Figure  3.3(3). 

As  before,  we  start  further  analysis  of  this  function  by  taking 


» ir'V., ' kr  ' , » ' «- 
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its  Laplace  transform  which  may  be  shown  to  be: 


L{s ( t) } = l(l+e’sT 


)[ 


(VT)2 

s(s2+(tt/T)2 


3.3-10 


Also: 


£<4-  s(t)>  = l(l+e~st)  [-  ^-/D-  -21 
dt  2 s 2+  ( tt/T  ) 2 


and,  as  before, 


2 2 

£{i(t)*^s(t)}  = l(l+e~st)  [ 3.3-11 

dt  2 s2+(tt/T)2  (s+n)2+b2 

Ignoring  for  the  moment  the  translation  factor  e"s^  consider  the 

problem 

l,r— 1.  a2  s2  , 

2[  s2+a2  ’ (s+n)2+b2 


= * L'\ V-9>] 

s +a 


(s+n)2+b2 


where  a = /T.  Now  from  Laplace  transform  tables  we  have: 


1(t)*4  s(t)  = l[L'1{-Pp}*L'1{ 

dt  2 s2+a2  (s+n)2+b2 

= i[a2cosat*[e"ntcosbt  - Pe"ntsinbt]] 


= ifcosat*a2e~n tcosbt)  - (cosat*--r^e~ntsinbt)  ] 


2 | v,  u c v.vj  ju  u i \ uvju  v>  j^“ 

Finally,  from  integration  tables  and  some  algebra  we  have: 


[(n2+(a-b)2)-1 { [bcosat-nsinat]  • 

fn+e"nt{ (a-b)sin(a-b)t  - ncos(a-b)t}] 

+ |bsinat+ncosat| • | (a-b)-e’nt{nsin(a-b)t  + 
(a-b)cos(a-b)tl] 1 
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+ (n2+(a+b)2)_1{[bcosat  + sincat]  • 

[n+e”nt{(a+b)sin(a+b)t-ncos(a+b)t} ] 

+ [bsinat-ncosat]* [(a+b)-e‘nt(nsin(a+b)t  + 

(a+b)cos(a+b)t} ]} ] 3.3-12 

As  may  be  seen  in  Figure  3.3(4),  a plot  of  this  function  for  T = .25  has 
no  discontinuities  at  t = 0 and  t = T but  the  slope  is  discontinuous  at 
these  points.  If  we  include  the  factor  e~s3\  the  form  of  equation  3.3-12 
is  unchanged.  We  have  only  to  substitute  (t-T)  for  T whenever  it  occurs. 
The  total  source-instrument  term  is  then  the  difference  of  these  two 
factors.  It  is,  of  course,  this  difference  that  is  plotted  in  Figure 
3.3(4). 

It  is  now  of  interest  to  carry  the  smoothing  of  the  source  term 
one  step  farther. 

3.33  The  Source  as  a Finite  Ramp  Function  with  Continuous 
First  and  Second  Derivatives 
Consider  the  following  polynomial  source: 

s ( t ) = at5+bt4+ct3+dt2+et+f 

We  wish  to  determint  the  coefficients  by  requiring  that 

s(o)  = ^ s(o)  = s(o)  = ^ s(T)  = -X  s(T)  = 0 
dt  dt 

S (T)  = 1 

It  is  trivial  to  show  that  under  these  conditions 


;(t)  =6t5-l|t4  + ^t3 
r t r 


0 < t < T 


Of  course  s(t)  is  defined  as  being  equal  to  one  for  t > T so  that  we 
may  write 


v.7-  * 
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s(t)  = {[4  t5  - t4  + 14  t3][n(t)  - H(t-T)] 

1 r r 

+ H(t-T)}  3.3-13 

We  see  that  this  equation  represents  a displacement  on  the  fault  much 
like  that  of  the  last  source  term  we  considered.  Now,  however,  no  dis- 
continuities exist  through  the  second  derivative.  A plot  of  equation 
3.3-13  appears  in  Figure  3.3(5),  again  for  T = .25  second. 

We  may  rewrite  equation  3.3-13  in  a form  whose  Laplace  transform 
is  more  easily  evaluated. 

s(t)  = {[4  t5  - '4  t4  + 1°  t3]  H(t) 

r r rJ 

[4  (t-T)5  + (t-T)4  + 4 ( t-T)3]  H(t-T)} 

- jO  7^  TJ 


The  Laplace  transform  of  this  expression  is  simply: 


L{s( t) } = 


60e 


_60_ 
r3 .4 


[ 


12 


TV  TS 

•ST  ]2 


- 1] 


0 2 Ts 
6 


TJ~  ^ 2 ? + Ts  + 


T s 


2 2 

TV 


3.3-14 


Again  we  have  the  sum  of  two  essentially  identical  terms.  As  before, 
one  of  the  terms  has  been  shifted  by  the  e"s^  factor  but  we  know  that 
this  is  of  no  concern.  Here  there  is  also  a different  sign  in  the  sec- 
ond term.  However,  this  too  will  offer  no  particular  difficulty.  Con- 
sidering only  the  first  term  of  equation  3.3-14  we  have: 


[i(t)*ddt  s(t)] 


ifV 


60 


-1 


first  term 
12  6 


(s+n)2+b2 


sT3  l2s2 


Ts 


l]l 


3.3-15 


We  may  immediately  write  the  inverse  of  the  first  transform. 


-1 


60 


1 

(s*  )c+b 


60  “Tit  . . 

e si n bt 
b 


1 Vvi*'  .a*.-.  , 
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The  inverse  of  the  second  transform  is  also  easily  found  to  be: 

1 -l,lr  12  6 . n.  J r6t2  6t  n 

^3  L ' Ts  + 1]}  " * T ‘ 1] 


[1(t)*As(t)]  = sin  bt) 

first  term  T 


* - T + 1)]  3‘3' 

After  several  integration  by  parts  schemes  and  some  algebra  we  find  for 
the  inidividual  terms 


60  r -nx  • k . 60  r.  -ntA 

-sr—  / e sinbxdx  = — * 5 — 5—  {b-e  A 

TJb  0 rb(n  +b  ) 


/ e’rilsinbx-(t-x)dx  = 


T4b(n2+b2) 


+ (b-e'ntA)  - ( +e’ntB) } 

n4+b^  n +b 


360  , . -ni_ _•  1 j.  \2j  _ 360  f.*.2k 

— r — / 6 s i nb  [ • ( t-x ) di  = r r\  o it  b + 

TJb  0 T5b(n  +b2) 

-5^-5-  [te'nt(2nA-bB)  - 3bnt  + 4^4  (n+e'r|tB) 
n2+b2  n +b^ 


+ (b-e“ntA)]} 
where  A = nsinbt  + bcosbt 
B = bsinbt  - ncosbt 

We  may  now  add  these  terms  to  derive  the  solution  to  equation  3.3-15. 

The  second  half  of  the  complete  solution  for  t > T may  easily  be  written 
by  substituting  (t-T)  for  t in  the  above  three  convolutions  and  chang- 
ing the  sign  of  the  second. 
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The  plot  of  the  total  solution  to  s ( t ) is  seen  in  Figure 

3.3(6). 

Note  that  now  there  is  a significant  non-zero  contribution  to  this 
curve  between  the  two  peaks  corresponding  to  t = 0 and  t = T.  We  shall 
see  later  that  this  property  of  this  last  source  term  has  important 
ramifications  when  we  begin  to  generate  model  accelerograms. 

3.34  Comparison  of  the  Source  Terms 

So  far,  as  the  various  sources  have  been  proposed,  we  have  ob- 
served briefly  their  differences  in  the  time  domain.  These  differences 
are  also  meaningfully  studied  in  the  frequency  domain. 

If  the  same  rise  time  is  used  for  each  source  and  if  each  is  nor- 
malized so  that  all  have  an  amplitude  of  one  for  t > T we  may  evaluate 
their  spectra  and  compare  them.  This  is  done  graphically  in  Figure 
3.3(7)  for  T - .25  second. 

It  may  be  seen  immediately  that  the  low  frequency  contribution  of 
all  the  source  terms  falls  off  as  1/s  with  increasing  frequency.  Thus 

| 

for  frequencies  in  our  model  less  than  ~1/T  the  choice  of  source  terms 
is  unimportant.  Indeed,  for  ground  motions  whose  frequencies  are  below 
1/T  the  simple  Heaviside  step  function  would  serve  as  well  as  any  of  the 
more  complex  sources  proposed.  Thus,  if  we  are  attempting  to  match  a 
displacement  record  where  the  dominant  frequency  content  could  very  well 
be  below  1/T  for  some  reasonable  rise  time,  any  source  here  postulated 
will  suffice. 

However,  as  higher  frequencies  are  approached  the  choice  of  source 
term  becomes  significant.  Said  another  way,  if  we  wish  to  distinguish 

' between  these  various  sources  when  modeling  earthquake  ground  motion,  we 

must  look  at  motion  with  high  enouqh  frequency  content  relative  to  the 

f 
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rise  time.  Thus,  if  an  accelerogram  can  be  convincingly  fitted,  there 
is  the  additional  benefit  that  the  form  of  the  source  term  may  be  more 
accurately  determined.  This  could  be  important  for  problems  such  as 
determining  theoretical  near-field  accelerations  which  are  of  great 
engineering  significance. 

3.4  Additional  Parameters 

There  are  three  remaining  parameters  needed  to  define  the  model 
episode  of  faulting.  As  mentioned  in  the  introduction  to  this  chapter 
these  are: 

1)  fault  geometry, 

2)  fault  propagation  velocity,  and 

3)  fault  displacement  amplitude. 

All  three  of  these  parameters  may  be  functions  of  position  along  the 
fault  plane.  The  correct  values  of  these  variables  are  neither  known, 
as  are  the  constants  of  the  Pacoima  Dam  accelerometer,  nor  amenable  by 
arbitrary  postulation  to  a small  number  of  cases,  as  is  the  source  time 
history  term.  However,  we  may  see  how  each  of  these  parameters  of  real 
faulting  may  be  modeled  by  the  use  of  line  dislocation  sources. 

Suppose  that  we  have  selected  appropriate  P-  and  S-wave  velocities 
in  our  half  space.  We  wish  to  find  the  displacement  at  some  surface 
point  u(x,o,t)  due  to  a point  dislocation  at  (x',z',t‘).  The  form  of 
the  Betti-Rayleigh  representation  theorum  needed  for  the  problem  at 
hand  is: 

co 

u|n(s,o,t)  = / / vjf[-i](x,,z,,t')  C^pq  ( x ' , z ' ) 

G.nn(x,°,t;x,,zl,t')  d>.dt ' 3.4-1 

where  is  the  plane  that  will  be  identified  with  the  fault  surface, 
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is  the  unit  vector  perpendicular  to  E pointing  from  the  positive  to  the 
negative  side  of  the  dislocation  on  E.  Here  [5^]  is  defined  as: 

[fM]  - < - q 

For  the  case  of  a thrust  fault  these  quantities  may  be  thought  of  as 
they  appear  in  Figure  3.4(1). 

From  this  figure  we  see  that  v2  = 0.  Let  us  further  assume  that 
the  space  in  which  the  dislocation  exists  is  homogeneous  and  isotropic. 
Then  the  elastic  constants  may  be  written  as 


C. . = A6.  . 6 + u[6 • 6 • +6.  5 . ] 

tjpq  ij  pq  L ip  jq  iq  jpj 

where  A and  p are  Lamp's  constants  and  the  6 ' s are  Kronecker  delta 
functions  as  defined  previously. 

For  m = 1 equation  3.4-2  then  becomes: 


3.4-2 


u,  = 


1 { {vl(’-5l]CHpqGlp,q)  + V3(^l]C13pqGlp,q) 
+v1([C2]C21pqGlp>q)  + v3([^2^C23pqGlp,q) 

+V1  (C^3]C3lpqGlP»q)  + v3^^3^C33pqGlp,q^d5dt' 


3.4-3 


Where  we  have  suppressed  the  functional  notation  and  where,  by 
the  summation  convention  we  have  such  identities  as 

vl(^l^CllpqGlp,q)  = V1^1^C1111G11  ,1  + C1 1 1 2G1 1 ,2 

+ C1 1 1 3G1 1 ,3  + C1121G12,1  + C1122G12,2 

+ C1 1 23G1 2,3  + C11 31G1 3, 1 + C1 1 32G1 3,2 

+ C1 1 33G1 3,3) 

But  for  our  problem  G^  2 = ^1 2 1 = G12  2 = G12  3 = G13  2 = G and  ^rom 
equation  3.4-2 


. • * 
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Cllll  = X+  2m 

cn  1 3 = ci  131  = 0 

C1 133  = x 

Thus 

V^CllpqGlp,q>  = vl[5l]({x+2u)Gllil  + XG13>3) 

and  similarly 

v3(  ^l1C13pqGlp,q)  = ^^1^^11,3^13,1^ 

V[^C31pqGlp,q)  = VltC3^MGn,3+MG13,l) 

v3(^31C33pqGlp,q)  = V^^ll  ,1  + {X+2m)G13,3) 
and  since  [C2 1 = 0 

^1  ( rC2] C21pqGlp,q)  = V3(tC2]C23pqGlp,q^  = 0 
Substituting  these  relations  into  equation  3.4-3  and  deriving  similar 
expressions  for  u3  we  may  write: 

ui(s,o,t)  = 7 / {v,[C1][(A+2y)G11  ^AG^  3] 

Z ’ ’ 

+ v3[^i  1 [T*Gii  ,3+yGi3,l  )1  + V1  ^3^[yGil  ,3+Vj6i3,ll 
+ v3U3]  [AGi]  >1  + (A+2u)Gi3j3]}dZdt'  3.4-4 

Now  let  us  consider  a point  dislocation  oriented  in  some  direction  n. 

Then 

fCj]  = 6(x-x'  )6 ( 2-z ' )6( t-t ' )n 

We  will  call  n the  direction  of  the  fault  plane  at  the  point  ( x 1 , z ' ) - 
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r^i  = [cTjx  + ic3]z  = 

5(x-x' )<5(z-z' )6 ( t-t ' ){sinOx  - cosOz} 

Also  from  figure  3.4(1)  we  see  that 

v-j  [C-j  ] = \ sin20[6(x-x'  )6(z-z*  )6(t-t‘  )1 

V3 [ ] = sin^0[5(x-x,)6(z-z’)6(t-t')] 

v-|[C3]  = - cos^0[6(x-x,)6(z-z,)6(t-t')] 

v3 [^3]  = \ S1* n2e [<5(x-x'  )6(z-z'  )6(t-t'  )J  = - Vjf&j] 

Substituting  these  expressions  into  equation  3.4-4  gives  us  the  dis- 
placement at  (x,o,t)  due  to  a point  dislocation  of  orientation  n at 
( x ’ , z ’ , t ' ) . We  now  have: 

u.j(x,o,t)  = u([sin20Gi1  1 cos20G^  3] 

- {sin28G.3  3 + cos20Gi3  ^]}  3.4-5 

where  G^  ^ = jp-  G^  (x,o,t;x‘  ,z'  ,t' ) 

and  similarly  for  the  other  terms.  We  have  already  calculated  these 
quantities  in  Chapter  2. 

The  notation  of  equation  3.4-5  allows  the  remaining  parameters  of 
the  model  to  be  handled  easily.  That  is,  to  represent  a fault  of  what- 
ever geometry  we  may  string  a series  of  point  sources  together.  Under 
this  scheme  a simple  planar  fault  would  be  represented  by  some  number  of 
point  dislocation  sources  oriented  so  as  to  line  up  on  the  imaginary 
fault  plane. 

Once  the  fault  geometry  has  been  fixed  the  fault  propagation 
velocity  is  easily  modeled  by  delaying  the  origin  time  at  each  point  in 
a prescribed  manner.  Thus,  if  uniform  velocity  is  desired  each  succes- 
sive point  dislocation  source  would  be  made  to  rupture  At  seconds  after 


■ HHr  . . iV>v-T  • 


the  one  before.  The  smaller  the  velocity  wished  to  be  modeled  the 
greater  the  At. 

Finally,  if  we  wish  to  model  varying  dislocation  amplitude  along 
the  fault,  it  is  a simple  matter  to  scale  individual  point  dislocations 
in  a manner  that  will  accomplish  this. 

The  extreme  flexibility  of  this  method  of  modeling  fault  motion 
is  at  once  apparent.  Changes  in  various  fault  parameters  are  easily 
implemented.  In  the  next  chapter  we  will  look  at  some  model  accelero- 
grams based  on  various  parameter  values  suggested  by  many  other 
investigators. 
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Chapter  4 

Comparison  of  Model  and  Actual  Accelerograms 
4.1  Introduction 

In  Chapter  1 were  discussed  several  previous  attempts  to  model  the 
dynamic  ground  motion  in  the  near-field  of  the  San  Fernando  earthquake 
of  February  9,  1971.  It  was  noted  that  several  improvements  in  the 
modeling  procedure  were  now  within  reach.  Primarily,  an  increase  in  the 
sophistication  of  the  medium  transfer  function  was  deemed  desireable  and 
possible.  All  previous  quantitative  near-field  dynamic  studies  have  used 
the  whole-space  transfer  function  of  Maruyama  (1963)  and  Haskell  (1969) 
to  desribe  the  modification  of  the  wave-form  between  the  source  and  the 
point  at  which  the  ground  motion  was  recorded.  This  mathematical  ideali- 
zation precludes  the  explicit  consideration  of  waves  whose  existence  de- 
pends upon  more  complex  medium  characteristics  than  those  of  a perfectly 
elastic,  homogenerous , isotropic  whole-space.  It  also  makes  analysis 
of  model  amplitudes  uncertain.  The  development  of  a half-space  medium 
transfer  function  was  proposed  as  the  most  logical  extension  of  this 
whole-space  theory.  Thus,  in  Chapter  2 the  necessary  ground  work  for 
such  a formulation  was  carried  out  for  the  two-dimensional  case  using 
the  method  of  Cagniard  (1939)  as  modified  by  de  Hoop  (1960). 

Secondly,  as  was  seen  in  Chapter  3,  several  apparently  minor  modi- 
fications of  the  form  of  the  source  time  history  may  have  important  ram- 
ifications to  the  high  frequency  content  of  the  ground  motion.  Since 
the  differences  in  these  time  histories  become  observable  only  at  fre- 
quencies high  relative  to  the  inverse  of  the  rise  time,  representations 
of  ground  motion  that  emphasize  these  high  frequencies  must  be  consider- 
ed. The  accelerograms  written  in  the  near  field  are  ideal  for  this 
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purpose.  In  view  of  these  considerations  it  was  decided  that  the  under- 
standing of  the  seismic  source  could  best  be  increased  by  the  considera- 
tion of  model  accelerograms  written  in  the  near-field  on  the  surface  of 
a half-space. 

The  available  records  most  ideally  suited  for  a useful  comparison 
of  theoretical  to  observed  near-field  strong  motion  within  the  limita- 
tions imposed  by  the  form  of  the  mathematical  development  of  Chapter  2 
are  those  which  were  generated  at  Pacoima  Dam  during  the  San  Fernando 
earthquake.  In  Section  4.2  we  shall  consider  these  records  and  the  form 
into  which  they  must  be  cast  for  the  most  direct  comparison  to  the 
theoretical  records  to  follow. 

Having  at  last  disposed  of  this  preliminary  worlc generation  of 
some  actual  model  accelerograms  may  be  done.  This  is  undertaken  in 
Section  4.3,  first  for  single  point  dislocation  sources  by  the  proper 
combination  of  the  separate  elements  considered  in  previous  chapters, 
and  then  for  a number  of  point  dislocations  taken  together  in  such  a way 
as  to  simulate  propagating  fault  motions  of  various  types.  Accelero- 
grams are  generated  using  both  a planar  fault  whose  geometry  relative 
to  the  Pacoima  Dam  station  is  essentially  that  of  Canitez  and  Toksoz 
(1972)  and  a hinged  fault  similar  to  that  of  Turnbull  and  Battis  (1973). 
A variety  of  fault  rupture  propagation  velocities  and  source  rise  times 
are  used  in  conjunction  with  the  several  forms  of  source  time  history 
explored  in  Chapter  3.  Throughout  Section  4.3  a discussion  is  present- 
ed of  the  implications  of  these  model  accelerograms  as  they  are  com- 
pared to  each  other  and  to  the  Pacoima  Dam  records. 

Finally,  in  Section  4.4  is  presented  a fairly  quantitative  analy- 
sis of  the  uncertainties  in  the  several  model  parameters  chosen  as  best 
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achieving  a satisfactory  match  between  real  and  generated  records  with- 
in the  limits  of  the  method  employed. 

4.2  The  Pacoima  Dam  Accelerograms 
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In  this  section  we  wish  to  consider  the  strong-motion  record 
written  at  Pacoima  Dam  at  the  time  of  the  San  Fernando  earthquake  and 
how  it  may  be  modified  to  allow  comparison  to  model  accelerograms  to  be 
presented  below.  In  Figure  4.2(1)  these  accelerograms  are  presented  in 
their  uncorrected  form  (Brady  et  al , 1971).  As  indicated  in  this  figure 
the  components  shown  are  the  correct  ones  of  S15°W,  N75°W,  and  Down  and 
not  S16°E,  S74°W,  and  Down  as  originally  reported  (Trifunac,  1974).  As 
will  be  the  case  whenever  15  seconds  of  time  history  are  presented  the 
abscissa  scale  is  0.5  sec/div.  The  ordinate  scale,  whose  significance 
will  vary  from  figure  to  figure,  will  be  discussed  in  greater  detail  in 
Section  4.3.  Here  it  is  approximately  1.67  g/in. 

As  previously  indicated,  the  mathematical  development  has  been 
carried  out  only  for  a two-dimensional  line  source  in  a half-space. 

Such  a formulation  may  naturally  be  expressed  in  terms  of  vertical  dis- 
placements and  displacements  in  the  horizontal  plane  perpendicular  to 
the  line  source.  Since  we  wish  to  model  the  faulting  history  of  the 
San  Fernando  earthquake  by  summing  a series  of  such  line  sources  lying 
in  the  plane  of  the  idealized  fault,  theory  dictates  that  the  most 
straightforward  horizontal  direction  to  consider  is  that  perpendicular 
to  the  strike  of  the  fault.  This  strike  is  quite  complicated  in  actu- 
ality but  if  generalized  as  the  average  trend  of  the  fault  zone  it  is 
approximately  N70°W  (Kahle,  1975).  Thus,  we  wish  to  look  at  the  * 
component  of  horizontal  ground  motion  at  Pacoima  Dam  and  since  t- 
not  conform  to  the  actual  disposition  of  the  horizontal  in<-v  -v 
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slight  numerical  resolution  is  necessary.  The  situation  is  summarized 
in  Figure  4.2(2).  Here  the  positive  or  up  directions  are  shown  for  the 
horizontal  components  of  the  Pacoima  Dam  strong-motion  instruments  as 
well  as  for  the  horizontal  component  perpendicular  to  the  general  trend 
of  the  San  Fernando  earthquake  fault  zone.  This  is  labeled  U1  in  con- 
formity with  the  theoretical  development  of  previous  chapters.  As  may 
be  seen,  U-j  and  the  S15°W  Pacoima  component  are  virtually  the  same, 
being  separated  by  an  angle,  6,  of  only  5°.  Specifically 
U^t)  = SI 5°W(t)  cose  + N75°W(t)  sine 

To  carry  out  the  resolution  implied  by  the  above  expression,  it  is 
necessary  to  interpolate  the  S15°W(t)  and  N75°W(t)  time  histories  to 
their  equal  time  interval  equivalents  or  at  least  to  interpolate  one 
time  history  to  the  times  of  the  other.  Such  a process  inevitably  has 
the  effect  of  reducing  the  peak  amplitudes.  The  magnitude  of  this  re- 
duction may  be  accurately  estimated  for  approximately  sinusoidal  wave 
motions  as  a function  of  the  ratio  of  the  interpolation  interval  to  the 
period  of  ground  motion  (Nigam  and  Jennings,  1968)  but  the  problem  be- 
comes more  complex  when  the  ground  motion  is  complex  and  the  interpola- 
tion interval  is  not  discreet  but  varies  from  point  to  point. 

Since  it  was  considered  important  not  to  degrade  the  peak  accelera- 
tion value  at  Pacoima  Dam  it  was  decided  to  use  the  times  of  the  S15°W 
component  for  the  interpolation  of  the  N75°W  component.  In  this  way, 
the  peak  value  of  1.25  g that  occurs  on  the  SI 5°W  component  must  survive 
the  interpolation  process  without  loss  to  round-off  error.  To  estimate 
the  error  introduced  by  this  procedure  another  interpolation  using  the 
N75°W  times  was  carried  out  and  a separate  resolution  done.  The  differ- 


Fig.  4.2(2)  Orientation  of  the  general  trend  of  the  San 
Fernando  fault  zone  relative  to  that  of  the  Pacoima 
• Dam  horizontal  strong-motion  instruments. 
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the  strike  of  the  fault  zone  trend  derived  using  these  two  methods  in- 
dicate an  introduction  of  error  of  less  than  2 percent. 

The  three  components—Down,  U-j  and  N70°W  (parallel  to  and  in  the 
direction  of  the  strike  of  the  fault  trend)— are  shown  in  Figure  4.2(3). 
Interpolation  using  S15°W  times  was  used  in  the  generation  of  this 
figure.  The  abscissa  and  ordinate  scales  are  the  same  as  used  in 
Figure  4.2(1).  The  positions  of  the  horizontal  traces  on  these  plots 
have  been  interchanged.  It  is  the  top  and  bottom  traces  of  Figure  4.2(3) 
that  are  to  be  reproduced  by  some  faulting  process  model.  They  are 
hereinafter  referred  to  as  the  Down  and  U1  components  of  recorded  accel- 
eration at  Pacoima  Dam.  These  appear  in  Figure  4.2(4).  The  format  of 
this  fiqure  will  serve  as  the  model  for  all  subsequent  figures  in  which 
the  model  equivalents  of  the  Down  and  components  are  presented.  The 
abscissa  scale  is  0.5  sec/divasin  Figures  4.2(1)  and  4.2(3).  The  or- 
dinate scale  is  such  that  the  maximum  zero-to-peak  acceleration  is  one 
linear  inch  in  this  figure.  In  this  case,  where  the  maximum  accelera- 
tion is  approximately  1.25  g,  this  implies  an  ordinate  scale  of 
1.25  g/in. 

Several  useful  pieces  of  information  may  be  derived  from  a com- 
parison of  Figures  4.2(1)  and  4.2(3).  First,  we  may  check  the  relative 
importance  of  resolution  and  interpolation  since  the  center  time  history 
of  4.2(3)  which  is  to  be  compared  to  the  bottom  time  history  of  4.2(1) 
has  been  properly  resolved  but  has  undergone  no  practical  interpolation. 
The  remaining  horizontal  record  of  Fiqure  4.2(3)  has  been  both  resolved 
and  interpolated  to  significantly  different  times.  As  we  would  expect 
the  effect  of  resolution  over  such  a small  rotation  is  very  small.  The 
effect  of  interpolation  is  somewhat  greater  but  still  is  not  important 


4.2(4)  Reproduction  of  the  top  and  bottom  traces  of  Fig.  4.2(3). 
These  are  the  records  with  which  all  subsequent  model  accelerograms 
will  be  compared.  The  ordinate  scale  is  approximately  1.25  g/inch. 
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for  the  purposes  of  this  study.  Secondly,  there  is  significant  ground 
motion  in  the  direction  parallel  to  the  fault  trace.  The  implication 
is  that  the  faulting  mechanism  must  have  included  a significant  strike 
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slip  component  which  we  shall  make  no  attempt  to  model  in  this  study. 
This  will  have  no  effect  on  the  model  solution  to  the  remaining  two 
components  presented  in  Figure  4.2(4).  We  may  now  generate  some  model 
accelerograms  which  may  be  directly  compared  to  this  figure. 


Presentation  of  Model  Accelerograms 


All  of  the  elements  needed  to  construct  a series  of  model  accelero- 


grams representative  of  those  generated  on  the  surface  of  a half-space 
due  to  a propagating  fault  have  now  been  considered.  It  remains  but  to 
combine  them  in  the  proper  way. 

In  Section  2.1  and,  more  specifically,  in  Section  3.4  we  saw  how 
to  derive  the  surface  motion  due  to  a point  dislocation  of  arbitrary 
orientation.  The  theory,  which  is  developed  for  point  forces  whose 
time  history  is  that  of  the  Dirac  delta  function,  6(t),  may  be  easily 
generalized  to  include  other  dislocation  time  histories  as  was  shown  in 
Section  3.3.  Ultimately  it  is  the  combination  of  these  point  sources  in 
quite  elementary  configurations  that  is  to  simulate  the  source  of  the 
desired  record.  Thus,  it  is  worthwhile  to  begin  by  considering  the 
record  time  history  of  a single  point  dislocation  for  a series  of  dis- 
location time  histories.  In  this  way  those  effects  due  primarily  to 
the  variation  of  source-time-history  may  be  separated  from  those  due  to 
variation  of  fault  geometry  and  fault  propagation  velocity. 

In  Figure  4.3(1)  are  presented  the  vertical  and  horizontal  dis- 
placements u1(-1.38,0,t)  due  to  a point  dislocation  at  (0,2. 5,0).  The 
time  history  of  this  dislocation  is  that  of  the  Heaviside  unit  step 


Fig.  4.3(1)  Vertical  and  horizontal  displacements  at 

(-1.38,0,t)  due  to  a point  dislocation  at  (0,2. 5,0). 

The  time  history  of  the  dislocation  is  that  of  the 
Heaviside  unit  step  function,  H(t).  For  this  figure 
the  ordinate  scale  is  22.7  cm/division/m  of  dislocation/ 
km  of  fault. 
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function,  H(t).  The  seismic  wave  propagation  velocities  used  are  5.58 
km/sec  and  2.90  km/sec  for  the  compressional  and  shear  waves,  respec- 
tively. This  compressional  wave  velocity  is  in  substantial  agreement 
with  that  employed  in  many  previous  studies  such  as  Hanks,  et  al  (1971), 
Hanks  (1974),  and  Bolt  and  Gopalakrishnan  (1975)  to  name  but  a few. 

The  shear-wave  velocity,  however,  is  significantly  lower  than  has  been 
typically  used.  A more  common  value  would  be  3. 2-3. 3 km/sec  as  expected 
from  the  above  P-wave  velocity  and  the  usual  assumption  that  Poisson's 
ratio  is  0.25.  The  lower  S-wave  velocity  that  has  been  used  for  this 
and  all  subsequent  calculations  is  the  result  of  a desire  to  consider  a 
higher  Poisson's  ratio  for  shallow  crustal  rocks  and  to  empirically 
better  match  the  S-P  interval  of  some  early  model  fault  rupture  schemes 
relative  to  the  observed  S-P  interval  on  the  Pacoima  records.  A study 
of  model  accelerograms  using  shear-wave  velocities  of  both  2.90  and  3.26 
km/sec  indicated  that  the  shapes  of  the  waveforms  were  not  significantly 
affected  by  this  velocity  variation  although  S-phase  arrival  times  ob- 
viously were  affected. 

There  are  several  valid  ways  to  present  the  ordinate.  If  we  in- 
terpret the  point  nature  of  the  source  strictly  so  that  the  actual  geo- 
metrical extent  of  the  source  is  infinitesimal  in  the  two-dimensional 
sense,  the  maximum  displacement  due  to  one  meter  of  dislocation  dis- 
placement at  the  source  may  be  shown  to  be  9.10  x 10"^  cm  at  the  above 
receiver  location.  This  plot  is  constructed  so  that  the  maximum  zero- 
to-peak  amplitude  is  one  inch  or  four  scale  divisions.  In  fact,  all 
succeeding  plots  of  ground,  or  more  commonly,  record  amplitudes  will 
follow  this  convention.  Thus,  on  the  interpretation  that  we  wish  the 
displacement  of  a geometrical ly  infinitesimal  source  the  ordinate  scale 


96 


-A 

may  be  quoted  as  2.27  x 10  cm/division/meter  of  source  dislocation. 
This  number  is  a little  small  for  convenient  quotation  so  another  in- 
terpretation is  worth  considering.  Let  the  source  not  be  infinitesimal 
but  extend  over  a finite  width  of  dislocation  surface  of  the  same  orien- 
tation. Then  the  effect  of  a dislocation  of  equal  amplitude  occurring 
over  a strip  of  dislocation  surface  may  be  approximately  modeled  by 
multiplying  the  above  amplitude  by  the  width  of  the  strip.  Under  this 
approximation,  the  ordinate  scale  of  Figure  4.3(1)  may  be  quoted  as 
2.27  x 10^  cm/division/meter  of  dislocation/km  of  fault  width.  We  have 
referred  to  the  width  of  a fault  strip  rather  than  to  a unit  fault  area 
since  this  is  a two-dimensional  solution. 

The  most  noteworthy  aspect  of  the  ground  motion  shown  in  Figure 
4.3(1)  is  its  unimpressive  appearance.  The  only  large  amplitude  arrival 
is  that  of  the  S-wave  although  the  P-wave  arrival  is  discernible.  For 
this  source-receiver  geometry  the  Rayleigh-wave  phase  and  any  permanent 
ground  offset  are  not  significant  relative  to  the  body-phase  amplitudes. 

We  are  now  in  a position  to  investigate  the  effect  of  the  convolu- 
tion of  this  ground  motion  with  various  source- instrument  terms.  As 
previously  shown,  the  source-instrument  part  of  this  convolution  may  be 
done  analytically  for  a number  of  source  time  histories.  Basically 
then,  we  shall  be  numerically  convolving  curves  such  as  those  shown  in 
Figures  3.3(2),  3.3(4),  and  3.3(6)  with  ground  displacements  such  as 
those  shown  in  Figure  4.3(1)  to  generate  the  instrumental  records  of  an 
accelerometer  with  specified  instrument  constants  at  the  specified  re- 
ceiver position  due  to  a dislocation  of  specified  time  history  and 
orientation  at  the  specified  source  position.  This  is  the  only  step  in 
which  a numerical  rather  than  an  analytical  operation  must  be  performed 


97 


to  derive  model  accelerograms  from  a two-dimensional  point  source  on 
the  surface  of  an  elastic  half-space. 

The  source  time  histories  considered  are  various  forms  of  the 
finite  ramp  dislocation  displacement  function  starting  at  zero  amplitude 
at  time  t = 0 and  stopping  at  unit  amplitude  at  t = T,  where  T is  the 
rise  time.  The  differences  between  the  three  source  time  histories  of 
Section  3.3  may  be  most  meaningfully  stated  in  terms  of  the  smoothness 
with  which  they  start  and  stop.  The  most  elementary  finite  ramp  is 
continuous  for  all  time  but  has  discontinuous  first  and  second  deriva- 
tives at  t = 0 and  t = T.  The  second  finite  ramp  source  is  continuous 
for  all  time  through  its  first  derivative  but  has  discontinuities  in 
its  second  derivative  at  the  usual  points,  t = 0 and  t = T.  Extending 
this  progression  logically,  the  final  finite  ramp  source  is  continuous 
for  all  time  through  its  second  derivatives.  For  brevity,  these  source 
time  histories  will  be  referred  to  hereinafter  as  sources  one,  two,  and 
three,  respectively. 

The  instrument  constants  used  are  those  of  an  . ccel erometer  with 
n = 68.49  and  wn  = 120.86  radians/sec,  where  n is  the  damping  coeffi- 
cient and  u)n  the  natural,  undamped  angular  frequency.  The  static  mag- 
nification used  is  113.2  which  corresponds  to  an  acceleration  sensi- 
tivity of  7.59  cm/g  for  the  above  value  of  go  . All  these  constants 

n 

correspond  to  the  average  value  of  the  equivalent  constants  of  the  three 
Pacoima  Dam  strong-motion  instruments. 

Experience  has  shown  that  a meaningfully  complete  set  of  record 
time  histories  may  be  generated  by  considering  rise  times  of  0.1,  0.25, 
and  0.5  seconds  for  each  of  the  three  source  time  histories  considered 
above.  Thus,  nine  plots  are  necessary  to  adequately  investigate  the 
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effect  of  source  time  history  variation  on  the  accelerometer  records 
pertinent  to  a given  source-receiver  geometry  and  dislocation  orienta- 
tion. For  the  specific  source-receiver  geometry  and  orientation  of  the 
ground  motion  of  Figure  4.3(1)  these  nine  accelerograms  are  presented 
in  Figures  4.3(2)  through  4.3(10).  Figures  4.3(2)  through  4.3(4)  are 
the  accelerograms  using  source  one.  Figures  4.3(5)  through  4.3(7)  are 
the  accelerograms  using  source  two,  and  Figures  4.3(8)  through  4 . 3 ( 1 0 ) 
are  those  using  source  three. 

Several  important  observations  may  now  be  made  from  a comparison 
of  these  three  sets  of  figures.  The  single  ground  displacement  pulses 
at  t = r/3  and,  to  a lesser  extent,  at  t = r/a  in  Figure  4.3(1)  have 
become  double  pulses  separated  by  times  equivalent  to  the  appropriate 
source  rise  time.  These  two  pulses  may  be  thought  of  as  the  starting 
and  stopping  phases  associated  with  the  dislocation  process  at  the 
source  point.  This  pulse  interval  is  obviously  a sensitive  barometer 
of  source  time  history  rise  time  for  all  sources  considered. 

In  contrast,  the  general  shapes  of  the  individual  pulses  are 
greatly  different  for  different  source  terms.  The  frequency  content  of 
the  starting  and  stopping  phases  of  source  one  is  apparently  unaffected 
by  the  rise  time.  For  source  two  the  impulsive  nature  of  these  two 
phases  is  still  apparent  but  more  rise  time  dependent  motion  is  seen  to 
occur  between  t = r/c  and  t - r/c  + T,  where  c is  either  the  congres- 
sional or  shear  wave  velocity.  For  source  three  the  pulse  width,  as 
well  as  the  pulse  separation,  is  closely  related  to  the  rise  time. 

Finally,  as  shown  in  Figure  3.3(7),  the  less  continuous  the  source 
the  more  high-frequency  energy  it  contains.  Therefore,  we  would  expect 
higher  amplitude  accelerations  to  be  associated  with  the  less  continuous 
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Fig.  4.3(2)  Model  accelerogram  at  (-1.38,0,t)  due  to 
a point  dislocation  at  (0,2. 5,0).  The  time  history 
of  the  dislocation  is  that  of  a finite  ramp  (source 
one)  with  a rise  time  of  0.1  sec.  The  ordinate 
scale  is  18.53  g/division/m  of  dislocation/km  of 
fault. 
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Fig.  4.3(3)  Model  accelerogram  generated  under 
conditions  identical  to  those  of  Fig.  4.3(2) 
except  now  the  rise  time  is  0.25  sec.  The 
ordinate  scale  is  7.36  g/division/m  of  dis- 
location/km of  faul t. 


Fig.  4.3(4)  Model  accelerogram  generated  under 
conditions  identical  to  those  of  Fig.  4.3(2) 
except  now  the  rise  time  is  0.5  sec.  The 
ordinate  scale  is  3.6ff  g/division/m  of  dis- 
location/km of  faul t. 


■ • % 


* Pig-  4.3(5)  Model  accelerogram  at  ( -1 . 38 ,0 , t ) due 

to  a point  dislocation  at  (0,2. 5,0).  The  time 
history  of  the  dislocation  is  that  of  a finite 
ramp  with  continuous  first  derivatives  (source 
two).  The  rise  time  is  0.1  sec  and  the  ordi- 
nate scale  is  7.40  g/division/m  of  dislocation/ 
km  of  fault. 
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Fig.  4.3(6)  Model  accelerogram  generated  under 
conditions  identical  to  those  of  Fig.  4.3(5) 
except  now  the  rise  time  is  0.25  sec.  The 
ordinate  scale  is  1.23  g/division/m  of  dis- 
location/km of  faul t. 


Fig.  4.3(7)  Model  accelerogram  generated  under 
conditions  identical  to  those  of  Fig.  d.3(5) 
except  now  the  rise  time  is  0.5  sec.  The 
ordinante  scale  is  0.30  g/division/m  of  dis- 
location/km of  faul t. 
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Fig.  4.3(8)  Model  accelerogram  at  (-1.38,0,t)  due 
to  a point  dislocation  at  (0,2. 5,0).  The  time 
history  of  the  dislocation  is  that  of  a finite 
ramp  with  continuous  first  and  second  deriva- 
tives (source  three).  The  rise  time  is  0.1  sec 
and  the  ordinate  scale  is  9.68  g/division/m  of 
dislocation/km  of  fault. 
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Fig.  4.3(10)  Model  accelerogram  generated  under 
conditions  identical  to  those  of  Fig.  4.3(8) 
except  now  the  rise  time  is  0.5  sec.  The 
ordinate  scale  is  0.19  g/division/m  of  dis- 
location/km of  fault. 
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sources,  all  other  things  being  equal  such  as  rise  time  and  absolute 
dislocation  amplitude.  Similarly,  since  the  corner  frequency  is  in- 
versely proportional  to  the  rise  time,  T,  for  a given  source,  as  was 
seen  in  Chapter  3,  more  high-frequency  energy  is  present  at  smaller 
rise  times.  Both  of  these  rather  obvious  contentions  are  amply  sup- 
ported by  the  maximum  accelerations  given  in  Figures  4.3(2)  through 
4.3(10).  Again  normalizing  to  those  accelerograms  given  for  one  meter 
of  dislocation  over  a one  kilometer  strip  of  fault  we  find  peak  accelera 
tions  of  14.7  g,  1.2  g,  and  .75  g for  sources  one,  two,  and  three, 
respectively,  using  a rise  time  of  T = 0.5  sec.  And,  normalizing  in  an 
identical  way,  we  find  peak  accelerations  of  38.6  g,  4.3  g,  and  again 
.75  g for  source  three  using  rise  times  of  0.1  sec,  0.25  sec,  and  0.5 
sec,  respectively.  It  is  clear  from  these  figures  that  very  large  peak 
accelerations  are  deterministically  obtainable  from  the  formulation  of 
the  problem  employed  here  if  no  mechanical  constraint  imposed  by  the 
allowable  behavior  of  real  elastic  material  is  appealed  to.  It  is  also 
clear  that,  since  the  waveforms  and  their  amplitudes  may  be  derived  in- 
dependently by  this  method,  peak  acceleration  amplitude  will  furnish  a 
valuable  cross  check  on  the  validity  of  any  final  solution.  A model 
accelerogram  that  matches  the  Pacoima  Dam  resolved  records  of  Figure 
4.2(4)  with  respect  to  general  appearance  of  the  time  history  must  also 
give  reasonable  peak  acceleration  values  for  reasonable  dislocation 
amplitudes  on  the  model  fault  surface. 

We  may  now  summarize  the  effect  of  varying  the  source  dislocation 
time  history.  More  continuous  sources  of  greater  rise  time  lead  to 
accelerograms  of  lower  frequency  content  and  lower  amplitude  assuming 
that  the  dislocation  amplitude  and  source-receiver  geometry  remain  fixed 
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As  will  be  shown  in  the  following  sections  these  conclusions  are  also 
applicable  to  the  sum  of  a series  of  such  point  sources  combined  in 
such  a manner  as  to  simulate  propagation  of  a dislocation  over  a fault 
surface  of  arbitrary  configuration. 

4.31  Model  Accelerograms  from  a Planar  Fault 

With  three  source  functions,  each  specified  for  three  rise  times, 
it  is  easy  to  see  that  a completely  unmanagable  number  of  model  accel- 
erograms may  be  generated  from  a relatively  few  dislocation  schemes. 
Practically,  we  can  consider  only  a very  few  fault  model  geometries 
since  the  source-receiver  geometry  is  the  most  fundamental  parameter 
that  must  be  considered.  The  evaluation  of  equation  3.4-5  for  both 
components  of  ground  displacement  at  two  thousand  time  points  requires 
approximately  12  seconds  on  the  CDC  6400  computer.  This  number  of  time 
points  allows  generation  of  10  seconds  of  ground  motion  time  history  at 
a time  increment  of  0.005  seconds.  For  the  two  basic  fault  models  to 

be  considered  below  two  hundred  line  sources  will  be  properly  arrayed 

j 

to  simulate  a continuous  fault  of  the  specified  geometry.  Thus,  the 
computation  of  what  may  be  termed  the  geometrical  part  of  the  solution 
involves  computer  evaluation  of  equation  3.4-5  eight  hundred  thousand 
times  at  an  expenditure  of  almost  forty  minutes  of  CDC  6400  computer 
- time.  By  contrast,  computation  that  allows  numerical  convolution  of 

this  geometrical  term  with  a given  source-time-history/instrument  term 
and  which  simulates  fault  propagation  velocity  variation  requires  gen- 
erally less  than  one  minute.  It  is  clear  that  fault  geometry  is  not 
well  suited  to  solution  by  random  experimentation. 

' Fortunately,  fault  geometry  may  be  inferred  from  many  types  of 

j data  and  consequently  is  one  of  the  better  determined  of  the  problem 
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parameters.  As  was  seen  in  Chapter  1 most  early  studies  found  that  the 
data  was  adequately  explained  by  a planar  fault  dipping  at  about  50°. 
Specifically,  by  analysis  of  the  polarity  of  teleseismic  P-  and  S-waves, 
of  the  spectra  of  surface  waves,  and  of  the  static  displacements, 
strains,  and  tilts  in  the  source  region,  Canitez  and  Toksoz  (1972)  de- 
rived a model  fault  geometry  similar  in  all  essentials  to  that  shown  in 
Figure  4.3(11).  Also  shown  is  the  position  of  the  surface  point  at 
which  acceleration  records  are  to  be  calculated.  Since  this  point  is 
to  represent  the  position  of  Pacoima  Dam  relative  to  the  real  fault 
surface  it  is  so  labeled  in  the  figure.  The  distance  between  the  zone 
of  surface  faulting  of  the  San  Fernando  earthquake  and  the  strong-motion 
instrument  site  at  Pacoima  Dam  is  an  easily  determinable  quantity  since 
it  may  be  measured  directly.  From  Kamb,  et  al  (1971)  this  distance  is 
slightly  in  excess  of  5 kilometers.  The  equivalent  distance  in  Figure 
4.3(11)  is  a little  greater  than  this  because  it  was  originally  defined 
to  be  one-half  the  distance  between  the  epicenter  and  surface  fault 
rupture  of  a planar  fault  originating  at  a depth  of  14  kilometers  and 
dipping  at  52°  to  the  horizontal.  Defined  in  this  way  the  pertinent 
distance  becomes  5.47  kilometers. 

To  model  the  ground  motion  at  Pacoima  Dam  due  to  a motion  along 
the  entire  fault,  the  ground  motion  from  individual  sources  along  the 
fault  plane  are  calculated  and  summed.  Thus,  the  next  step  is  to  carry 
out  these  calculations  and  store  them  for  later  summation.  As  mention- 
ed above  it  was  decided  somewhat  arbitrarily  to  model  the  fault  using 
two  hundred  two-dimensional  point  sources  evenly  distributed  over  its 
surface  because  it  was  found  empirically  that  this  spacing  of  individual 
sources  allowed  the  generation  of  model  accelerograms  whose  appearance 
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Fig.  4.3(11)  Planar  model  of  the  fault  associated  with  the 
San  Fernando  earthquake  of  February  9,  1971  (after 
Canitez  and  Toksoz,  1972).  The  hypocenter  of  the  earth 
quake  is  at  the  lower  end  of  the  fault.  The  sense  of 
fault  motion,  the  relative  position  of  Pacoima  Dam,  and 
the  seismic  velocities  used  in  subsequent  accelerograms 
are  also  shown. 
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did  not  depend  upon  any  reasonable  change  in  the  interval  between 
sources.  These  accelerograms  could,  thus,  be  considered  more  a property 
of  the  fault  model  than  of  the  particular  method  of  evaluation.  The 
ground  motion  at  the  receiver  is  then  stored  on  magnetic  tape  for  each 
source  point  and  for  both  components  so  that  the  numerical  equivalent 
of  two  hundred  figures  similar  to  Figure  4.3(1)  are  now  available. 

The  summation  of  these  two  hundred  ground  motions  is  now  desired 
in  such  a way  as  to  simulate  propagation  of  the  dislocation  episode 
along  the  fault.  It  is  obvious  that  this  may  be  done  quite  simply  for 
certain  discrete  values  of  fault  propagation  velocity.  For  example,  let 
the  fault  length  be  L and  the  incremental  length  between  point  sources 
along  the  fault  be  A£.  From  the  above  requirement  that  there  be  two 
hundred  of  these  point  sources  we  see  that  = L/200.  Further,  let  the 


displacement  at  Pacoima  Dam  at  time  t due  to  source  j be  U.(t)  where  U 

J 

may  be  either  the  horizontal  or  vertical  component  of  ground  motion. 

Now  assume  that  the  (j+l)t*1  source  begins  at  At  seconds  later  where  At 
is  the  time  increment  at  which  ground  and  record  motions  have  been  cal- 
culated. From  above  At  = 0.005  seconds.  The  combined  ground  motion  at 


the  receiver  is  then  Uj(t+At)  = Uj+, 


(t). 


This  process  effectively 


models  a fault  propagation  velocity  of  Ail/At  or,  for  the  geometry  of 
Figure  4.3(11),  a velocity  of  17.77  km/sec.  This  is  the  largest  finite 


velocity  obtainable  from  the  above  incremental  process  which  may  be 


represented  by  the  following  series  of  equations: 

U(0)  = 11,(0) 


U(At)  = U, (At)  + U2(0) 

U(?At)  = (J,  (2At)  + U 2 ( A t ) + U3(0) 

U(3At)  = U,(3At)  + U 2 ( 2At ) + U3( At)  + U4(0)  4.3-1 


In  these  expressions  the  time  incremental  nature  of  the  evaluation 
of  displacement  time  history  has  been  made  explicit.  Slower  fault 
propagation  velocities  may  be  modeled  by  allowing  multiple  values  of  the 
time  increment  between  source  points.  In  this  case  the  fault  propaga- 
tion velocity  is  A£/nAt,  n = 2,3,4 For  example,  with  n = 2 we  have 

a propagation  velocity  of  8.88  km/sec  and,  in  analogy  to  equation  set 
4.3-1 

U(0)  = U}(0) 

U(At)  = Uj ( At) 

U(2At)  = U](2At)  + U2(0) 

U(3At)  = U] (3At)  + U2(At) 

U(4At)  = U](4At)  + U2(2At)  + U3(0)  4.3-2 

The  method  now  is  to  retrieve  the  individual  source  time  histories 
from  the  magnetic  tape  and  sum  them  as  indicated  in  equation  sets  4.3-1 
or  4.3-2  or  their  equivalent  for  some  other  value  of  n.  This  sum  now 
represents  the  ground  displacement  at  Pacoima  Dam  due  to  two  hundred 
point  sources  of  infinitesimal  geometric  extent  evenly  distributed  over 
the  fault  surface  of  Figure  4.3(11)  and  whose  time  history  of  dislocation 
is  that  of  a Heaviside  unit  step  function  propagating  over  the  fault 
surface  at  a velocity  of  AH/nAt.  It  remains  but  to  convolve  this  summed 
ground  motion  with  the  various  source-time-history/instrument  terms  dis- 
cussed in  Chapter  3 and  used  to  generate  accelerograms  4.3(2)  through 
4.3(10).  Since  convolution  of  the  sum  of  a series  of  functions  is  dis- 
tributive we  need  only  convolve  the  above  sum  with  an  appropriate  source- 
instrument  term  to  derive  an  accelerogram  identical  to  that  generated  by 
the  summation  of  individual  point  source  accelerograms.  This,  of  course, 
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saves  a great  deal  of  time  and  expense. 

Previous  dynamic  studies  (Mikuino,  1973;  Trifunac,  1974;  Bolt, 

1972;  Hanks,  1974;  Turnbull  and  Battis,  1973)  have  generally  indicated 
fault  rupture  propagation  velocities  of  between  2.0  and  3.0  km/sec. 
Therefore,  it  was  decided  to  put  most  emphasis  on  this  range  of  values. 
In  fact,  preliminary  efforts  of  this  study  were  suggestive  of  velocities 
in  the  lower  half  of  this  range.  Thus,  three  values  were  used.  These 
were  2.538  km/sec,  2.221  km/sec,  and  1.974  km/sec  as  dictated  by  the 
incremental  considerations  above. 

An  exhaustive  use  of  all  the  variables  considered  would  entail  the 
generation  of  twenty-seven  model  accelerograms  for  each  model  fault 
geometry.  This  was,  in  fact,  not  necessary  since  it  was  discovered 
quite  early  that  the  effect  of  fault  propagation  velocity  was  not  sig- 
nificant in  altering  those  features  of  an  accelerogram  which  might  be 
characterized  as  general  frequency  content.  Thus,  if  a model  accelero- 
gram did  not  "look  like"  the  Pacoima  record  because  of  a difference  in 
predominant  period  of  significant  phases  the  alteration  of  the  propaga- 
tion velocity  would,  in  general,  not  improve  or  degrade  the  match.  This 
will  become  clear  when  later  figures  are  compared  representing  model 
accelerograms  whose  sources  differ  only  in  propagation  velocity. 

We  are  now  in  a position  to  present  some  model  accelerograms.  Be- 
ginning with  the  most  simple  cases  in  Figures  4.3(12)  through  4.3(14) 
are  presented  the  model  accelerograms  for  the  Cani tez-Toksoz  fault  ge- 
ometry, a truncated  ramp  source-instrument  term  (type  one  source)  and  a 
fault  rupture  propagation  velocity  of  1.974  km/sec.  The  source  rise 
times  are  0.1,  0.25,  and  0.5  seconds,  respectively.  The  significance 
of  the  ordinate  scale  is  similar  to  that  of  Figures  4.3(2)  through 


118 


4.3(10).  Each  point  source  is  thought  of  as  acting  over  a segment  of 
the  fault  whose  width  is  ML.  With  this  assumption  we  may  quote  the 
ordinate  scale  in  terms  of  centimeters  of  record  displacement/division/ 
meter  of  dislocation  along  the  entire  fault.  We  have,  of  course,  as  yet 
not  considered  the  effect  of  allowing  dislocation  amplitude  to  vary  with 
position  on  the  fault.  Alternatively,  since  the  static  magnification  of 
the  instrument  used  to  derive  expressions  for  the  source-instrument  term 
is  such  that  the  acceleration  sensitivity  is  7.59  centimeters  of  record 
displacement/g  of  ground  acceleration  (where  g is,  as  usual,  the  accel- 
eration of  gravity)  we  may  express  the  ordinate  scale  in  terms  of  g's/ 
division/meter  of  dislocation  along  the  fault.  It  is  this  second  con- 
vention that  will  be  followed  throughout  the  remainder  of  this  study. 
Thus,  the  ordinate  scales  of  Figures  4.3(12)  through  4.3(14)  are  8.81, 
3.42,  and  1.73  g's/division/meter  of  fault  displacement,  respectively. 

The  similarity  of  these  first  three  model  accelerograms  to  the 
Pacoima  records  as  presented  in  Figure  4.2(4)  is  disappointing  in  many 
respects.  First  as  might  have  been  anticipated  from  the  single  point 
source  accelerograms  of  Figures  4.3(2)  through  4.3(4)  which  employed  the 
same  source- instrument  term,  the  frequency  content  of  these  model  accel- 
erograms is  much  too  high  compared  to  that  of  the  Pacoima  records.  The 
pulse  width  of  each  individual  spike  of  acceleration  is  simply  not  ade- 
quate to  allow  constructive  or  destructive  interference  to  be  noticeable 
when  the  motions  of  two  hundred  sources  are  superposed  in  the  prescribed 
manner.  Secondly,  the  constraint  imposed  on  the  magnitude  of  fault  dis- 
location by  these  accelerograms  is  unacceptable.  If  these  peak  values 
of  acceleration  are  to  match  the  approximately  1.25  g value  of  the 
Pacoima  maximum  acceleration,  fault  dislocations  over  the  entire 
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dislocation  surface  of  0.04,  0.09,  and  0.18  meters  are  implied  for  rise 
times  of  0.1,  0.25,  and  0.5  sec,  respectively.  These  derived  values  of 
average  dislocation  amplitude  are  far  too  small  as  may  be  seen  from 
previously  derived  estimates  for  this  quantity  as  presented  in  Table 
1-1.  Finally,  it  is  obvious  that  some  long  period  noise  has  entered 
the  calculation  during  the  numerical  convolution  of  the  source- 
instrument  term  with  the  ground  displacement  term.  It  is  intuitively 
unacceptable  to  have  a model  accelerogram  with  a final  constant  step  in 
acceleration.  After  the  last  wave  from  the  faulting  process  has  passed 
the  receiver  position  all  acceleration  must  return  to  zero. 

In  spite  of  these  shortcomings  the  accelerograms  of  Figures  4.3(12) 
through  4.3(14)  are  not  completely  devoid  of  interest.  It  can  be  seen, 
for  example,  that  using  the  1.974  km/sec  fault  propagation  velocity  of 
these  particular  model  accelerograms  a total  record  of  approximately 
8.5  sec  is  written  between  the  arrival  of  the  first  P-wave  at  2.7  sec 
and  that  of  the  last  surface  wave  at  about  11.2  sec.  Thus  using  similar 
fault  propagation  velocities  and  fault  model  geometries  we  may  have 
reasonable  hopes  of  modeling  about  8 sec  of  the  actual  Pacoima  record 
without  recourse  to  more  complicated  medium  transfer  functions  or  dis- 
location time  history.  It  may  also  be  seen  that  surface  waves  are  quite 
apparent  on  all  records  and  that  their  peak  amplitudes  are  about  as 
large  as  those  of  earlier  body-wave  phases.  The  importance  of  surface 
waves  to  near-field  accelerograms  has  often  been  discounted  since  their 
prominence  at  teleseismic  distances  is  due  to  their  less  rapid  geometric 
attenuation  relative  to  body  waves.  Here  we  see,  however,  that  where 
the  fault  rupture  actually  approaches  the  surface  of  the  earth  these 
surface  waves  are  likely  to  play  an  important  role  in  near-field  ground 


motions.  One  other  feature  of  these  very  preliminary  accelerograms  is 
worthy  of  note.  That  is,  the  duplication  of  pulses  at  t = r/c  and  t = 
r/c  + T observed  before  in  connection  with  the  single  point  source 
accelerograms  has  survived  the  superposition  process.  This  is  particu- 
larly evident  for  the  surface  wave  train  near  the  end  of  these  records 
but  may  also  be  seen  for  earlier  body-wave  pulse  pairs.  This  phenomenon 
may  be  very  useful  in  determining  the  proper  rise  time  to  use  if  suit- 
able distinct  pulse  pairs  may  be  observed  in  the  real  record. 

We  have  yet  to  explore  the  effect  of  variation  of  the  propagation 
velocity  of  the  dislocation  along  the  fault.  In  view  of  the  disagree- 
able aspects  of  the  accelerograms  in  Figures  4.3(12)  through  4.3(14)  we 
will  defer  this  more  complete  exploration  of  all  the  parameters  avail- 
able to  a given  geometry  until  a more  satisfactory  source-instrument 
term  is  employed.  We  now  undertake  this  task  using  source  two. 

Figures  4.3(15)  through  4.3(23)  represent  all  possible  combina- 
tions of  the  three  rise  times  T = 0.1,  0.25,  and  0.5  sec  with  the  three 
fault  rupture  propagation  velocities  2.538,  2.221,  and  1.974  km/sec 
using  source  two  and  the  planar  fault  geometry  of  Figure  4.3(11).  These 
particular  velocities  are  dictated  by  the  discrete  nature  of  the  problem 
formulation  as  we  have  already  seen.  Let  us  begin  by  considering  the 
three  model  accelerograms  associated  with  the  fastest  propagation  ve- 
locity-^.538  km/sec.  These  appear  in  Figures  4.3(15)  through  4.3(17) 
for  rise  times  of  0.1,  0.25,  and  0.5  sec,  respectively. 

Many  of  the  features  evident  in  these  accelerograms  are  the  logi- 
cal extensions  of  their  source  one  counterparts  as  presented  in  the 
previous  three  figures.  For  example,  the  S-P  interval  of  approximately 
2.5  sec  has  not  changed  since  this  is  a function  only  of  the  hypocenter- 
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Fig.  4.3(20)  Model  accelerogram  generated  under  conditions  identical  to 
those  of  Fig.  4.3(18)  except  now  the  rise  time  is  0.5  sec.  The 
ordinate  scale  is  0.08  g/division/m  of  dislocation. 


Fig.  4.3(21)  Model  accelerogram  generated  using  the  fault  geometry  of 
Fig.  4.3(11)  and  the  source-instrument  term  of  Fig.  3.3(4)(source 
two)  with  a rise  time  of  0.1  sec.  Fault  rupture  begins  at  the  hypo- 
center  and  propagates  along  the  fault  at  1.974  km/sec.  The  ordinate 
scale  is  1.73  g/division/m  of  dislocation. 
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receiver  geometry  and  the  medium  seismic  velocities  both  of  which  re- 
main the  same.  Conversely,  the  total  length  of  the  record  of  ground 
motion  appearing  on  these  model  accelerograms  is  much  shorter  since  this 
is  very  definitely  a function  of  fault  rupture  propagation  velocity  as 
well  as  of  fault  geometry  and  seismic  wave  velocities.  The  interval  be- 
tween the  arrival  of  the  first  or  hypocentral  P-wave  at  2.7  sec  and  the 
last  Rayleigh  wave  at  approximately  9.5  sec  is  now  6.8  sec.  This  is  the 
maximum  record  length  we  may  hope  to  match  to  the  Pacoima  Dam  records 
using  the  transfer  function  of  a homogeneous  half-space,  the  geometry 
of  Figure  4.3(11),  and  a fault  rupture  velocity  of  2.538  km/sec.  The 
pulse  duplication  that  we  have  seen  to  occur  under  favorable  circum- 
stances in  Figures  4.3(12)  through  4.3(15)  at  t = r/c  and  t = r/c  + T 
is  also  evident  in  these  source  two  accelerograms,  again  especially  in 
the  surface  wave  coda  near  the  end  of  the  records.  Finally,  the  impor- 
tance of  the  surface  wave  is  again  clearly  distinguishable. 

As  might  have  been  anticipated  from  the  single  source  point  rec- 
ords of  source  one  and  source  two  significant  differences  exist  between 
their  respective  two  hundred  source  point  equivalents.  Host  notably, 
some  constructive  interference  of  the  accelerations  of  individual  sources 
is  evident  in  Figures  4.3(15)  through  4.3(17)  especially  in  association 
with  P-  and  S-wave  arrivals  from  the  hypocentral  region.  The  durations 
of  these  P and  S pulses  are  obviously  related  to  the  rise  time  of  the 
dislocation  source  function  since  in  each  case  the  period  of  the  pulse 
is  equal  to  the  pertinent  rise  time.  Indeed,  a comparison  of  the  ini- 
tial S-wave  motions  of  these  model  accelerograms  with  those  of  the 
Pacoima  Dam  records  of  Figure  4.2(4)  indicates  a marked  similarity  for 
the  longer  model  rise  times.  We  have  for  the  first  time  in  this  study 
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made  a model  acceleration  to  look  something  like  a small  part  of  a real 


I record.  However  there  remain  many  points  of  dissimilarity.  The  over- 

all  frequency  content  of  the  model  records  is  still  too  high  in  that 
individual  record  peaks  are  more  impulsive  than  their  observed  counter- 
parts. We  have  not  matched  the  general  appearance  of  the  actual  Pacoima 
Dam  accelerograms.  Nor  have  any  significant  phases  after  the  initial 
S-wave  been  satisfactorily  duplicated. 

Any  contention  of  minimal  success  for  any  of  the  model  accelero- 
grams appearing  in  Figures  4.3(5)  through  4.3(17)  depends  wholly  on  the 
similarity  between  the  observed  and  generated  horizontal  S-wave  arrival. 
This  similarity  is  best  attained  by  Figure  4.3(17)  where  a rise  time  of 
0.5  sec  is  used.  Indeed,  it  appears  that  a somewhat  better  fit  in  this 
regard  could  be  obtained  by  using  a somewhat  longer  rise  time.  In  view 
of  the  over-all  failure  of  the  model  in  other  respects  we  shall  not 
attempt  the  use  of  a longer  rise  time  but  shall  instead  continue  our 
exploration  of  the  effect  of  the  variation  of  propagation  velocity  of 

I I 

fault  rupture.  Before  leaving  the  records  of  Figures  4.3(15)  through 
4.3(17)  several  other  of  their  features  are  worth  mentioning.  First, 
some  long  period  noise  is  present  in  the  form  of  final  acceleration 

steps.  This  is  again  caused  by  numerical  imperfections  in  the  numerical 

a 

convolution  of  source-instrument  with  geometric  term  that  was  seen  with 
respect  to  Figures  4.3(12)  through  4.3(14).  This  final  non-zero  value 

for  acceleration  is  as  intuitively  unacceptable  now  as  before.  Secondly, 

I 

the  model  peak  accelerations  again  impose  somewhat  prohibitive  limits  on 
allowed  fault  dislocation  amplitude.  Following  previous  convention  the 
- ordinate  scales  of  Figures  4.3(15)  through  4.3(17)  may  be  quoted  as 
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peak  accelerations  of  these  models  to  be  in  accord  with  the  1.25  g ob- 
[ served  peak  acceleration,  fault  dislocation  amplitudes  of  0.14,  0.80, 

and  2.23  m,  respectively,  are  required.  These  are  in  more  reasonable 
agreement  with  observation  but  are  still  rather  small.  Thus,  the  im- 
provements in  model  and  observed  accelerogram  fit  made  by  consideration 
of  source  two  rather  than  source  one  may  be  summaried  as  a better  ini- 
tial S-wave  match  and  more  realistic  peak  accelerations. 

The  opportunity  to  fit  later  segments  of  the  Pacoima  Dam  record 
has  been  degraded  by  employing  the  faster  fault  rupture  propagation 
velocity,  however,  since  the  model  records  are  now  significantly  shorter 
than  their  observed  counterpart.  It  is  certainly  not  clear  whether  this 
is  due  to  inadequacies  in  the  medium  transfer  function  of  Chapter  2 upon 
which  this  study  is  based  or  to  inappropriate  fault  geometry  or  even  to 
errors  in  the  seismic  velocities  employed.  What  is  clear  is  that  longer 
model  accelerograms  will  be  written  by  faulting  schemes  employing  slower 
rupture  propagation  velocities.  Will  these  slower  velocities  also  im- 
prove model  fit  in  ways  other  than  merely  increasing  model  record  length? 
As  we  shall  now  see  the  answer  to  this  question  is  generally  "no". 

In  Figures  4.3(18)  through  4.3(20)  are  presented  accelerograms 
identical  to  those  immediately  preceeding  except  now  a fault  rupture 
propagation  velocity  of  2.221  km/sec  is  used  as  part  of  the  model  fault- 

- ing  process.  To  best  see  the  effect  of  varying  this  velocity  let  us 

I | 

compare  two  model  accelerograms  differing  only  in  this  respect.  In 
particular  Figures  4.3(15)  and  4.3(18)  both  employ  source  two  disloca- 
tion time  histories  with  a rise  time  of  0.1  sec  and  are  otherwise  iden- 
tical  except  the  first  is  associated  with  a fault  rupturing  at  2.538  km/ 

T 

j sec  and  the  second  with  one  rupturing  at  2.221  km/sec.  The  single  most 
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important  point  to  be  made  from  such  a comparison  is  at  once  evident. 
Other  than  with  respect  to  record  length,  these  two  accelerograms  look 
remarkably  similar.  It  is  almost  as  if  we  had  merely  expanded  the  time 
scale  since  the  relative  shapes,  positions,  and  amplitudes  of  all  im- 
portant phase  arrivals  have  remained  unchanged.  This  phenomenon,  which 
will  turn  out  to  be  more  or  less  generally  observed  in  all  succeeding 
similar  model  pairs,  may  be  called  conservation  of  accelerogram  char- 
acter under  a general  transformation  of  uniform  fault  rupture  propaga- 
tion velocity. 

To  the  extent  that  this  observation  is  univerally  applicable  it 
appears  that  once  a given  source-instrument  term  and  rise  time  have  been 
selected  the  remaining  distinctive  character  of  a model  accelerogram  is 
determined  almost  exclusively  by  the  source-receiver  geometry.  Nor  must 
we  look  far  for  the  cause  of  this  phenomenon.  It  is,  perhaps,  inevi- 
table since  so  little  constructive  or  destructive  interference  is  taking 
place  between  the  phase  arrivals  of  individual  point  sources.  Thus,  it 
is  in  the  character  of  the  initial  S-wave  pulse  where  we  have  already 
seen  constructive  interference  to  be  in  operation  that  this  conservation 
principle  is  weakest  and  where  the  variation  of  rupture  propagation 
velocity  has  its  greatest  effect.  This  is  better  seen  in  comparisons 
between  Figures  4.3(16)  and  4.3(19)  and  again  between  Figures  4.3(17) 
and  4.3(20)  where  some  roughness  in  the  initial  S pulse  is  introduced 
at  the  slower  rupture  velocity.  Even  here  the  principle  of  conserva- 
tion of  character  is  not  seriously  violated. 

It  might  be  contended  that  this  phenomenon  makes  the  modeling  of 
accelerograms  subtly  superior  to  the  modeling  of  ground  motion  repre- 
sentations with  lower  frequency  content.  This  argument  is  premature 


here  but  a more  elaborate  discussion  of  this  point  will  be  undertaken 
at  the  conclusion  of  this  study. 

The  ordinate  scales  of  Figures  4.3(18)  through  4.3(20)  may  be 
quoted  as  1.71,  0.32,  and  0.08  g/division/m  of  dislocation.  Under  the 
usual  constraint  of  the  observed  Pacoima  Dam  accelerogram's  peak  accel- 
eration this  implies  fault  dislocations  of  0.18,  0.98,  and  3.91  meters, 
respectively.  These  values  are  in  essential  agreement  with  those  deter- 
mined for  the  previous  three  model  accelerograms  again  demonstrating  the 
power  of  the  conservation  of  character  principle.  The  longer  rise  time 
models  are  obviously  approaching  realistic  fault  dislocation  values. 

The  final  three  model  accelerograms  of  this  series  are  presented 
in  Figures  4.3(21)  through  4.3(23).  These  represent  the  model  records 
generated  under  conditions  identical  to  those  of  the  previous  six  except 
now  a fault  rupture  propagation  velocity  of  1.974  km/sec  is  used.  The 
ordinate  scales  of  these  figures  are  1.73,  0.27,  and  0.07  g/division/m 
of  dislocation.  Again  under  the  constraint  of  the  observed  peak  accel- 
eration this  implies  fault  dislocations  of  0.18,  1.16,  and  4.58  meters, 
respectively.  There  is  little  in  these  records  to  alter  or  add  to  the 
conclusions  reached  from  a study  of  their  six  predecessors. 

Thus  far  our  attempts  to  model  the  Pacoima  Dam  records  with  a 
model  faulting  process  must  be  regarded  as  essential  failures.  The 
effort  has  been  far  from  wasted,  however,  since  very  definite  gains  have 
been  made  in  an  understanding  of  the  effects  of  model  parameter  varia- 
tion upon  the  accelerogram  ultimately  generated.  Therefore,  let  us 
briefly  ignore  the  pertinence  of  Figures  4.3(15)  through  4.3(23)  to  the 
actual  faulting  process  at  San  Fernando  on  February  9,  1971  and  reca- 
pitulate observed  model  accelerogram  changes  expected  in  connection  with 
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variation  of  those  parameters  at  our  disposal  as  evidenced  by  these  last 
nine  figures. 

The  three  major  model  determinants,  in  order  of  increasing  impor- 
tance, are  fault  rupture  propagation  velocity,  source-instrument  term 
time  history,  and  fault-receiver  geometry.  The  first  of  these  has  littl 
effect  on  the  character  of  the  model  accelerogram  but  merely  contracts 
or  expands  the  total  record  along  the  time  axis  in  such  a way  that  the 
frequency  content  remains  unaltered.  Because  a small  amount  of  inter- 
ference among  phase  arrivals  from  individual  point  sources  does  occur 
(most  notably  in  the  source  two  accelerograms)  variation  of  rupture 
propagation  velocity  has  some  effect  on  the  character  of  the  record. 

This  effect  is  very  small,  however,  for  all  cases  so  far  studied.  The 
principal  features  of  a model  accelerogram  determined  by  the  rupture 
velocity  are  total  record  duration  and  separation  of  significant  phases 
for  a given  fa u 1 t- receiver  geometry. 

The  primary  effect  of  the  source-instrument  term  is  to  specify  the 
shape  and  frequency  content  of  each  individual  peak  of  the  model  accel- 
erogram. This  point  has  been  adequately  explored  previously  and  is  well 
demonstrated  in  Figures  4.3(2)  through  4.3(10).  It  will  not  be  belabor- 
ed here.  Again,  due  to  the  essential  lack  of  interference  between 
arrivals  from  individual  point  sources  conclusions  reached  with  regard 
to  single  and  composite  source  records  are  the  same. 

Finally,  the  effect  of  the  source-receiver  geometry  is  to  deter- 
mine the  relative  amplitude  and  position  of  distinctive  phase  arrivals. 
To  a lesser  extent  their  shapes  are  also  affected  especially  where  in- 
terference is  taking  place.  The  geometry  controls  the  most  important 
and  recognizable  features  of  a model  record. 


In  this  content  we  may  expect  to  improve  the  models  thus  far  pro- 
duced by  considering  both  a better  source  term  and  a better  fault  ge- 
ometry relative  to  the  receiving  location.  If  the  summarization  just 
presented  is  valid,  neither  alteration  will  be  sufficient  in  itself  to 
achieve  a good  match  to  the  observed  Pacoima  records  but  more  realistic 
formulation  of  either  parameters  should  cause  a notable  improvement  in 
fit  independent  of  the  other.  Therefore,  let  us  consider  a number  of 
models  generated  using  source-instrument  term  three  in  conjunction  with 
the  familiar  fault  geometry  and  rupture  propagation  velocities.  These 
are  presented  in  Figures  4.3(24)  through  4.3(29). 

In  Figures  4.3(24)  through  4.3(26)  a fault  rupture  propagation 
velocity  of  2.538  km/sec  is  used.  The  source  three  rise  times  are  0.1, 
0.25,  and  0.5  sec,  respectively,  and  the  corresponding  ordinate  scales 
are  2.32,  0.43,  and  0.16  g/division/m  of  dislocation  implying  fault 
dislocations  of  0.13,  0.73,  and  1.95  meters  if  the  peak  record  accelera- 
tion is  to  be  1.25  g as  the  Pacoima  records  imply  it  should  be.  In 
Figures  4.3(27)  through  4.3(29)  are  presented  similar  accelerograms  de- 
rived using  a rupture  velocity  of  1.974  km/sec.  The  three  source  rise 
times  are  again  0.1,  0.25,  and  0.5  sec,  respectively,  and  the  corres- 
ponding ordinate  scales  are  1.63,  0.23,  and  0.06  g/division/m  of  dis- 
location implying  fault  dislocations  of  0.19,  1.36,  and  5.21  meters 
under  the  usual  constraint. 

A comparison  of  these  records  with  their  source  two  equivalents 
indicates  that  in  all  cases  we  have  succeeded  in  generating  somewhat 
smoother  versions  of  these  model  accelerograms.  Several  longer  period 
features  (in  the  range  of  two  to  four  Hertz)  have  become  accentuated 
and  most  of  the  shorter  period  features  (in  the  range  of  ten  to  twenty 


137 


Fig.  4.3(24)  Model  accelerogram  generated  using  the  fault  geometry  of 
Fig.  4.3(11)  and  the  source-instrument  term  of  Fig.  3.3(6)(source 
three)  with  a rise  time  of  0.1  sec.  Fault  rupture  begins  at  the 
hypocenter  and  propagates  along  the  fault  at  2.538  km/sec.  The 
ordinate  scale  is  2.32  g/division/m  of  dislocation. 


4.3(25)  Model  accelerogram  generated  under  conditions  identical  to 
those  of  Fig.  4.3(24)  except  now  the  rise  time  is  0.25  sec.  The 
ordinate  scale  is  0.*3  g/division/m  of  dislocation. 


Fig.  4.3(26)  Model  accelerogram  generated  under  conditions  identic 
those  of  Fig.  4.3(24)  except  now  the  rise  time  is  0.5  sec.  Th 
ordinate  scale  is  0.16  g/division/m  of  dislocation. 


141 


- *Y 


Fig.  4.3(28)  Model  accelerogram  generated  under  conditions  identical  to 
those  of  Fig.  4.3(27)  except  now  the  rise  time  is  0.25  sec.  The 
ordinate  scale  is  0.23  g/division/m  of  dislocation. 
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Hertz)  have  been  de-emphasized.  The  final  acceleration  values  which 
have  proved  annoying  before  have  now  been  reduced  to  essentially  zero. 
Otherwise,  remarkable  similarity  exists  between  all  source  two-source 
three  model  pairs. 

A comparison  of  the  accelerograms  of  Figures  4.3(24)  through 
4.3(29)  with  the  resolved  Pacoima  Dam  records  of  Figure  4.2(4)  indi- 
cates that  these  last  models  are  apparently  capable  of  matching  the 
frequency  character  of  the  observed  accelerograms  at  least  in  gross 
detail.  Thus,  one  of  the  features  of  previous  models  that  was  found 
objectionable  has  been  seemingly  overcome.  As  was  the  case  with  all 
previous  accelerograms,  however,  a true  match  of  significant  arrivals 
is  lacking.  Again,  if  the  summary  of  the  effect  of  parameter  variation 
given  above  is  valid,  as  these  last  six  records  would  seem  to  substan- 
tiate, this  situation  will  be  improved  by  a consideration  of  a more 
accurate  fault  geometry.  This  is  attempted  in  the  next  section. 

4.32  Model  Accelerograms  from  a Hinged  Fault 

As  discussed  in  Chapter  1 the  planar  fault  model  employed  by 
Canitez  and  Toksoz  (1972),  Mikumo  (1973),  and  Trifunac  (1974)  to  explain 
various  static  and  dynamic  consequencies  of  the  faulting  process  associ- 
ated with  the  San  Fernando  earthquake  of  February  9,  1971  is  unsatis- 
factory in  several  respects.  As  has  been  shown  by  Turnbull  and  Battis 
(1973)  and  by  Hanks  (1974)  some  of  these  deficiencies  may  be  qualita- 
tively and  quantitatively  reduced  by  considering  a fault  surface  that 
is  hinged  in  some  manner  so  that,  generally,  the  dip  of  the  fault  plane 
increases  with  depth.  We  now  wish  to  see  whether  model  strong-motion 
records  may  also  be  improved  by  the  use  of  such  a fault  model. 

Any  proposed  fault  must  conform  to  certain  constraints.  Thus,  the 
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hypocentral  position,  or  point  of  initial  rupture,  must  lie  on  the 
fault  surface,  the  dip  of  the  fault  surface  must  be  near  52°  at  the 
hypocenter  and  must  conform  to  observed  surface  faulting,  and  the  ge- 
ometry of  the  surface  fault  rupture  relative  to  the  epicenter  (and  to 
Pacoima  Dam)  must  be  maintained.  A model  incorporating  all  of  these 
constraints  is  shown  in  Figure  4.3(30).  This  fault  model  ir  a somewhat 
simplified  version  of  the  one  derived  by  Alewine  and  Jordan  (1973)  using 
the  closed  form  analytic  expressions  for  displacement  in  an  elastic 
half-space  due  to  an  inclined  rectangular  fault  given  by  Mansinha  and 
Smylie  (1971).  These  displacements  are  calculated  for  each  of  forty- 
seven  appropriately  oriented  fault  elements  constrained  to  conform  to 
other  observed  geophysical  properties  of  the  fault  as  outlined  above 
and  then  summed  to  give  the  final  displacements  for  the  model.  Not  con- 
sidered here,  but  included  in  the  model  of  Alewine  and  Jordan  are  a 
system  of  splay  faults  near  the  surface  and  a component  of  normal  fault 
motion  on  the  San  Gabriel  fault.  The  first  feature  is  an  unnecessary 
complication  at  this  level  of  the  study  since  it  can  only  affect  the 
very  end  of  any  model  accelerogram  generated  in  the  manner  of  those  we 
have  previously  considered  where  rupture  starts  at  the  hypocenter  and 
propagates  to  the  surface.  The  second  feature  is  omitted  since  it  is 
not  clear  whether  any  motion  on  the  San  Gabriel  fault  is  better  associ- 
ated with  regional  stress  adjustments  after  the  main  earthquake  or  with 
the  main  earthquake  itself.  The  part  of  the  model  used  in  this  study 
is  in  general  agreement  with  that  employed  by  Hanks  (1974). 

We  now  proceed  much  as  before.  That  is,  to  model  the  ground  motion 
at  Pacoima  Dam  due  to  motion  along  the  entire  fault,  the  ground  motions 
from  two  hundred  individual  sources  evenly  distributed  over  the  fault 
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Fig.  4.3(30)  Hinged  model  of  the  fault  associated  with 
the  San  Fernando  earthquake  of  February  9,  1971 
(after  Alewine  and  Jordan,  1973).  The  hypocenter 
of  the  earthquake  is  at  the  lower  end  of  the  fault. 
The  sense  of  fault  motion,  the  relative  position  of 
Pacoima  Dam,  and  the  seismic  velocities  used  in 
subsequent  calculations  are  also  shown. 
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surface  are  calculated  and  stored  on  magnetic  tape.  The  summation  of 
these  ground  motions  is  then  done  in  such  a way  as  to  simulate  propa- 
gation of  the  dislocation  episode  along  the  fault  at  some  fault  rupture 
velocity.  Only  discrete  values  of  this  velocity  are  possible  if  addi- 
tional interpolation  between  points  is  to  be  avoided  and  since  the  total 
fault  length  is  somewhat  greater  in  this  case  while  the  number  of  point 
sources  has  remained  the  same  these  discrete  velocity  values  will  now 
be  somewhat  different.  From  Figure  4.3(30)  we  see  that  the  total  fault 
length  is  now  21.66  km.  The  allowable  rupture  propagation  velocities 
within  the  range  of  interest  are  now  1.969,  2.166,  and  2.407  km/sec 
corresponding  to  time  steps  between  points  of  0.055,  0.050,  and  0.045 
sec,  respectively.  It  remains  but  to  convolve  this  summed  ground  motion 
with  the  various  source- instrument  terms  of  Chapter  3 in  the  familiar 
manner. 

Using  this  new  fault  model  let  us  begin  by  considering  the  sim- 
plest source.  Model  accelerograms  employing  the  source-instrument  term 
one,  the  model  geometry  of  Figure  4.3(30),  and  a fault  rupture  propaga- 
tion velocity  of  1.969  km/sec  are  presented  in  Figures  4.3(31)  through 
4.3(33)  for  rise  times  of  0.1,  0.25,  and  0.5  sec,  respectively.  The 
ordinate  scales  are  15.36,  6.11,  and  2.91  g/division/m  of  dislocation 
implying  fault  displacements  of  0.02,  0.05,  and  0.11  meters  under  the 
constraint  imposed  by  the  observed  Pacoima  peak  acceleration. 

These  model  accelerograms  are  disappointing  in  that  they  do  not 
match  the  observed  records  of  Figure  4.2(4)  much,  if  any,  better  than 
their  planar  fault  counterparts . Their  frequency  content  is  too  high 
and  the  constraints  imposed  on  the  amplitude  of  fault  dislocation  are 
unacceptable.  Also,  with  the  re-introduction  of  source-instrument  term 
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one,  numerical  noise  is  again  present  in  the  form  of  a step  offset  in 
acceleration  remaining  after  all  elastic  energy  has  arrived  from  the 
source.  Differences  do  exist,  however,  arising  from  the  different  ge- 
ometry. Several  could  easily  have  been  anticipated  such  as  the  length- 
ening of  the  S-P  wave  interval  from  the  point  of  initial  rupture  re- 
sulting from  the  increased  hypocenter-station  distance  and  the  increased 
total  record  length  due  to  the  longer  fault  surface  of  Figure  4.3(30). 
Perhaps  also  to  be  expected  is  the  increase  in  peak  acceleration  since 
the  closest  approach  of  fault  to  receiver  is  reduced  for  the  geometry 
here  employed.  Other  changes  might  not  have  been  so  easily  anticipated, 
however,  such  as  the  appearance  of  an  important  phase  arrival  near  the 
middle  of  the  vertical  component  time  history  which  was  not  previously 
apparent.  The  distinct  features  of  this  geometry  will  be  considered  in 
greater  detail  later  when  a more  satisfactory  source-instrument  term 
will  be  used. 

Although  Figures  4.3(31)  through  4.3(33)  have  not  obviously  ad- 
vanced the  search  for  a fruitful  model  of  the  faulting  process  associ- 
ated with  the  San  Fernando  earthquake  in  any  significant  manner,  they 
have  admirably  fulfilled  their  principal  aim  which  was  to  demonstrate 
the  invariance  of  what  has  been  termed  the  conservation  of  character 
under  a transformation  of  model  geometry.  These  three  figures  are 
similar  to  their  planar  predecessors  exactly  as  we  would  expect  and  are 
also  dissimilar  in  precisely  the  expected  manner. 

Let  us  now  go  on  to  consider  one  final  series  of  model  accelero- 
grams. That  is,  we  wish  to  generate  those  three  sets  of  records  with 
fault  rupture  propagation  velocities  of  2.407,  2.166,  and  1.969  km/sec 
using  the  fault  geometry  after  Alewine  and  Jordan  (1973)  as  shown  in 


151 


Figure  4.3(30)  and  source-instrument  term  three.  Since  each  source- 
instrument  term  is  evaluated  for  three  rise  times  a complete  series  of 
this  type  consists  of  nine  model  accelerograms.  Actually,  the  record 
of  this  set  employing  a rupture  propagation  velocity  of  2.407  km/sec  and 
a rise  time  of  0.1  sec  was  inadvertently  omitted.  The  remaining  eight 
accelerograms  are  shown  in  Figures  4.3(34)  through  4.3(41).  Since  the 
procedure  by  which  these  model  records  are  presented  is  somewhat  un- 
avoidably repetitious,  we  shall  list  their  fundamental  properties  here 
in  a very  straightforward  way  and  defer  more  elaborate  considerations 
until  all  eight  records  have  been  treated. 

In  Figures  4.3(34)  and  4.3(35)  are  shown  the  two  models  of  this 
series  evaluated  for  a rupture  velocity  of  2.407  km/sec.  The  source 
function  rise  times  are  0.25  and  0.5  sec,  respectively,  and  the  ordi- 
nate scales,  again  respectively,  are  0.43  and  0.14  g/division/m  of  dis- 
location implying  total  fault  dislocation  amplitudes  of  0.72  and  2.19 
meters  if  the  maximum  model  accelerations  are  to  agree  with  the  1.25  g 
value  observed  at  Pacoima  Dam. 

Model  accelerograms  employing  the  source-instrument  term  three, 
the  model  geometry  of  Figure  4.3(30),  and  a fault  rupture  propagation 
velocity  of  2.166  km/sec  are  presented  in  Figures  4.3(36)  through 
4.3(38)  for  rise  times  of  0.1,  0.25,  and  0.5  sec,  respectively.  The 
ordinate  scales  are  4.18,  0.27,  and  0.09  g/division/m  of  dislocation 
implying  fault  displacements  of  0.07,  1.14,  and  3.64  meters  under  the 
constraint  imposed  by  the  observed  Pacoima  Dam  peak  acceleration. 

Finally,  model  accelerograms  of  this  series  employing  a fault 
rupture  propagation  velocity  of  1.969  km/sec  are  presented  in  Figures 
4.3(39)  through  4.3(41)  for  rise  times  of  0.1,  0.25,  and  0.5  sec. 
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Fig.  4.3(34)  Model  accelerogram  generated  using  the  fault  geometry  of 
Fig.  4.3(30)  and  the  source-instrument  term  of  Fig.  3.3(6)  (source 
three)  with  a rise  time  of  0.25  sec.  Fault  rupture  begins  at  the 
hypocenter  and  propagates  along  the  fault  at  2.407  km/sec.  The 
ordinate  scale  is  0.43  g/division/m  of  dislocation. 
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Fig.  4.3(36)  Model  accelerogram  generated  using  the  fault  geometry  of 
Fig.  4.3(30)  and  the  source-instrument  term  of  Fig.  3.3(6)  (source 
three)  with  a rise  time  of  0.1  sec.  Fault  rupture  begins  at  the 
hypocenter  and  propagates  along  the  fault  at  2 . 1 66  km/sec.  The 
ordinate  scale  is  4.18  g/division/m  of  dislocation. 


Fig.  4.3(38)  Model  accelerogram  generated  under  conditions  identical  to 
those  of  Fig.  4.3(36)  except  now  the  rise  time  is  0.5  sec.  The 
ordinate  scale  is  0.09  g/division/m  of  dislocation. 
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respectively.  The  ordinate  scales  are  4.05,  0.22,  and  0.06  g/division/ 
m of  dislocation  implying  fault  displacements  of  0.08,  1.40,  and  5.58 
meters  under  the  usual  constrain  imposed  by  the  observed  Pacoima  Dam 
peak  acceleration. 

The  first  and  most  significant  contention  that  may  be  made  about 
the  records  of  this  series  is  that  they  look  rather  like  the  resolved 
Pacoima  Dam  records  we  wish  to  match.  Since  the  fundamental  difference 
between  the  input  parameters  of  these  models  and  all  previous  efforts  is 
in  the  fault  geometry  used,  it  is  obvious  that  the  current  fault  ge- 
ometry is  more  in  accord  with  the  actual  faulting  that  took  place  during 
the  San  Fernando  earthquake.  This  is  hardly  surprising  in  view  of  the 
wealth  of  independent  evidence  of  previous  studies.  What,  perhaps,  is 
surprising  is  the  support  by  high-frequency  model  motions  of  inferences 
derived  from  much  longer  period  data.  This  is  a result  of  far-reaching 
import. 

Let  us  begin  by  compiling  a list  of  features  of  the  model  records 
which  are  believed  to  be  in  substantial  agreement  with  the  resolved 
Pacoima  Dam  records  of  Figure  4.2(4).  r’ost  striking  is  the  manner  in 
which  the  relative  amplitudes  of  the  model  accelerograms  match  their 
observed  counterparts.  The  term  relative  amplitudes  is  used  here  to 
describe  both  the  variation  of  amplitude  with  time  of  either  component 
and  the  ratio  at  any  given  time  of  the  amplitude  of  one  component  to 
that  of  the  other.  Specifically  the  record  amplitudes  between  the 
arrival  of  the  first  P-wave  and  the  first  S-wave  are  quite  small,  the 
horizontal  component  motion  being  somewhat  smaller  than  the  vertical  in 
agreement  with  observation.  Analysis  of  all  the  eight  records  of  this 
series  shows  that  this  conclusion  is  generally  valid  although  details 
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differ  somewhat  from  record  to  record.  The  arrival  of  the  S-wave  energy 
precipitates  a 1 to  1.5  second  packet  of  much  greater  amplitude  motion 
which  is  somewhat  more  prolonged  and  energetic  on  the  horizontal  com- 
ponent. Again,  the  feature  of  the  model  records  which  is  quite  general 
is  in  general  agreement  with  observation.  Following  the  arrival  of  the 
shear  energy  from  the  hypocentral  region,  an  interval  of  several  seconds 
during  which  a marked  decrease  in  record  amplitude  occurs  is  seen.  This 
decrease  is  more  noticeable  on  the  vertical  record.  After  this  interval 
another  energetic  arrival  is  found  on  the  vertical  component  while  the 
horizontal  component  remains  quiet.  At  this  point  the  match  between 
model  and  observed  records  becomes  much  less  satisfactory.  However,  as 
we  shall  shortly  see  the  model  time  histories  between  the  first  P-wave 
arrival  and  the  first  prominent  arrival  after  the  hypocentral  region 
S-wave  are  generated  by  motions  along  the  fault  from  the  hypocenter  to 
a point  less  than  two  kilometers  below  the  ground  surface.  Thus,  con- 
siderations of  relative  amplitudes  alone  would  indicate  that  over  the 
initial  80  percent  of  its  length  the  model  fault  of  Figure  4.3(30)  is 
adequate  to  explain  the  observed  Pacoima  Dam  records.  This,  in  turn, 
leads  to  several  important  conclusions.  Since  the  amplitude  configura- 
tion described  above  is  generally  true  for  a variety  of  fault  rupture 
propagation  velocities  and  source  rise  times,  it  must  be  approximately 
independent  of  both,  at  least  over  the  ranges  of  these  parameters  that 
have  been  investigated.  This  is,  of  course,  completely  equivalent  to 
stating  that  the  relative  amplitudes  of  model  accelerograms  are  power- 
fully dependent  on  the  source-receiver  geometry.  The  same  is  presum- 
ably true  of  real  accelerograms. 

Since  all  the  model  records  of  this  series  were  generated  assuming 
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a constant  fault  rupture  propagation  velocity  and  a constant  dislocation 
amplitude  over  the  entire  fault  surface,  it  is  clear  that  rapid  or  large 
variations  of  these  parameters  are  not  necessary  to  explain  the  first 
six-to-seven  seconds  of  the  Pacoima  record.  This  is  certainly  a pleas- 
ing result  from  a modeling  standpoint.  More  importantly,  it  specifi- 
cally argues  against  the  conceptul ization  of  strong  high-frequency  mo- 
tions as  the  result  of  an  essentially  random  process,  at  least  in  the 
particular  case  of  the  San  Fernando  earthquake. 

At  this  point  we  must  consider  a problem  that  has  been  previously 
ignored  primarily  because  model  accelerograms  generated  had  not  achiev- 
ed a fit  to  the  observed  data  of  a degree  to  warrant  its  consideration. 
As  mentioned  above,  even  the  best  models  thus  far  generated  are  success- 
ful  only  in  matching  the  first  6.5  seconds  of  observed  record  motion. 
However,  the  actual  peak  acceleration  of  1.25  g,  which  has  been  used  to 
calculate  the  implications  of  any  given  model  to  average  fault  disloca- 
tion ampltiude,  occurs  after  this  time.  The  peak  acceleration  in  that 
part  of  the  actual  Pacoima  Dam  records  for  which  a good  model  fit  is 
obtained  is  much  less  than  this  value  being  approximately  0.6  g.  It  is 
this  value  which  should,  more  properly,  be  used  to  calculate  the  aver- 
age fault  dislocation  amplitudes.  This  observation  only  strengthens 
the  conclusion  of  the  general  unacceptability  of  early  model  motions. 

The  implication  of  this  refinement  to  the  average  fault  dislocation 
amplitudes  to  be  inferred  from  any  model  is  to  approximately  halve 
previously  quoted  values. 

Are  there  other  features  of  the  observed  records  that  may  be 
duplicated  by  some  model?  The  answer  is  apparently  yes.  For  several 
model  accelerograms,  most  notably  for  that  employing  a fault  rupture 
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propagation  velocity  of  1.969  km/sec  and  a source  rise  time  of  0.25 
sec,  the  character  of  the  vertical  component  initial  shear  wave  arrival 
on  the  Pacoima  Dam  records  is  rather  well  duplicated.  This  motion  is 
a distinctive  series  of  approximately  4 Hertz  oscillations  of  approxi- 
mately one  second  duration.  The  next  really  distinctive  arrival  on  the 
vertical  component  is  a single  pulse  of  about  the  same  frequency  occur- 
ring some  three  seconds  after  the  onset  of  the  S-wave  motion.  This 
pulse  is  faithfully  duplicated  on  virtually  all  model  accelerograms  but 
has  the  Droper  frequency  content  only  on  those  with  0.25  sec  rise  times. 

There  are  also  notable  failures  of  any  of  the  artificial  records 
to  match  some  features  of  the  observed  motion.  The  comparison  between 
model  and  real  horizontal  motion  is  generally  much  less  satisfactory 
than  a similar  comparison  between  vertical  motions.  This  is  difficult 
to  explain  in  terms  of  minor  model  parameter  variation  especially  for 
the  first  horizontal  S-wave  pulse  which  apparently  have  significantly 
different  frequency  content  on  the  vertical  and  horizontal  Pacoima  Dam 
records.  For  the  moment  let  it  merely  be  contended  that  the  best  model 
found  of  those  tried  is  that  employing  the  Alewine  and  Jordan  derived 
fault  geometry,  a fault  rupture  propagation  velocity  of  approximately 
2 km/sec,  and  a rise  time  of  0.25  sec.  This  model  rather  well  dupli- 
cates the  first  6.5  seconds  of  the  recorded  Pacoima  Dam  vertical  com- 
ponent both  with  respect  to  amplitude  and  phase  character  but  does  a 
considerably  less  impressive  job  of  matching  the  recorded  horizontal 
component  with  respect  to  the  character  of  the  individual  phases  al- 
though the  overall  amplitude  variation  with  time  relative  to  the  verti- 
cal component  is  fairly  successful. 
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4.4  Discussion  of  the  Results 

The  purpose  of  this  section  is  essentially  to  set  reasonable  lim- 
its of  uncertainty  on  the  basic  model  parameters.  The  achievement  of 
this  goal  is  aided  by  the  deterministic  nature  of  the  modeling  process. 
That  is,  we  have  access  to  the  origin  of  each  feature  of  the  model  mo- 
tion either  explicitly  through  consideration  of  the  equations  governing 
it  or  implicitly  by  observation  of  the  effect  of  parameter  variation. 

It  is  possible  to  pick  a distinctive  phase  and  unambiguously  state  that 
it  is  due  to  a given  type  of  elastic  wave  arrival  from  a specific  part 
of  the  model  fault.  Insofar  as  we  may  confidently  pair  such  an  arrival 
with  one  on  the  observed  record  the  same  conclusions  hold  for  the  actual 
record  motion. 

Let  us  begin  in  the  most  straightforward  way  imaginable.  That  is, 
let  us  calculate  for  some  particular  model  the  arrival  times  of  body 
wave  phases  from  significant  parts  of  the  fault.  In  this  analysis, 
which  is  wholly  an  exercise  in  geometry,  we  will  use  the  model  accelero- 
gram of  Figure  4.3(40)  which  was  found  to  be  a fairly  successful  attempt 
to  reproduce  the  observed  accelerations.  As  we  have  seen  this  model 
employs  a fault  derived  from  that  of  Alewine  and  Jordan  (1973),  a rup- 
ture propagation  velocity  of  1.969  km/sec  and  a rise  time  of  0.25  sec 
for  source-instrument  term  three.  The  points  of  interest  along  this 
fault  from  which  phase  arrivals  will  be  calculated  are,  in  order  of 
origin  along  the  fault  with  decreasing  depth: 

i)  the  hypocenter  - point  a 

ii)  the  point  at  which  the  slope  of  the  fault  changes 
from  52°  to  40°  - point  b 

iii)  the  point  at  which  the  slope  of  the  fault  changes 
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from  40°  to  35°  - point  c 

iv)  the  lower  point  of  intersection  of  the  fault  plane 

with  the  wedge  of  0 £ 0cr  where  6cr  is  the  critical 

angle  as  discussed  in  Chapter  2 and  where  0 and  0cf. 

are  measured  from  the  vertical  at  the  receiving 

station  - point  0L 
cr 

v)  the  upper  point  of  intersection  of  the  fault  plane 
with  the  wedge  of  0 £ 0cr  - point  0^r 

vi)  the  point  of  closest  approach  of  the  fault  plane  to 
the  receiving  station  - point  CA 

vii)  the  point  at  which  the  slope  of  the  fault  plane 
changes  from  35°  to  30°  - poind  d 

viii)  the  point  of  intersection  of  the  fault  plane  with 
the  surface  of  the  half  space  - point  e 
The  coordinates  of  these  points  as  well  as  of  the  receiving  sta- 
tion are  given  in  Table  4-1.  The  reasons  for  selecting  these  particular 
points  is  quite  obvious.  In  the  case  of  the  fault  hinge  points  we  ex- 
pect a reordering  of  the  importances  of  the  various  Green's  function 
terms  to  occur  as  may  be  seen  from  a consideration  of  equation  3.4-5. 
Between  and  0^  no  energy  in  the  form  of  SP  motion  may  arrive  at 
the  receiving  station.  Thus  it  is  of  interest  to  look  for  differences 
in  the  character  of  the  record  in  this  interval  that  might  distinguish 
it  from  the  remainder  of  the  record.  Finally,  the  point  of  closest 
approach  is  a superficially  obvious  place  from  which  to  expect  large- 
ampl i tude  motion. 

As  in  Table  4-1  let  us  designate  the  model  equivalent  of  the  re- 
ceiving station  at  Pacoima  Dam  as  point  PD.  Then,  where  ab  is  the 
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length  of  the  line  segment  between  points  a and  b and  where  v^,  v$,  and 
are  the  compressional  wave,  shear  wave,  and  fault  rupture  propagation 
velocities,  respectively,  the  arrival  of  the  P-wave  from  point  c is 
simply 

ab/Vf  + bc/v.p  + cPD/Vp 

The  time  interval  between  the  arrival  of  this  phase  and  the  arrival  of 
the  hypocentral  P-wave  is  the  difference  between  the  above  quantity  and 
aPD/Vp.  In  this  way  a number  of  relative  or  absolute  arrival  times  may 
be  worked  out.  The  computation  of  the  arrival  times  of  SP  phases  are 
somewhat  more  complex  but  hardly  more  difficult. 

In  Figure  4.4(1)  are  presented  the  arrival  times  of  the  compres- 
sional phases  from  each  of  the  above  eight  points  super  -sed  on  the 
model  accelerogram  of  Figure  4.3(40).  Only  one  arrival  is  shown  for  P- 
waves  from  points  c and  6 . This  is  due  to  the  extreme  proximity  of 

these  two  points  for  the  particular  fault  geometry  used  as  may  be  seen 
by  comparing  their  coordinates  in  Table  4-1.  In  Figures  4.4(2)  and 
4.4(3)  are  presented  similar  plots  for  S-  and  SP-wave  arrivals.  Again, 
in  both  cases  arrival  times  from  points  c and  6*"  are  essentially 
indistinguishable. 

Perhaps  the  single  most  impressive  feature  of  these  three  plots 
is  the  complete  unconcern  of  the  model  motions  for  my  seismol ogical  pre- 
conceptions. With  a few  exceptions  to  be  discussed  below,  not  a single 
important  feature  of  the  model  records  can  be  attributed  to  phase  arri- 
vals from  any  of  the  points  proposed  as  significant. 

The  most  notable  exception,  of  course,  is  the  initial  S-P  interval. 
This  interval,  which  is  not  dependent  upon  the  fault  rupture  propagation 
velocity  but  only  upon  the  hypocenter-receiver  separation  and  the  shear- 


4.4(1)  P-wave  arrivals  from  the  significant  points  of  Table  4-1  for  a fault 
rupture  propagation  velocity  of  1.969  km/sec.  The  model  accelerogram  upon 
which  these  arrival  times  are  superposed  is  that  of  Fig.  4.3(40). 
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and  compressional -wave  velocities,  is  2.91  seconds.  That  is  a somewhat 
higher  value  than  observation  would  indicate.  Assuming  that  the  Pacoima 
Dam  accelerometers  were  tripped  very  near  the  initial  P-wave  arrival  and 
that  the  S-wave  arrival  is  signaled  in  actuality  by  the  onset  of  long- 
period  motion  on  the  horizontal  records,  a more  realistic  S-P  interval 
for  the  observed  case  is  approximately  2.5  seconds. 

The  model  and  observed  intervals  may  be  reconciled  by  decreasing 
the  model  hypocenter-recei ver  distance,  by  decreasing  the  model  compres- 
sonal  wave  velocity,  or  by  increasing  the  model  shear-wave  velocity. 

The  last  alternative  seems  the  most  satisfactory  since  the  S-wave 
velocity  employed  is  lower  than  that  of  other  studies--2.9  as  opposed 
to  3.2  or  3.3  km/sec.  A shear-wave  velocity  of  3.1  km/sec  for  this  ge- 
ometry and  a P-wave  velocity  of  5.58  km/sec  produces  an  S-P  interval  of 
2.52  sec  in  much  better  agreement  with  observation.  We  shall  keep  this 
point  in  mind  as  a source  of  uncertainty  in  other  faulting  parameters 
which  are  implicitly  dependent  on  the  S-wave  velocity. 

Another  minor  exception  to  the  failure  of  the  generated  motion  to 
obey  our  preconceptions  is  in  the  case  of  the  P-wave  arrival  from  point 
b.  A small  but  definite  increasein  P-wave  amplitude  is  seen  to  occur 
at  the  arrival  of  compressional  energy  from  this  point.  Can  we  explain 
the  general  unimportance  to  model  motion  of  the  rapid  changes  in  fault 
dip?  It  turns  out  that  we  can  and  that  the  explanation  is  surprisingly 
simple. 

Consider  from  equation  3.4-5  the  fault  slope  deoendence  of  the 
total  solution 

u^ (horizontal ) = p{sin20[G^  ^ - G-^  3] 
c°s20[G11>3  + G13j]} 
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^(vertical)  = p{sin20[G31  ^ - G33  3) 

- cos20[G33>1  + G31 >3]} 

In  these  equations  0 is  the  dip  of  the  fault.  Thus  sin26  and  cos26  may 
be  thought  of  as  the  coefficients  due  to  the  fault  orientation  which 
determine  which  of  the  spatial  derivatives  of  the  Green's  functions-- 
G^  p — wi 1 1 be  of  dominant  importance.  Unlike  the  source-receiver  angu- 
lar dependence  buried  within  the  Green's  functions  (also,  somewhat 
unfortunately,  symbolized  by  0 in  Chapter  2)  which  is  varying  constantly 
as  the  fault  rupture  process  propagates  from  the  hypocenter  to  the  half- 
space surface  these  26  dependent  coefficients  remain  constant  for  con- 
stant fault  plane  orientation.  Thus,  in  going  from  the  fault  segment 
between  points  a and  b to  that  between  points  b and  c,  20  changes  from 
104°  to  80°  and  the  factors  sin20  and  cos20  change  from  0.970  and  -0.242 
to  0.985  and  0.174,  respectively.  This  is  simply  not  a large  effective 
difference.  Apparently  it  is  virtually  insignificant  when  compared  to 
the  variation  of  the  Green's  functions  themselves  over  rather  small  dis- 
tances. Similarly,  the  variations  of  sin20  and  cos26  between  the  fault 
segments  be  and  cd  are  0.985  to  0.940  and  0.174  to  0.342.  Finally,  be- 
tween segments  cd  and  de  the  variation  in  sin26  is  from  0.940  to  0.866 
while  that  of  cos20  is  from  0.342  to  0.500.  These  changes,  if  they  are 
visible  at  all  are  obscured  by  the  complexity  of  motions  in  the  latter 
part  of  the  record  where,  unlike  the  case  of  compressional  arrivals  from 
point  b,  different  elastic  wave  types  are  simultaneously  arriving  from 
different  parts  of  the  fault. 

Figure  4.4(3)  graphically  shows  that  very  little  of  the  character 
of  the  model  accelerogram  can  be  attributed  to  the  SP  phase  as  defined 
in  Chapter  2.  This  may  very  well  be  due  to  the  small  amplitude  of  this 
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phase  relative  to  the  other  body-wave  phases,  especially  the  S-wave. 

In  particular,  that  section  of  the  model  record  during  which  no  SP  mo- 
tion is  allowed  (between  the  SP  arrivals  from  e^r  and  6^r  in  Figure 
4.4(3))  is  not  apparently  different  from  those  portions  directly 
adjacent. 

This  leaves  only  the  arrivals  from  the  point  of  closest  approach 
of  fault  to  receiver.  Any  contention  that  this  point  is  significant 
must  apparently  be  regarded  as  unconvincing.  Simply  very  little  is 
happening  at  the  appropriate  time.  This  must  be  due  to  the  fact  that 
a point  dislocation  does  not  radiate  energy  independently  of  azimuth 
but  prefers  some  directions  to  others.  Thus  the  amplitude  at  the  re- 
ceiver is  greatly  affected  by  the  source-receiver  angular  relationship 
and  not  merely  by  source-recei ver  separation.  The  necessary  angular 
condition  for  large-amplitude  motion  is  lacking  in  this  case.  Con- 
versely, they  must  be  fulfilled  in  the  case  of  the  aforementioned 
single  pulse  visible  on  all  model  accelerograms  of  the  last  series. 

This  pulse  represents  one  of  the  most  convincing  equivalencies  between 
model  and  observed  motion  so  it  is  worthwhile  to  dwell  on  its  implica- 
tions for  a moment. 

Assuming  that  this  pulse  is  indeed  equivalent  to  that  on  the 
Pacoima  Dam  vertical  component  at  approximately  5.8  seconds  after  the 
onset  of  the  record  motion,  what  may  we  infer  about  the  actual  fault 
rupture  history?  It  seems  most  likely  that  the  model  pulse  is  due  to 
an  S-wave  arrival  from  some  point  on  the  fault  near  its  upper  end.  This 
is  supported  by  Figure  4.3(1)  which  shows  the  ground  displacement  from 

, a single  point  dislocation  with  the  Heaviside  unit  step  source  function. 

; In  this  figure,  in  which  the  source-receiver  geometry  is  fairly 
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comparable  to  one  appropriate  for  a source  near  the  upper  end  of  the 
fault  of  Figure  4.3(30)  the  amplitude  of  the  S-wave  motion  is  seen  to 
be  completely  predominant.  Given,  then,  that  this  pulse  on  the  model 
accelerograms  is  due  to  some  shear-wave  arrival,  from  what  point  on  the 
fault  does  it  come? 

From  Figure  4.4(2)  we  see  that  it  must  be  from  a point  between 
that  of  closest  approach  and  point  d.  Indeed,  once  we  have  selected  an 
appropriate  arrival  time  for  this  model  pulse  it  is  trivial  to  construct 
equations  of  constraint  and  solve  for  the  source  point  uniquely.  The 
exact  time  of  pulse  onset  is  somewhat  problematical  but  7.2  seconds  after 
the  first  P-wave  arrival  seems  a good  estimate.  The  S-wave  arriving  at 
this  time,  given  a fault  rupture  propagation  velocity  of  1.969  km/sec 
and  model  compressional  and  shear  velocities  of  5.58  and  2.90  km/sec 
can  only  come  from  the  point  on  the  fault  whose  coordinates  are  (13.32, 
2.08).  The  important  thing  to  note  is  that  alteration  of  any  of  these 
three  velocities  will  change  the  arrival  time  of  the  model  pulse  but 
not  its  character.  At  least  this  is  what  all  our  previous  experience 
with  parameter  variation  would  indicate.  We  may  bring  the  model  inter- 
val of  7.2  seconds  into  closer  agreement  with  the  observed  interval  of 
5.8  seconds  either  by  increasing  the  rupture  propagation  velocity  or 
the  S-wave  velocity  or  by  decreasing  the  P-wave  velocity.  As  before  in 
discussing  the  S-P  interval,  the  last  alternative  seems  least  satis- 
factory. Changing  the  shear  velocity  alone  is  not  likely  to  be  ade- 
quate since  the  source-receiver  separation  is  rather  small  (3.43  km). 

This  implies  that  a large  change  in  the  S-wave  velocity  would  be  neces- 
sary. In  fact  no  increase  in  shear  velocity  alone  will  suffice  since 
at  the  time  the  pulse  in  question  appears  on  the  observed  record  the 


.3*-'  , ‘ 


175 

model  has  not  yet  propagated  to  the  critical  point  on  the  fault.  There- 
fore, it  must  be  that  the  difference  between  model  and  observed  inter- 
vals is  primarily  due  to  the  difference  between  model  and  actual  fault 
rupture  propagation  velocity. 

If  the  entire  difference  is  to  be  so  explained  the  implied  actual 
rupture  velocity  is  2.32  km/sec.  If  the  shear  velocity  is  increased 
somewhat  to  3.1  km/sec  as  we  have  seen  is  'indicated  by  the  observed  S-P 
interval  the  pertinent  fault  rupture  propagation  velocity  falls  to  2.30 
km/sec.  Thus,  it  would  seem  that  if  we  have  properly  matched  model  and 
observed  phases  the  average  rupture  velocity  of  the  faulting  process 
over  a very  substantial  part  of  the  fault  may  be  given  as  2.3  km/sec 
within  very  restrictive  limits  of  uncertainty.  In  the  absence  of  con- 
fident pairing  of  other  model  and  actual  phases  nothing  can  be  said 
about  possible  variation  of  fault  rupture  velocity  along  the  fault. 

There  is  nothing  in  this  study  to  either  support  or  contradict  such  a 
possibility. 

Finally,  aside  from  the  overall  satisfactory  appearance  of  models 
of  the  last  series  employing  a source  rise  time  of  0.25  seconds  the 
width  of  the  particular  S-wave  pulse  to  which  we  have  been  referring 
also  supports  this  rise  time.  As  may  be  seen  by  examining  Figures 
4.3(34)  through  4.3(41)  this  particular  arrival  is  a sensitive  indica- 
tor of  model  source  rise  time  and  allows  us  to  contend  with  some  assur- 
ance that  0.25  seconds  is  the  best  to  be  found.  There  is  nothing  in 
the  character  of  the  vertical  S-wave  motion  coming  from  the  hypocentral 
region--and  thus  from  a completely  different  part  of  the  fault--to  con- 
tradict this  conclusion. 

In  terms  of  the  parameters  inferred  from  model  studies  to  be  most 
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appropriate  for  the  faulting  episode  associated  with  the  San  Fernando 
earthquake  of  February  9,  1971  we  may  conclude  that  rupture  occurred 
along  a surface  substantially  similar  to  that  of  Figure  4.3(30),  that 
the  form  of  the  dislocation  time  history  was  closely  approximated  by  a 
smoothed  ramp  whose  rise  time  was  near  0.25  seconds,  that  this  disloca- 
tion propagated  along  the  fault  surface  at  a velocity  near  2.3  km/sec, 
and  that  it  had  an  average  amplitude  of  between  one  and  two  meters. 
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Chapter  5 

Conclusions  and  Areas  for  Future  Study 
The  primary  goal  of  this  study  was  fairly  limited:  to  show  whe- 

ther or  not  a near-field  strong-motion  record  could  be  modeled  using 
the  mathematical  representation  of  motion  arising  from  dislocations  em- 
bedded in  an  elastic  homogeneous  isotropic  half-space  and,  thus,  to  add 
to  our  knowledge  of  the  seismic  source.  It  is  felt  that  this  goal  has 
been  realized;  if  not  perfectly,  then  at  least  substantially. 

To  gain  knowledge  of  real  phenomena  from  model  studies,  of  course, 
the  model  must  duplicate  the  actuality  in  some  way.  The  more  satisfac- 
tory the  duplication  the  more  far-reaching  and  likely  are  the  inferences 
that  may  be  drawn.  The  question  of  drawing  valid  inferences  from  an 
imperfect  model  is  a delicate  one  largely  defying  generalization  but  it 
is  believed  that  the  following  conclusions  may  be  reached  from  this 
study: 

1)  It  is  clear  both  from  the  discussion  in  Chapter  3 and  from 

I 

the  many  different  accelerograms  presented  in  Chapter  4 
that  in  order  to  adequately  characterize  the  high-order 
behavior  of  real  dislocation  source  functions  frequencies 

must  be  considered  that  are  high  relative  to  the  inverse 

P 

of  the  source  rise  time.  Given  most  proposed  rise  times 
and  the  instrument  response  spectra  of  modern  seismographs 

I 

this  implies  that  it  is  necessary  to  consider  accelerome- 
ter records.  This  is  far  from  a trivial  point  since  the 
accelerations  to  be  expected  from  different  seismic  sources 


178 


2)  The  characterization  of  accelerations  associated  with  any 
earthquake  as  being  the  result  of  an  essentially  random 
process  no  longer  seems  quite  as  necessary.  This  concept 
is  valid  only  so  long  as  it  is  maintained  that  small  por- 
tions of  the  source  act  independently  of  each  other.  The 
evidence  of  this  thesis  weakens  the  necessity  of  this  no- 
tion for  the  San  Fernando  earthquake.  The  general  charac- 
ter and  many  individual  features  of  an  actual  accelerogram 
may  evidently  be  duplicated  by  considering  a faulting  pro- 
cess that  is  perfectly  homogeneous  over  the  entire  fault 
surface.  It  is  not  necessary,  for  the  first  6.5  seconds 
of  the  Pacoima  Dam  record  at  least,  to  call  upon  special- 
ized situations  to  occur  at  particular  points.  This  is 
most  notably  true  for  those  portions  of  the  record  associ- 
ated with  the  peak  recorded  accelerations  in  that  part  of 
the  observed  accelerogram  for  which  a good  model  match  was 
obtained.  This  conclusion  suggests  many  interesting  rami- 
fications. For  example,  although  it  is  mechanically  simple 
to  postulate  situations  with  a dislocation  model  in  which 
long-  and  short-period  phenomena  are  virtually  unrelated, 
it  appears  that  this  is  not  always  the  case.  This,  in 
turn,  implies  that  the  vast  scatter  observed  in  empirical 
relations  between  magnitude  and  acceleration  are  not 
necessarily  a fundamental  property  of  the  seismic  source 
but  may  be  due  to  other  complexities  such  as  those  that 
undoubtedly  exist  in  the  medium  between  source  and  re- 
ceiver. This  argument  is  made  less  convincing  by  the 
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the  failure  of  any  model  to  match  the  observed  motion  in 
that  part  of  the  Pacoima  Dam  records  in  which  the  peak 
acceleration  of  1.25  actually  occurred.  It  is  also  less 
convincing  as  very  large-magnitude  events  are  considered 
where  some  material  inhomogeneity  must  occur  over  the 
fault  surface.  However,  the  implication  of  this  study  is 
that  aspects  of  fault  motion  of  substantially  disparate 
frequencies  are  more  closely  related  than  was  previously 
bel ieved. 

3)  A practically  important  principal,  termed  previously  the 
conservation  of  character,  has  emerged  fo^  strung-motion 
accelerometer  model  studies.  This  principal  states  that 
near  the  source  the  effects  of  variation  of  source  rise 
time,  source  rupture  propagation  velocity,  and  source- 
receiver  geometry  are  quite  separate  and  distinguishable 
for  high-frequency  ground  motions. 

Areas  for  future  study  are  many  but  they  can,  perhaps,  be  placed 
into  two  broad  categories.  The  first  involves  better  model  refinement 
so  that  better  knowledge  of  the  processes  occurring  at  the  seismic 
source  may  be  obtained.  For  example,  the  two-dimensional  source  may  be 
fundamentally  objectionable  in  that  it  requires  the  simultaneous  rupture 
of  greatly  separated  parts  of  the  fault.  This  is  not  intuitively  pleas- 
ing and  it  becomes  less  tractable  as  the  total  size  of  the  source  in- 
creases. It  is  clearly  desirable  to  extend  the  two-dimensional  method 
of  this  thesis  to  a three-dimensional  method.  Such  an  extension  will 
lead  to  a major  increase  in  computational  complexity  but  in  principal 
there  is  no  reason  why  it  cannot  be  accomplished.  The  groundwork  for 
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such  an  extension  has  already  been  done  (Johnson,  1974). 

The  second  major  area  for  future  research  is  the  converse  to  the 
first.  If  we  may  successfully  model  the  motion  on  a fault  from  a given 
accelerogram  then  we  may  also  successfully  produce  the  acceleration  to 
be  expected  from  a given  fault  motion.  The  implications  to  earthquake- 
resistant  design  of  structures  in  a given  geologic  setting  would  be 
significant.  Certainly  the  gross  characterization  of  the  strong  motion 
to  be  expected  at  a given  site  by  duration,  peak  acceleration,  general 
frequency  content,  and  amplitude  variation  with  time  would  be  an  impor- 
tant step  towards  rational  engineering  design. 

Finally,  as  with  any  study,  there  are  important  questions  left  to 
be  resolved  by  future  investigations..  The  failure  to  duplicate  the  ob- 
served motion  after  the  first  6.5  seconds  is  disappointing  particularly 
in  view  of  the  ■'"act  that  the  peak  acceleration  occurred  in  the  later 
part  of  the  Pacoima  record.  The  failure  of  any  model  horizontal  com- 
ponent to  equal  the  success  of  its  vertical  counterpart  is  a fundamental 
problem  that  remains  unexplained.  As  more  data  becomes  available  from 
future  earthquakes  it  is  hoped  that  these  difficulties  may  be  overcome. 

Perhaps  it  is  best  merely  to  contend  that  many  of  the  properties 
of  strong  motion  that  are  now  obscure  are  indeed  fathomable.  In  a small 
way  it  is  believed  that  this  work  supports  this  contention. 
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